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Abstract 
(N 

^ We consider branching Brownian motion on the real line with the following selection mcch- 

£L| anism: Every time the number of particles exceeds a (large) given number N, only the A^ 

•^ right-most particles are kept and the others killed. After rescaling time by log 3 N, we show that 

_i the properly recentred position of the [ct./Vj-th particle from the left, a e (0,1), converges in 

law to an explicitly given spectrally positive Levy process. This behaviour has been predicted 

to hold for a large class of models falling into the universality class of the FKPP equation with 

weak multiplicative noise [Brunet et al., Phys. Rev. E 73(5), 056126 (2006)] and is proven here 

Ph for the first time for such a model. 
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1 Introduction 

1.1 Definition of the model and statement of the main result 

In this article, we consider an instance of branching Brownian motion (BBM) with selection, dubbed 
the A-BBM and defined as follows: Given a reproduction law {q{k))k^o with m = ^(k — l)q{k) > 
and finite second moment, particles diffuse according to standard Brownian motion and branch at 
rate q{k)/{2m) into k particles, for every k ^ 0. We fix a (large) parameter TV" and select particles 
according to the following simple mechanism: Each time the number of particles exceeds A, we keep 
only the A right-most and instantaneously kill the others. In Section [l.2[ we explain the motivation 
behind this system. We first state our result. 

For a finite counting measure v on R, a e (0, 1) and A e N, we define 

qu^(i^) = inf{x e R : u([x, oo)) < aN}. 



For JVeN large enough, we then define a^ = log N + 3 log log N and 



i 7T 2 =1 7T 2 | 37r 2 loglogiV | o 

^ Y a^ 2 log 2 TV log 3 iV °Mog 3 iV 

We then set M~f(t) = qu^(u^) — fiNt, where vf is the counting measure formed by the positions 
of the particles of iV-BBM at time t. Our main theorem is the following: 

Theorem 1.1. Suppose that at time there are N particles distributed independently according to 
the density proportional to sin(Trx/aN)e~ x l( . ajv )(x) . Let a 6 (0, 1) and define x a by J x ye~ y dy = 
a. Then the finite- dimensional distributions of the process [M^ {tlog A")) converge weakly as 
N — ► oo to those of the Levy process (L t + x a ) t ^o with Lq = and 

rCO 

log E[e iXLl ] = i\c + 7T 2 \ e iXx -l-i\xl (x<1) A(dx). (1.1) 

Jo 

Here, A is the image of the measure (x~ 2 lr x> Q\)dx by the map x *-*■ log(l + x) and c £ R is a 
constant depending on the reproduction law q(k). 

1.2 Motivation and related work 

The iV-BBM and the results of this article are related to several other mathematical models, a 
few of which we wish to outline in this section. We are not going to elaborate on the applications 
to other sciences such as physics or biology, for which we refer to |17t [T9] and the review articles 
[561E2]. 

Most importantly, the A r -BBM is a prototype for noisy traveling waves of FKPP type, i.e. 
travelling waves of the stochastic partial differential equation dtu = d 2 u + u{l — u) + ^Ju/N W, where 
W is space-time white noise. This equation is known as the FKPP equation with weak multiplicative 
noise and is believed to share the same phenomenology as the A r -BBM. The connection between 
BBM without selection and the FKPP equation without noise is indeed well known, at least since 
McKean's work [48J, and has had many fruitful applications, one of the most important being 
Bramson's study of the law of the right-most particle of BBM [14] . The relation between the Af- 
BBM and the noisy FKPP equation is less explicit, but there is an exact duality relation between this 
equation and a system of branching and coalescing Brownian motions |55j . More generally, particle 
systems governed by the FKPP equation in the large population limit can often be modeled by the 
noisy FKPP equation, the parameter A^ having the meaning of an "effective population size" ([16 , 
[561 Chapter 7], [52]). 

Brunet, Derrida, Mueller and Munier [IT] give a semi-deterministic description of such systems 
and basically obtain our Theorem |1.1| Their arguments, which we outline in Section 1.3, work for 



very general models but rely on some fundamental assumptions which are intuitively convincing but 
far from easy to prove in particular settings. The author is indeed aware of only three mathematically 
rigorous articles in this area: Berard and Gouerc [6j prove the log - N correction to the linear speed 
of the A r -particle branching random walk (Af-BRW). Durrett and Remenik [27] study the empirical 
distribution of A^-BRW and show that it converges to a system of integro-differential equations with 
moving boundary. Finally, Mueller, Mytnik and Quastel [49] obtain an estimate for the speed of 
the traveling waves of the noisy FKPP equation, which partly confirms the physicists' predictions. 
To our knowledge, the present article is therefore the first example of a rigorous proof of the full 
statistics for a particular system. 

In order to prove Theorem |l.l[ we approximate the A^-BBM by BBM with absorption at a linear 
space-time barrier. This process is well-studied [261 E01 [32], [39] E] and is much more tractable 
than A r -BBM due to the greater independence between the particles and its connection with some 



differential equations |32] Rf)] 151] . Of paramount importance to us is the article by Berestycki, 
Berestycki and Schweinsberg, who study BBM with absorption at a barrier with slope fijy, starting 
from TV particles. If Nt denotes the number of particles at time t in their system, they show that the 
process (N t )t^o converges in law to Neveu's continuous-state branching process and the genealogy of 
the particles to the Bolthausen-Sznitman coalescent, which has been observed before [18]. Although 



we do not apply their results directly, our proof of Theorem 1.1 uses many of their ideas. We remark 



that the genealogy of the iV-BBM should converge to the Bolthausen-Sznitman coalescent as well, 



but our methods do not seem to be applicable (see the end of Section 1.3 for details). 

We finally mention the relation between the iV-BBM or iV-BRW and an instance of a Fleming- 
Viot-type process: Let N particles perform independent continuous-time nearest-neighbour random 
walks on the integers with drift towards the origin. Furthermore, as soon as a particle reaches the 
origin, let it jump at random onto one of the other particles. This type of model was introduced in 
order to yield a particle picture of quasi-stationary distributions [21], and it is indeed conjectured 
that the empirical distribution of this system under the stationary distribution converges, as N goes 
to infinity, to the quasi-stationary distribution of the random walk (see [2] |3] for recent progress 
on this subject). The relation with the iV-BBM is obvious and we conjecture that the empirical 
distribution of a general iV-BRW is an ergodic process whose stationary distribution converges to 
a deterministic measure, which is the quasi-stationary distribution of the underlying random walk 
with a certain drift towards the origin. 

1.3 Overview of the proof 

Since this article is fairly long, we give a detailed overview of our main results. We recall the 
heuristic semi-deterministic description of iV-BBM established in [17j : 

1. Most of the time, the particles are in a meta-stable state. In this state, the diameter of 
the cloud of particles (also called the front) is approximately logiV, the empirical density 
of the particles proportional to e~ x sin(7rx/log N), and the system moves at a linear speed 
^cutoff = 1 — tt 2 /(2 log 2 N). This is the description provided by the cutoff approximation from 
11511 mentioned above. 



2. This meta-stable state is perturbed from time to time by particles moving far to the right and 
thus spawning a large number of descendants, causing a shift of the front to the right after a 
relaxation time which is of the order of log N. To make this precise, we fix a point in the bulk, 
for example the barycentre of the cloud of particles, and shift our coordinate system such that 
this point becomes its origin. Playing with the initial conditions of the FKPP equation with 
cutoff, the authors of [IT] find that a particle moving up to the point log N + x causes a shift 
of the front by A = log I 1 + f J , for some constant C > 0. In particular, in order to have 
an effect on the position of the front, a particle has to reach a point near log N + 3 log log N. 

3. Assuming that such an event where a particle "escapes" to the point log N + x happens with 
rate Ce~ x , one sees that the time it takes for a particle to come close to logiV + 3 log log N 
(and thus causing shifts of the front) is of the order of log N » log N. 

4. With this information, the speed of the front is found to satisfy vn — v cuto s * 7r 2 , og °fj and 
the fluctuations are found to be given by the Levy process from Theorem 1 1.1] 

Berestycki, Berestycki and Schweinsberg [7J made a first attempt to put this description onto a 
rigorous foundation. They consider BBM with absorption at the origin and with drift —fJ-N- Their 
starting point is to introduce a second barrier at the point <ztv,a = a N — A for some large positive 
constant A and divide the particles at time t into two parts; on the one hand those that have stayed 
inside the interval (0, ajv,^), on the other hand those that have hit the point cln,a before hitting 0. 



This corresponds roughly to the division of the process into a deterministic and a stochastic part. 
Indeed, killing the particles at cln,a prevents the number of particles from growing fast and thus 
permits to calculate expectations and variances of various quantities. For example, if at time we 
have iV particles distributed according to the meta-stable density, then the variance of the number 
of particles at the time log N is of order of e~ A N 2 . For large A, the particles inside the interval 
(0, ojv,a) therefore behave almost deterministically at the timescale log N. Moreover, the leading 
term in the Fourier expansion of the transition density of Brownian motion (with drift — /U/v an d 
killed at the border of the interval (0, ajv,^)) is proportional to e~^ NX sin(7r x/cln,a), which explains 
the meta-stable density predicted by the physicists. 

As for the particles that hit a/v,Aj the authors of [7j find that 1) the number of descendants at 
a later time of such a particle is of the order of e~ NW, where W is a random variable with tail 
P(W > x) ~ 1/x, as x —> oo and 2) the rate at which particles hit the right barrier is of the order of 
e / log N. Putting the pieces together, they then show that the process which counts the number 
of particles of the system converges in the log N timescale to Neveu's continuous-state branching 
procesqjand its genealogy to the Bolthausen-Sznitman coalescent. 

Our proof of Theorem |1.1| builds upon the ideas of [7J presented above. The basic idea is to 
approximate the iV-BBM by a BBM with absorption at a random barrier, which is chosen in such 
a way that it keeps the number of particles almost constant. We call the resulting system the "B- 
BBM" (B stands for "barrier"), see Figure [I] for a graphical description. The B-BBM takes two 
(large) parameters a and A, which have similar purposes than aw,A an d A above. We then set 
\i = y/l — TT 2 /a 2 and start with BBM with drift —\x and an absorbing barrier at the origin. Now, at 
the beginning, this barrier stays at the origin and does not move. When and only when a particle 
hits a and spawns a lot of descendants do we increase the drift to the left. This increase is in order to 
kill particles and thus make the population size stay almost constant. Note that moving the barrier 
to the right is an equivalent operation, but increasing the drift is technically more convenient. After 
the system has relaxed (which takes a time of order a 2 ) the drift is set to —/j, again and the process 
is repeated. 

As in [7], an important quantity is Zt = ^w z(X u (t)) , where wz(x) = ae^ x ~ a ' sin™. Here, we 
sum over all the particles u alive at time t and X u (t) denotes the position of the particle u at time 
t. The process (Zt)t^o is important because it is a martingale for BBM with absorption at and 
a and drift —\x. Furthermore, it gives the approximate number of particles at a time t » a 2 . More 
precisely, set N = \2ire a _3 e /ia ] and suppose we kill particles at and a. The expected number of 
particles at a time t, with a 2 « t « a 3 , is then approximately Ne~ Zq. Moreover, the variance is of 
the order of N 2 e~ 2 Zq. Therefore, if Zq *s e then the number of particles is concentrated around 
its expectation for large A. 

When a particle hits the right barrier at the time r, say, we absorb its descendants at the space- 
time line Jz? T : x = a — y + (1 — fj){t — r), where y is a large constant depending on A only (this idea 
comes from |7J). In doing so, the number of particles absorbed at the barrier has the same law as in 
BBM with drift —1 and absorption at —y, starting from a single particle at the origin. Moreover, 
the time it takes for all the particles to be absorbed depends only on y, not on a. Now, if x±, . . . , x n 
are the positions of the absorbed particles and Z' = ^iWz(xi), then the number of descendants 
of this particle at a later time is of the order of e~ A NZ', provided that the drift stays constant. 
Consequently, we say that a breakout occurs, whenever Z' > ee , where e = A~ v for some large p 
to be chosen later. 

A continuous-state branching process (CSBP) (Z t )t^o is a time-changed Levy process without negative jumps: 
at time t, time is sped up by the factor Z t -. CSBPs are scaling limits of Galton- Watson processes and thus have an 
inherent notion of genealogy. Neveu's CSBP is the CSBP with Levy measure a; _2 l( a;> o) dx, whose genealogy is given 
by the Bolthausen-Sznitman coalescent [5j. As a wise reader, you have read the introductory chapter, such that you 
know what this coalescent is. 



In order to define a breakout properly, we classify the particles into tiers. Particles that have 
never hit the point a form the particles of tier 0. As soon as a particle hits a (at the time r, say) it 
advances to tier 1. Its descendants then belong to tier 1 as well, but whenever a descendant hits a 
and has an ancestor which has hit the line Jz? T after r, it advances to tier 2 and so on. Whenever a 
particle advances to the next tier, it has a chance to break out. We can then define the time T of 
the first breakout and will indeed show that T is approximately exponentially distributed with rate 
proportional to £ -1 a~ 3 . Interestingly, we will see that with high probability breakouts only occur 
from particles which are of tier or 1. In fact, the number of breakouts occuring from particles of 
tier 1 between the times and a 3 is approximately proportional to A (and the remaining « e^ 1 
breakouts occur from particles of tier 0). 




Figure 1: A caricatural graphical description of the B-BBM. The fugitive and its descendants are drawn with 
thick lines, the other particles with thin lines. A breakout happens at time T and the barrier is moved from 
the time T + on. The process starts afresh at the time O x . Note that technically we increase the drift to the 



left instead of moving the barrier. The three important timescales (1, 
that a sk log N. 



and a ) are shown as well. Recall 



After the breakout, we will then increase the drift to the left slightly, in order to kill more 
particles than usual (remember that increasing the drift to the left corresponds to moving the 
barrier to the right). For this, we first choose a family of increasing smooth functions (f x )x^o with 
/x(0) = and / x (+co) = x for all x ^ 0. Such a function will be called a barrier function. The 
drift after a breakout is then set to fit = fi + (d/dt)f/\(t/a 2 ), where A is the total amount by which 
we have to move the barrier. Thus, the only randomness in the choice of the barrier is in its total 
shift, not in its shape. Looking at the definition of Zt and the fact that Zq % e , one easily guesses 
that we have to choose A = log(l + e Z 7 ) in order to get Zt ultimately back to its initial value. 



This already explains the convergence of the barrier to the Levy process given by ( 1.1 ): On the one 
hand, we have Z' * ttW, where W is the random variable mentioned above. This implies that the 
law of e~ Z' conditioned on Z' > ee is approximately ex~ 2 lt x y> £ \ dx for large A and ar\ On the 
other hand, we will show that breakouts occur at a rate proportional to £ _1 a~ 3 . Together with the 



definition of A, this explains the Levy measure A(dx) in (1.1). 



The statement "for large A and a" means that we let first a, then A go to infinity, see Section 6.1 



As for the shape of the barrier, it is determined by the fact that we want the number of particles 
at each time to be approximately N. This is accomplished by choosing 

fA(t) = log (l + (e A - l)vr- 2 e- 2i / 2 ^(l,t) 



where the theta function 9(x,t) is defined in (2.1). 



We remark that in order to study the B-BBM up to the time T of the first breakout, we need 



to study it conditioned to break out at time t for every t ^ 0. We will see in Section 6.4 that this 
will lead to a decomposition of the particles into 1) a fugitive, which is conditioned to break out 
at t and which will effectively be a spine, and 2) the other particles conditioned not to break out 
before t. This is essentially a Doob transform of the process. Furthermore, we will see that the tier 
1 particles will have an essential role: At the timescale a 3 they will lead to an additional shift of 
the barrier by an amount of the order of A. This term will play the role of the linear compensation 
that is necessary in order to obtain in the limit the Levy process of infinite variation stated in 
Theorem 11.11 

We now describe how we use the results on the B-BBM in order to prove Theorem |1.1[ Initially, 
our plan was to couple the iV-BBM and the B-BBM, i.e. construct them on the same probability 
space. We would then assign a colour to each particle: blue to the particles which appear in the 
iV-BBM but not in the B-BBM, red to those that appear in the B-BBM but not in the A-BBM, and 
white to the particles that appear in both processes. Our aim was then to show that the number of 
blue and red particles was negligible after a time of order a 3 . This, unfortunately, did not work out, 
because we were not able to handle the intricate dependence between the red and blue particles. 

Instead, we couple the A-BBM with two different processes, the B - and the B"-BBM, which 
are variants of the B-BBM and which bound the position of the A-BBM in a certain sense from 
below and above, respectively. The B-BBM is defined as follows: Initially, all particles are coloured 
white and evolve as in the B-BBM. A white particle is coloured red as soon as it has N or more 
white particles to its right. Children inherit the colour of their parent. After a breakout and the 
subsequent relaxation, all the red particles are killed immediately and the process restarts with 
the remaining particles. It is intuitive that the collection of white particles then bounds the N- 
BBM from below (in some sense) because we kill "more" particles than in the A-BBM. Indeed, in 



Section 10.1, we show by a coupling method that the empirical measure of the white particles in 
B-BBM is stochastically dominated by the one of the A-BBM with respect to the usual stochastic 
ordering of measures. It then remains to show that the number of red particles in B-BBM is 
negligible when A and a are large. We do this through precise estimates on the number of particles 
in the interval [r, oo) for every r ^ 0. These allow us to estimate the expected number of particles 
which turn red at the point r. It turns out that this expectation is small enough, which permits to 
conclude. 

The definition of the B'-BBM, which is used to bound the ./V-BBM from above, is more intricate 
than the one of B-BBM. Again, we colour all initial particles white and particles evolve as in B- 
BBM with the following change: Whenever a white particle hits and has less than N particles to 
its right, instead of killing it immediately, we colour it blue and let it survive for a time of order a 2 . 
More precisely, we cut time into intervals I n = [t n ,t n+ i), with t n = Kna 2 for some large constant 
K. A particle which gets coloured blue during I n then survives until the time t n+ 2- At this time, 
all of its descendants to the left of the origin are killed and the others survive. It will turn out 
that this system bounds the A-BBM from above with high probability and that the number of blue 
particles will remain negligible during a time of order a 3 , as long as A and a are large. 

We note that although our technique of bounding the A-BBM from below and from above works 
well for the position of the particles, it does not give us information about the genealogy; the reason 
being that the coupling deforms the genealogical tree of the process. Thus, although it should not 
be difficult to show that the genealogy of the B-BBM (and of the B - and B'-BBM) converges to the 

7 



Bolthausen-Sznitman coalescent we do not know at present how one could transfer this information 
to the iV-BBM. 



1.4 Notation guide 

This article is quite long and uses a wealth of different notation. Below is a list of recurrent symbols, 
roughly in the order of their first appearance. Following this list are some further remarks about 
notational conventions. 



Symbol 

q(k) 
m 

/3o_ 
e,e 

w x 

Pt(x,y) 

Et 

I a {x,S),J a (x,S) 

w x ^ w x i ,y 

rr taboo' vy taboo 
N, N* 

u 



sui Oil , tt u 

X u (t) 

px W pi 



E p 



&se 

L 

m 2 

W 

Pt(x,y) 

P /' E / 
wz(x) 

wy{x) 

Zt,Y t 

N t 

Rt 

P,E 

A,a,e,ri,y,( 



A 



(0 



ce(u,t) 

7 {l) 



Meaning 

Reproduction law 

™ = Yi k ( k ~ l )l( k ) 
Branching rate, /3o = l/(2m) 

Theta functions 

Law of Brownian motion started at x 
Transition density of Brownian motion killed outside [0, a] 
Error term 

Integrals related to Brownian motion killed outside [0, a] 
Law of Brownian taboo process and its bridge 
The set of natural numbers including and excluding 0, resp. 
The space of individuals 
A realisation of a branching Markov process 
Lifetime, birth and death times of an individual u 
Set of individuals alive at time t 
Position of the individual u at the time t 
Law of and expectation w.r.t. a branching Markov process 
(later: BBM) started at x or with particles distributed ac- 
cording to a finite counting measure v 
cr-algebra with information up to time t 
cr-algebra with information up to the stopping line Jz? 
Stopping line generated by a stopping time T 
A random variable with law q{k) 
m2 = E[L(L-l)] 
Seneta-Heyde martingale limit of BBM with absorption 

/i = Vl-vr 2 /« 2 

Density of BBM killed outside an interval 

Law of and expectation w.r.t. BBM with varying drift 

w z {x) = asin(vrx/a)e^ :c " a )l( :re[0 . a]) 

wy(x) = e^- a h {x>0) 

Sums of wz(x),wy(x) over the positions of time t particles 

Number of particles at time t 

Number of particles hitting a up to time t 

Law/expectation of BBM weakly conditioned not to hit a 

Parameters of BBM before a breakout and B-BBM 

Stopping line of tier I particles at time t 

Stopping line of tier I particles hitting the critical line before 

t 

Stopping line of tier I particles hitting a before t 

Stopping line of descendants of (u, t) hitting the critical line 



Sum of wz (x) over particles from ,jV t 



ID 



Sect. 



2.1 



2.1 



2.1 



2.1 



2.1 



2.2 



3.1 



3.1 



3.1 



3.1 



3.1 



3.1 



3.1 



3.2 



3.2 



3.2 



5.1 



5.1 



5.1 



5.2 



5.2 



5.2 



5.3 



5.4 



5.5 



6.1 



6.1 



6.1 



6.1 



6.1 



6.1 



Z t ,Y t ,R t ,N t 

PB 

Q a 

TO) 

P,E(alsoP A ,E A ) 
7 {i) 

A,A t 

Zt, Z t 
Z t ,Y t ,N tl ... 



z t fug ,y/ ug ,7v/ ug ,. 



T,T n 

rp — rp — rp + rp + 



G n 

7o 

G<?/, Gf ug , G, 

G, G/\, G n bab 



Zt = Z\ (from Section 6 on). Same for Yt, Rt and Nt 
Probability of a breakout 
P a conditioned not to break out 
Time of first breakout from tier I 

Law/expectation of BBM conditioned not to break out before 
some fixed time 

Sum of wz(x) over particles from S^ 
The fugitive (the particle which breaks out) 
The descendants of the fugitive at time t (split into two 
groups) 

Value of Zt restricted to particles from A and A, resp. 
Value of Zt, Yt, Nt . . . restricted to particles not related to the 
fugitive 

Value of Zt,Yt,Nt ■ ■ ■ restricted to particles descending from 
the fugitive 

Time of the first/n-th breakout 
Random times related to the times T, T n 
The total shift of the barrier after a breakout 
The position of the first n pieces of the barrier at time t (bar- 
rier process) 

Beginning of the n + 1-th piece of B-BBM (i.e. time of the 
n-th breakout plus relaxation time) 
"Good event" related to the first n pieces of B-BBM 
7o = (vrp B e A ) _1 . 

Several good sets related to a piece of B-BBM 



6.1 



6.1 



6.1 



6.1 



6.3 



6.3 



6.4 



6.4 



6.4 



1A 
71 



7T 
7A 



7.1 



7.1 



7.1 



7.2 



The complement of an event G is denoted by G c . 

From Section [5] on, the symbol C stands for a positive constant, which may only depend on the 
reproduction law q and the value of which may change from line to line. Furthermore, if X is any 
mathematical expression, then the symbol 0{X) stands for a possibly random term whose absolute 
value is bounded by C|X|. 

In Section [6j we introduce two parameters A and a and will first let a then A go to infinity. 
This will be expressed by the statements "for large A and a we have. . . " or "as A and a go to 



infinity. . . " (see Section 6.1 for a precise definition). These phrases will become so common that in 



Sections [7] to [9] they will often be used implicitly, although they will always be explicitly stated in 
the theorems, propositions, lemmas etc. Section k3] furthermore introduces the notation o(l), which 
stands for a (non-random) term that only depends on the reproduction law q and the parameters 
A, a, e, rj, y and £ and which goes to as A and a go to infinity. 

Sections [8] and [9] each use special notation which only appears in those sections. This notation 
is defined at the beginning of both sections. Moreover, in both sections, we sometimes denote 
quantities which refer to descendants of the fugitive after a breakout by the superscript "fug" . 



2 Brownian motion in an interval 



In this section, we recall some explicit formulae concerning real- valued Brownian motion killed upon 
exiting or conditioned not to exit an interval. 



2.1 Brownian motion killed upon exiting an interval 

We define for x e R and t > the following function of Jacobi theta-type: 



;, t) — 2 + 2j n =l e 



oo —-Frnt 



cos(imx) 

(x-2n) 2 \ 
2t /' 



(2.1) 
(2.2) 



the two representations being related by the Poisson summation formula (see [5j, §9). We also define 



0(i) 



2 . 2 + d 



7T 



2 e2 'dt^ 1 '* ) ' 



which is a smooth function on R + . By (2.1) and (2.2), one can show that is stricly increasing^ 
with 0(0) = and 0(+oo) = 1. 

Various quantities of Brownian motion killed upon exiting an interval can be expressed by 6(x, t). 
For x e R, let W x be the law of Brownian motion started at x, let (Xt)t^o be the canonical process 
and let H y = inf{£ ^ : X t = y}. For a > and x e [0, a], denote by Wy lled a the law of Brownian 
motion started at x and killed upon leaving the interval (0, a). Let p%(x, y) be its transition density, 
i.e. 

Pt(x, y)dy = W^ lcdja (X t e dy) = W x {X t e dy, H AH a >t), x,ye [0, o]. 

Then (see [36], Problem 1.7.8 or [13J, formula 1.1.15.8), 



(2.3) 



Pt(x,y) = a l 



x — y t 
a a z 



x + y t 



a cf 



(2.4) 



Equation (2.1) then yields 



p t {x 



■v)-;2 



e 2a 2 



sinf-zm-) sin(vrn-). 

v a' v a' 



(2.5) 



ra=l 



This representation is particularly useful for large t: Define 

oo 

E t = vr 2 £ ^e-^ 2 ^ 2 - 1 )*/ 2 . 

n=2 



(2.6) 



By (2.5) and the inequality | sinn:r| =% nsinx, x e [0,7r], one sees that 



J2 9 

TV , £ 

e^ p t (x, y) sin(7rx/a) sin(7ry/a) 



sc 



E t i a 2 - sin(7rx/a) sin(7ry/a). 



The Green function (see e.g. |38j, Lemma 20.10, p379) is given by 



f 

Jo 



Pt{x,y)dt = 2a (x Ay)(a-x v y). 



(2.7) 



(2. 



Set H = H AH a and define rf(x) = W X (H e dfc, X H = a)/dt. Then (see [13], formula 1.3.0.6), 

r a t {x) = V (- - 1 * V (2.9) 

where 0' denotes the derivative of with respect to x. 



3 More precisely, by elementary computations, (2.1 1 gives to e R-+, such that is strictly increasing on (to, oo) and 
(|2.2| gives ti > to, such that # is strictly increasing on [0, ti). 
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The following two integrals are going to appear several times throughout the article, which is 
why we give some useful estimates here. For a measurable subset ScR, define 



I a (x,S) = W*(e£ H "l (Ho>HaeS) ) = J" 

1 (x,y,S)= f 

Js 



e2^ s r'i(x)ds, 



and 



5n(0,oo) 



)Sn(0,co) 

which satisfy the scaling relations 



e^ s p a s {x,y)ds, 



' x S 



i a ( x ,s)=i(-,^), r( x , y ,s) = aj(-,y,^) 



(2.10) 
(2.11) 

(2.12) 



with I = I 1 and J = J . The following lemma provides estimates on I(x,S) and J(x,y,S). It is 
easily proven from the above equations (see \i7\ Lemma 2.2.1]). 

Lemma 2.1. There exists a universal constant C, such that for every x e [0, 1] and every measurable 
S cz R + , we have 



\I(x, S) — tt\(S) sin(7rx)| < C(x a E- m fs(l a X(S)) sin(7rx) J, and 
\J(x,y,S) -2A(5')sin(7rx)sin(7ry)| < C{ [{x a y)(l - (x v y))] a E- mi ssin(7rx) sin(7ry) ) , 



where X(S) denotes the Lebesgue measure of S and -EW^ is defined in (2.6) 



2.2 The Brownian taboo process 

The Brownian taboo process on the interval (0, a) is the diffusion with infinitesimal generator 

\ (dx) + a co * ^ dx- ^he name °f this process was coined by F. Knight [40J who showed that 
it can be interpreted as Brownian motion conditioned to stay inside the interval (0, a). It is the 
Doob transform of Brownian motion killed at and a, with respect to the space-time harmonic 
function h(x,t) = sin(7rx/a) exp(-7r 2 t/(2a 2 )). As a consequence, its transition density is given for 
x,ye [0, a] by 

sin(7n//a) 



taboo(O.a)/ \ 

Pt {x,y) 



■t 2 2 

-e^ f p a t {x, y) = - sinVy/a)(l + 0(l)E t/a2 ), 
sm{Trx/a) a ' 



(2.13) 



by (2.7). For x e [0,o] we denote the law of the Brownian taboo process on (0,a) started from 



x by Wjabooa- Often we will drop the a if its value is clear from the context. We also denote by 
^ta'boo t ne la w of taboo bridge from x to y of length t. Note that the taboo process is self-dual 
in the sense that for a measurable functional F and t > 0, we have W^^ [F((X S ; < s < t))] = 

The following lemma will be needed in Sections [6] and [7j 

Lemma 2.2. Let c > 0. There exists a constant C ', depending only on c, such that we have for 
every x, y e [0, a], 



w? 



taboo 



r 

Jo 



and for t ^ a 2 , 
Ift^a 2 , we still have, 



W x '!' y 

taboo 



f 

Jo 



e- cXa ds 



W x,ty 
taboo 



< C(t/a 3 + (1 a x- 1 )), 
^c(t/a 3 + (lAx- 1 ) + (lAy- 1 )). 

f 

Jo 



e- cXfl ds 



^C. 



(2.14) 
(2.15) 
(2.16) 



11 



Proof. We write k{x) = e cx . We first show that (2.14) implies (2.15). By the self-duality of the 
taboo process, we have 



W x,ty 
taboo 



f k(X s ) ds 
Jo 

It therefore remains to prove that 



w x,t,y 

taboo 



rt/2 

k{X s 
JO 



ds 



w y,t,x 

taboo 



rt/2 

k(X s ) 
Jo 



ds 



r x,t,y 



rt/2 

k(X s 
Jo 



Conditioning on 0~(X S ; < t ^ t/2), this integral equals 

-p\^ 00 {X t/2 ,y) rt/2 



)ds\ ^C(t/a 3 + (l ax" 1 )). 



E{x, y) = W t x ahoo 



pf hoo (x,y) 



f 

Jo 



k{X s ) ds 



By (2.13), there exists a universal constant C, such that for t > a 



E(x,y)^CW t x abc 



rt/2 

k(X s ) ds 
Jo 



Equation (2.14) therefore implies (2.15). 



Heuristically, one can estimate the left side of (2.14) in the following way: Since k(x) is decreasing 



very fast, only the times at which X s is of order 1 contribute to the integral. When started from the 
stationary distribution, the process takes a time of order a 3 to reach a point at distance 0(1) from 
the origin |43j and it stays there for a time of order 1, hence the integral is of order t/a 3 . When 
started from the point x, an additional error is added, which is of order 1, when x is at distance of 
order 1 away from 0. Adding both terms gives the bound appearing in the statement of the lemma. 
The actual calculations can be performed in the following way: Let Y be a random variable 
with values in (0, a) distributed according to fh(dx) := 2/asin 2 (-7rx/a) dx, which is the stationary 
probability measure of the taboo process. Let Hy = inf{£ > : X s = Y}. We then have 



wz 



taboo 



Jo 



ds 



n; 



taboo 



f Y k(X s )ds+ f k(X s )ds 

Jo JH Y 



^W t x 



taboo 



rH Y 

k(X s 
Jo 



ds 



" taboo 



f k(X s 

Jo 



ds 



Both of these quantities can now be expressed through the scale function and speed measure of the 



taboo process and after some calculations (see |47l Lemma 2.2.2] for details), one obtains (2.14). 



In order to prove (2.16), a different method is needed. We may assume that x, y < a/2, otherwise 



we decompose the path at the first and/or last time it hits a/2 and bound the parts above a/2 
trivially by a 2 e~ ca ' 2 = 0(1). The transition density of the taboo bridge can be written 



W t 



x,t,y 



taboo 



I, e dz 



p a s (x,z)pt_ s (z,y) 
Pt(x,y) 



dz. 



If we denote by p® (x, y) = (2irt) 1 ' 2 exp( 



z z + x' 



)/2t)2sinh(zx/t) the transition density of Brow- 



nian motion killed at 0, then we have the trivial inequality p%(x,y) < P®{x,y) and furthermore 
pf(x, y) > Cpt(x, y) for x, y < a/2 and t =% a 2 by Brownian scaling. It follows that 



w x,tv 

taboo 



i 

Jo 



k(X s ) ds 



< CR x ^ y 



f k{X s 
Jo 



ds 
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where R x > l >y denotes the law of the Bessel bridge of dimension 3. This Bessel bridge is the Doob 
transform of the Bessel process started at x with respect to the space-time harmonic function 
h y (z, s) = Pt_ s (z, y)/Pt(x, y). By the standard theory of Doob transforms, this is the Bessel process 
with additional drift 



d z d zy z y zy 

-{\oghy{z,s)) = h — logsmh- = — \- coth 



dz '' ' t — s dz t — s t — s t — s t — s 

This in an increasing function in y and standard comparison theorems for diffusions (see e.g. [53j, 

Theorem IX. 3. 7) now yield that for y\ ^ y%, we have R x,t,y2 \k(X s )\ ^ R x ' t,yi [k(X s )], since A; is a 

decreasing function. This is true in particular for y\ = 0. Using the self-duality of the Bessel bridge, 

we can repeat the same reasoning with x. We thus have altogether 

r rt -, r rt 

W x,t,y 



taboo 



f k(X s )ds] s; CR°' tfi \ f k(X s )ds 
Jo - L J 



for any x,y ^ a/2. This calculation can be done explicitly and yields (2.16). □ 



3 Preliminaries on branching Markov processes 

In this section we recall the definition and basic properties of branching Brownian motion and 
branching Markov processes in general. 

3.1 Definition and notation 

A branching Markov process in a Polish space $ can formally be defined in many ways, e.g. as 
a Markov process in the space [J n S n [35] or as a measure- valued Markov process |25j . A more 
powerful definition which maintains the genealogy, uses Neveu's marked trees [501 [231 [22] . We will 
briefly describe this construction and refer to these references or to |47[ Section 3] for details. The 
basic idea is to define a tree t as a suitable subset of U = {0} u U n>1 N* n , the space of individuals 
(we use the notation N* = {1,2,3, . . .} and N = {0} u N*). The space U is endowed with the 
ordering relations < and < defined by 

u < v <^^> 3w e U : v = uw and u < v <^^> u < v and u # v. 

To each individual u in a tree t, we then associate a trajectory E u 6 D{S). Here, D(S') is the 
Skorokhod space of functions S : [0, oo) — * $ u {A} (A ^ $ being a cemetary symbol) which are 
right-continuous with left limits, with S(0) ^ A and for which S(t) = A implies 2(s) = A for all 
s > t. For a trajectory 3 e D(S'), we define £(H) = inf{£ 5= : 3(£) = A} and set Q u = C(H U ) for 
an individual u 6 U. The birth and death times of the individual u are the defined by b u = Yjv<u Cv 
and d u = b u + Q u . The set of individuals alive at time t (i.e. for which b u < t < d u ), is denoted by 
N(t). The position of an individual u at time t is defined by X u (t) = H„(£ — b v ), where v 6 U is 
such that v e JV(£) and v < u. If d u ^ t, then we set X u {t) = A. 

Now suppose we are given a defective strong Markov process X = (Xt)t^o on $, with paths 
in D((f>). The law of X started in x e $ will be denoted by P . For simplicity, we will assume 
that for every x s £, we have CPO < °°> P -almost surely. Furthermore, let ((q(x, k))keN)xeS be a 
family of probability measures on N, measurable with respect to x. Then we define the branching 
Markov process with particle motion X and reproduction law q as the (unique) family of probability 
measures (P x ) xe g on £1 which satisfies 

P x (du)=P X (dX )q(X (( -),k )Y[P X 0^-\dJ^. (3.1) 

Note that by looking at the space-time process (Xt, t)t^o, we can (and will) extend this definition 
to the time-inhomogeneous case. 
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3.2 Stopping lines 

The analogue to stopping times for branching Markov processes are (optional) stopping lines, for 
which several definitions exist. For branching Brownian motion, they have first been defined by 
Chauvin |22j , the definition there is however too restricted for our purposes. Jagers [37] has given a 
definition of more general stopping lines for discrete-time branching processes; our definition will in 
fact mix both approaches. Note also that Biggins and Kyprianou |9j build up on Jagers' definition 
of stopping lines and define the subclasses of simple and very simple stopping lines (again for 
discrete-time processes). Chauvin's definition then corresponds to the class of very simple stopping 
lines. 

We first define a (random) line (called "stopping line" in [37]) to be a set £ = £{oS) c U x [0, oo), 
such that 

1. u 6 jV(t) for all (u,t) e £ and 

2. (u, t) e £ implies (v, s) $ £ for all v < u and s < t. 

Note that a line is at most a countable set. For a pair (u, t) 6 U x [0, go) and a line £, we write 
£ < (u, t) if there exists {v, s) e £, such that v < u and s ^ t. For a subset A c U x [0, oo), we write 
£ < A if £ < (u, t) for all (u, t) e A. If £\ and £2 are two lines, we define the line £\ a £2 to be the 
maximal line (with respect to <), which is smaller than both lines. 
We now define for each u e U two nitrations on Q u by 

&u{t) = (n u n a(E u (s);0 sj s sj t - b u )) v \/ (Q v n <r(H„)) 
#-P re (i) = (ft u n <t(H u (s); < a < t - &„)) v Y(fi« n <r(H v )). 

Informally, ^ u (t) contains the information on the path from u to the root between the times and 
t, and ^u VC {t) contains this information and of all the other particles excluding the descendants of 
u. In particular, we have ^ u {t) cz ^E re (t) The filtration ^ u {t) is denoted by srf u {t) in Chauvin's 
paper [22] and JFE rc (t) corresponds to the pre-(u,t)-sigma- algebra as defined by Jagers [37] . 

We can now define a stopping line ("optional line" in [37]) if to be a random line with the 
additional property 

3a. V(n,t) eU x [0,oo) : {we!] u : Jz? < («, t)} e ^ re (t). 
The sigma-algebra ^_^ of the past of J£ is defined to be the set of events £e^, such that for all 
(u,t) e U x [0,oo), 

£n{wE!l tl :^<(u,t)}eJn^ 

For example, for any t > 0, the set *sV(t) x {t} is a stopping line. This permits us to define the 
filtration (^tjt^o by 

^t = ^>(t)x{t}- 

Following Biggins and Kyprianou [9], we now say that if is a simple stopping line or a very simple 
stopping line, if it satisfies the property 3b or 3c below, respectively: 

3b. V(n, t) e U x [0, 00) : {oj e Q u : ££ < (u, t)} e Q u n #j. 

3c. V(n, f)e[/x[0,oo):jwEO u :i*< («, t)} e & u (t). 
Then ^V{t) x [t] is obviously a very simple stopping line. Furthermore, if T : D($) —* R+ is a 
stopping time, then 

&r = {(u,t) e U x [0,oo) : ue/(t) and £ = T(X U )}} 

is a very simple stopping line as well. We recall that the definition of stopping lines in [22] is 
equivalent to the definition of very simple stopping lines given here. 
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The first important property of stopping lines is the strong branching property. In order to state 
it, we define the subtree rooted at (u, t) for t ^ and u £ ^V(t) by 



LO 



(u,t) _ (,{u) c^i ,(«) 



(&\(E' uv ;vet^)), 



with H^(-) = E u (- + t — b u ) and r! uv = E uv for v e v- u '\{0}. The strong branching property ( [22"1 
Proposition 2.1], |37| Theorem 4.14]) then states that for every stopping line Jif, conditioned on 
^2>, the subtrees uj^ u,t \ for (u,t) £ _Sf, are independent with respective distributions P X «M. 

3.3 Many-to-few lemmas and spines 

Another important tool in the theory of branching processes is the so-called Many-to-one lemma and 
its recently published extension, the Many-to-few lemma [33] along with the spine decomposition 
technique which comes along with it and has its origins in [55], although it appeared implicitly in 
the literature before that, see e.g. the references in the same paper. Here we state stopping line 
versions of these lemmas, which to the knowledge of the author have not yet been stated in this 
generality in the literature, although they belong to the common folklore. We will therefore only 
sketch how they can be derived from the existing literature. 

We assume for simplicity that the strong Markov process X admits a representation as a con- 
servative strong Markov process X with paths in D(S'), which is killed at a rate R(x), where 
R : $ —* [0, oo ) is measurable. The law of X started at x is denoted by P x and the time of killing 
by £. Given a stopping time T for X, we can then define a stopping time T for X by setting T = T, 
if T < £ and T = oo otherwise. For simplicity, we write J£t for J^p. Finally, for every x e $, define 

m ( x ) = Hk^o( k ~ iM^ k )> rni(x) = Xifc^o k( l( x > k ) and m 2{x) = Yik^o k ( k ~ l )l{ x ^ k )- 

We are now going to present the spine decomposition technique, following [31] . They assume 
that q(x, 0) = 0, but this restriction is actually not necessary, as noted in [33]. Given a tree t, a 
spine of i is an element of the boundary of i, i.e. it is a line of descent £ = (£o = 0, £,1, £2, • • •) from 
the tree, which is finite if and only if the last element is a leaf of the tree. We augment our space 
£1 to the space Q* by 

0* = {(w, £) : to 6 CI, £ is a spine of the tree underlying oj} 

We are going to denote by £t the individual u £ U that satisfies u £ ^V(t) and u e £ if it exists, and 
it = otherwise. Instead of the redundant X^ t (t), we write X^{t). We also note that the definition 
of stopping lines can be extended to f2* by projection. 

Now, for every x £ $ , one can define a probability measure ~P*' X on £1* in the following way: 

• Initially, X^(0) = x. 

• The individuals on the spine move according to the strong Markov process X and die at the 
rate mi(y)R(y), when at the point y £ $. 

• When an individual on the spine dies at the point y £ S, it leaves k offspring at the point where 
it has died, with probability (mi(x))^ 1 kq(x, •) (this is also called the size-biased distribution 
of </(*,- fl. 

• Amongst those offspring, the next individual on the spine is chosen uniformly. This individual 
repeats the behaviour of its parent (started at the point y). 

• The other offspring initiate independent branching Markov processes according to the law ~P y , 
independently of the spine. 

This decomposition first appeared in [23] , 

If we start with n initial particles 1, . . . , n at positions xi, . . . , x n , we can extend this definition by 
defining a (non-probability) measure P*, which is the sum of n probability measures P* + • • • + P* , 



4 The size-biased distribution of the Dirac-mass at is again the Dirac-mass at 0. 
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where under P*, the particle i follows the law P*> Xi and the remaining particles j ¥= i follow the 
law P X J . 

We now have the important 

Lemma 3.1 (Many-to-one). Let Jz? be a simple stopping line. Define T by (£r,T) 6 j£f if it exists, 
and T = oo otherwise. Let Y be a random variable of the form Y = 2(ut)eif ^u^-(ueQj where Y u an 
^■j^ -measurable random variable for every ueU . Then 

~ Ye llR(X i {t))m(X i (t))dt : 



E a 



L (T<oo) 



(3.2) 



E[ 2 *," 

(«,£)ejz? 

Proofs of this result can be found for fixed time in [32], [3]] or [33]. With simple stopping lines, 
it has been proven in the discrete setting [9\ Lemma 14.1] and their arguments can be used to adapt 
the proofs in [3]] and [33] to yield the result stated here. 

Often, we will use a simpler version of the Many-to-one lemma, which is the following 

Lemma 3.2 (Simple Many-to-one). Let T = T{X) be a stopping time for the strong Markov process 
X which satisfies P X (T < oo) = 1 for every x e <f. Let f : $ — * [0, oo) be measurable. Then we have 

E x \e& R ( x t) m ( x t) dt f(X T ) 



E*[ 2 f{X u (t)) 
(u,t)e3f T 



The next lemma tells us about second moments of sums of the previous type. To state it, we 
define for a stopping time T for X the density of the branching Markov process before J£t-, by 



p T {x,dy,t) = E 



[S 



l (X u (t)edy, t<T(X u )) 



(3.3) 



" ue^V(t) 

Lemma 3.3. Let H be the hitting time functional of a closed set F a $ on D{$) which satisfies 
P X (H < oo) = 1 for every x e <§ . Let f : $ — * [0, oo) be measurable. Then we have 

2n 



E J 



2 f(Xu(t)) 
(u,t)eJ? H 



E*[ £ (f(X u (t))f 

"00 



+ r f p H (x,dy,s)R(y)m 2 (y)M £ /(X u (i))j) 2 ds (3.4) 

Remark 3.4. The intuitive explanation of Lemma |3.3| (see, e.g. the proof of Proposition 18 in [7J) 
is that a particle at the point y spawns an expected number of R(y)m2{y) ds of ordered pairs of 
particles during the time interval [s, s + ds], and by the strong branching property, the two particles 
in a pair evolve independently. This lemma can be proven using the Many-to-few lemma from [33] 



(which is valid for stopping lines as well by the same argument as the one above) or with Lemma 3.1 
by noting that 

2 f{X u (t))f= 2 {f{X u {t))) 2 + 2 (/(*,(*)) S K x ^ 

(u,t)e3f H ' (u,t)e^ H (u,t)e^ H (v,s)e£f H , v^u 

Taking for X the space-time process (Yt,t)t^o of a possibly non-homogeneous strong Markov 
process (Yt)t^>o with paths in D(S') and the closed set F = $ x {t}, for some t ^ 0, we obtain the 
following useful corollary, which appeared already in [57] and |54] in the homogeneous case. 

Lemma 3.5. Let f : t§ x R + — > [0, oo) be measurable and let t ^ 0. Then we have 

e ( *' 0) [( 2 /(y u (t),t)) 2 j=E^°)[ 2 {f{Y u {t),t)f 



ue^V(t) 



' ue^V(t) 



f f P (x,dy,s)R(y,s)m 2 (y,s)(E ( -y^\ £ f(Y u (t),t)]Y ds (3.5) 



MG^(t) 
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4 BBM with absorption at a critical line 

From this section on, q{k) will denote a law on {0, 1, 2, . . .} and L a random variable with law q{k). 
We define m = E[L — 1] and rri2 = E[L(L — 1)] and suppose that m > and 771,2 < 00. We 
study the branching Markov process where, starting with a single particle at the origin, particles 
move according to standard Brownian motion with drift — 1 and branch at rate j3q = l/(2m) into k 
particles according to the reproduction law q{k). At the point —y, we add an absorbing barrier to 
the process, i.e. particles hitting this barrier are instantly killed. Formally, we are considering the 
process up to the stopping line J2f#_ y , where H-y is the hitting time functional of the point — y. It 
is well-known since Kesten [39] that this process gets extinct almost surely. As a consequence, the 
number of particles absorbed at the barrier, i.e. the random variable 

is almost surely finite. By the strong branching property and the translational invariance of Brow- 
nian motion, one sees that the process {N y ) y ^,Q is a continuous-time Galton-Watson process, a fact 
which was first noticed by Neveu [51] (see [I], Chapter III or |34j . Chapter V for an introduction 
to continuous-time Galton-Watson processes). Let u{s) be its infinitesimal generating function. 
Neveu stated that u = ifi' o tp~ , where ip is a so-called travelling wave of the FKPP (Fisher- 
Kolmogorov-Petrovskii-Piskounov) equation: Write f(s) = ^ i i c s k q(k). Then tp is a solution of the 
equation 

i<-^ = A>ty-/oV0, (4.1) 

with ip(— 00) = 1 and ifj(+cc) is the extinction probability of the process, i.e. the smaller root of 
f(s) = s. For a proof of these results, see [Ml Section 3]. 

In the same paper [51] , Neveu introduced his multiplicative martingales, which he used to derive 
the Seneta-Heyde norming for the martingale e~ y N y . He proved that in the case of binary branching, 
one has 

W y := ye~ y N y — ► W almost surely as y -* 00, (4.2) 

where W > almost surely. His proof relied on a known asymptotic for the travelling wave ip, 
namely that 

1 — ip{— x) ~ Kxe~ x , as x —*■ 00, (4.3) 

for some constant K > 0. It was recently shown |59| that this asymptotic is true if and only if 
-E"[Llog L] < co and the proof of (4.2) works in this case as well. We also still have in this case, 



for every ieR, 

E[e- eXW ]=i;(x), (4.4) 

a fact which was already proven by Neveu |51j for dyadic branching. 

In [7], further properties of the limit W have been established under the hypothesis of dyadic 
branching, namely 

P(W > x) , as x -* 00, (4.5) 

x 

and 

E[W r l (w ^ a .)] - log x -» qo], as x -* 00, (4.6) 



for some constant qj^] £ R. Equation (4.5) has been proven in Propositions 27 and 40 of [7], and 
(4.6) appears in the proof of Proposition 39 of the same paper. Their arguments were very ingenious 
but indirect and although they could be extended to general reproduction laws with finite variance, 
we will reprove them here directly under (probably) minimal assumptions, based on methods of 
51. The main result in this section is 



Proposition 4.1. If E[L log L] < 00, then (4.5) holds. If E[L log 6 L] < 00, then (4.6) holds. 
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See also |20| for a proof of (4.5) in the case of branching random walk. Before proving this 
result in the next subsection, we state a lemma which is immediate from (4.2) and the fact that N y 
is almost surely finite (see also Corollary 25 in |7J): 

Lemma 4.2. Suppose -E[Llog L] < oo. For any rj > 0, there exist y and Q, such that y > r/ _1 and 
P(\W y -W\ > v)+P(N y > C)+P(^H_ a $ C/x[l,C])+P( sup #{«/(*) : (u,t) < J? H _ y } >Q<V- 



4.1 Proof of Proposition |4.1 

Define xO) = E[e~ xw ] for A > 0. Our first result is: 

Lemma 4.3. Suppose E[L log 2 L] < oo. T7ienx"(A) ~ A -1 as A — » 0+. Furthermore, E[Llog 3 L] < 
oo if and only ifr(X) = A -1 - x"W > for A > 0, wwt/t Jo r (^) ^ < °°- 



Proof. Define <p(x) = 1 — ip(— x), such that u(s) 
E[Llog 2 L] <oo, 

6(x) ~ Kxe~ x , 



as a; — »• oo 



1 (s)). By (4.3) and the hypothesis 

(4.7) 



Furthermore, by (4.1), we have 



1 



<t>"(x) + tf/(x) = A)(/(l " 0(x)) - (1 - <Kx))). 



(4. 



Setting 5 (s) = 2/3ol/(l - a) - 1 + /'(l)a] ^ and p 
Po = l/(2m) and /'(l) = m + 1 by definition, 



6', we get from (4.8), and the fact that 



p'{x) = -p{x)+g{<j>{x)). (4.9) 

As in the proof of Theorem 1.1 of |46| . we will study the function p through the integral equation 
corresponding to (4.9), namely 

p(x) = e-( /3 (O) + J X e ^(0(y))dy 



P(0) + 

Jd>(x 



4>(o) p^-Hs) 



9(*) 



ds . 



)0 ' ' J<p(x) "«(*) 

Now, by Theorem B of [11] (see also Theorem 8.1.8 in [12]) we have for every d > 



(4.10) 



f 

Jo 



— ^ 5 ( S )ds<00 ^ — - ^ 



</(s) ds < oo 



£[Llog i+d L] <oo. (4.11) 



Furthermore, by Proposition 3.2 of [46], we have — u(s) ~ s as s — > 0, and by (4.7), we have 



o4>~\s) 



(log l/s)/s as s — > 0. By the hypothesis i?[L log -L] < oo, this gives 



Jo 



^ 1(s) g(^) 

-«(s) 



ds < oo, 



whence, by (4.10), 



p(x) 



Ke 



as x — > oo, 



(4.12) 



where the constant if is actually the same as the one in (4.7), see again the proof of Theorem 1.1 



log A), whence, by ( |49| ) and ( |4.10[ ), 

00 



of [46]. Now, from ( [44 ), we get x(A) = 1 — 0( 

x"(A) = -^p / (-logA) = x + ^(-Aj e^(0(i/))di/-^(-Io g A))J 

f eV(y)5 / (0(y))dy), 



if 1 

A + A 



(4.13) 



-log A 
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where the last equation follows from integration by parts. This proves the first statement, with the 
constant K instead of 1, since the last integral vanishes as A — > 0. Now, setting 

rW = ~\(T e y cf>'(y)g'(<t>(y))dy] 



log A 



we first remark that r(A) > 0, since the integrand is negative for y e R. By the Fubini-Tonelli 
theorem, we then have 

rl rl i rCO rCO 

r(A)dA = -- eV(y)</(<My))dydA= -ye y 4>'(y)g'(<j>(y))dy 

JO JO A J- log A JO 

r<t>(°) _-, 

e^^c/ ) - 1 (y)g'(y)dy, 
Jo 



which is finite if and only if E[Llog 3 L] < oo, by ( |4.11[ ) and the fact that e^ ^<f>~ 1 (y) ~ (log 2 l/s)/s. 
This proves the second statement, again with the constant K instead of 1. 

The previous arguments worked for every travelling wave ip. In order to show that that the 
constant K is equal to 1 in our case, we use Neveu's multiplicative martingale (this idea was also 
used in |44| . Theorem 2.5). It was observed by Neveu [51] (see also p2j for a rigorous proof), that 



((1 — 4>(x + y)) Ny ) y ^o is a martingale for every ieR with values in [0, 1]. By (4.2) and (4.7), we 
then get by dominated convergence, for every x € R, 



X (Ke x ) = lim E[e- KyeX yJV «] = lim E[(l - <f>(y - x)) N «] = 1 

y— >oo y—>cc 



X{e 



This yields K = 1. 



□ 



Remark 4.4. Choosing arbitrary initial points xq,x\ e R instead of in (4.10), one sees that 

r.Q0 



J '00 pOO 

e y g{(t){y)) dy = e xi p(xi) + e y g((f)(y)) 
XQ Jx± 



dy. 



In particular, since p is bounded, letting xo — > — oo and x\ — * +oo yields J_ e y g(4>(y)) dy = 1. One 
could hope (see the proof of Proposition 4.1 below) that this helps in determining the constant qjj], 
but apparently this does not seem to be the case. 



Proof of Proposition J^.l, For n e N, we define the function 

V n (x) = f y n P(W e dy) = E[W n l {w ^ x) ], 
Jo 

such that with x denoting the n-th derivative of x, we have for A > 0, 



x (n) (A ) = (_!)« 

Jo 



-Ax 



; dVJx). 



If i?[L log L] < oo, Proposition 4.1 and Karamata's Tauberian theorem ([29J, Theorem XIII. 5. 2 or 
|12j . Theorem 1.7.1) now yields 

V2(x) ~ x, as x —* oo. (4.14) 

By an integration by parts argument (see also [29], Theorem VIII. 9. 2 or [12) . Theorem 8.1.2), we 
get (4.5). Now suppose that i?[L log L] < oo. By Lemma 4.3, we have x'(A) — logA — * c e R, as 
A — > 0. By Theorem 3.9.1 from [12] (with £(x) = 1), this yields 



V\(x) — logx — ♦ 7 



as x — »• oo, 



where 7 is the Euler-Mascheroni constant. This is exactly (4.6). 



□ 
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5 BBM in an interval 

In this section we study branching Brownian motion killed upon exiting an interval. Many ideas in 



this section (except for Section 5.5 and parts of Section 5.3) stem from Sections 2 and 3 of [7] and 
for completeness, we will reprove some of their results with streamlined proofs. However, we will 
also extend their results to the case of Brownian motion with variable drift. 

5.1 Notation 

During the rest of the paper, the symbol C stands for a positive constant, which may only depend 
on the reproduction law q. Its value may change from line to line. If a subscript is present, then 
this subscript is the number of the equation where this constant appears for the first time (example: 
C J5~36l) . In this case, this constant is fixed after its value has been chosen in the corresponding 
equation. If X is any mathematical expression, then the symbol 0{X) stands for a possibly random 
term whose absolute value is bounded by C|X|. 

Recall the definition of q{k), m, rri2 and /3o from Section Hand the hypotheses on m and 7712. In 
this section, we let a ^ tt and set 



M = \/l--o- (5-1) 



From (5.1), one easily gets the basic estimate 



a 2 



7T 2 



0<l-/i<— ^. (5.2) 

2/icr 

We then denote by P 1 the law of the branching Markov process where, starting with a single 
particle at the point x e R, particles move according to Brownian motion with variance 1 and drift 
— /j, and branch at rate 1 into k particles according to the reproduction law q{k). Expectation with 
respect to P x is denoted by E 31 . On the space of continuous functions from R + to R, we define Hq 
and H a to be the hitting time functionals of and a. We further set H = Hq a H a . Then note that 



the density of the branching Brownian motion before Jz?h, as defined in (3.3), has a density with 
respect to Lebesgue measure given for t > and x, y e (0, a) by 

p t (x,y) = e^ x - y)+ ^ t p?(x,y), (5.3) 



where pi was defined in (2.3) 



Now, let / : R + — ► R + be non-decreasing, with /(0) = 0, continuous and such that the left- 
derivative /' exists everywhere and is of bounded variation. Such a function will be called a barrier 
function. We define 

11/11 = maxjll/Hoo, Unioo, liril^, J°° Id/'Cs)!}, (5.4) 

where || - 1| 00 is the usual supremum norm. Furthermore, we set Err(/, t) = ||/|| ( ^ + \ J . Now define 

Vt = V + 4/(i/a 2 ) = \i + ^f'(t/a 2 ), (5.5) 

at a z 

such that fj,Q = fi and /j,t ^ \i for all t > 0. We denote by P x r the law of the branching Brownian 
motion described above, but with infinitesimal drift — /xj. Expectation with respect to P^ is denoted 

by E^ and the density of the process is denoted by p{(x,y). 

The above definitions can be extended to arbitrary initial configurations of particles distributed 
according to a counting measure v on (0, a). In this case the superscript x is replaced by v or simply 
omitted if v is known from the context. 
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5.2 The processes Z t and Y t 

Recall from Section K^ that the set of particles alive at time t is denoted by .jV(t). We define 

^6,«(t) = {u S Jf{t) : H{X U ) > t}, 

where H was defined in the previous subsection. Now set wz(x) = ae^ x ^ a > sm(irx/a)l( xe [ 0a ^ an d 
wy(x) = e^ x ~ a ^lf x ^Q-\ and define 

Z t = 2 w z {X u {t)) and Y t = £ u* (*«(*)) • 

ue,jV ,a{t) «6^8, (t) 

Then ijj is a martingale under P x , the proof of which is standard and relies on the branching 



property, the Many-to-one lemma (Lemma 3.2) and the fact that e t ' 2 wz{Bt) is a martingale for a 
Brownian motion with drift — jjl killed at and a, which is easily seen by Ito's formula, for example. 
Furthermore, it is easy to see as well that Zt is a supermartingale under P?. 

The following lemma relates the density of BBM with variable drift to BBM with fixed drift. 

Lemma 5.1. For all x,y e [0, a] and t #= 0, we have p((x, y) = pt(x,y)e^^ t ' a )+°( En W)). 

Proof. By the Many-to-one lemma and Girsanov's theorem, we have 

p{(x,y) = e ^ mt W^ t (B t e dy, H > t) 

= exp (t/2 - f ^^ds\ W% f exp ( - f fi s - fx d£? s ) ,B t edy,H>tY 

By the integration by parts formula J jx s dB s = fifBt — [iBq — J B s d/x s and some easy estimates, 
we get with (5.3), 

p{(x,y) = p t (x,y)e-^ t ^ +0 ^ E ^M\ 



The lemma now follows from ( |5.2| . D 

Proposition 5.2. Under any initial configuration of particles, for every t ^ 0, we have 

B f [Z t ] = Z e-fW^ + °( E "^\ (5.7) 

and if in addition \i ^ 1/2, then 

Var / (Z t )^Ce-^/ a2 ) +0 ( Err ^)(^Zo + y ). (5.8) 

Furthermore, we have for every t > (without hypothesis on /i), 

B f [Y t ] sj Ce-M/^+^^Yo. (5.9) 

and for t > a 2 , 

E f [Y t ] ^ Ce-^ i / a2 ) + °( Err ^)) Zo . (5.10) 

Moreover, for every a 2 < t ^ a 3 , we have 

Vax f (Y t ) < Ce-^l^ ^^. (5.11) 
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Proof. Equation (5.7) follows from Lemma 5.1 and the fact that Zt is a martingale under P x . In 



order to show (5.9) and (5.10), it suffices by Lemma 5.1 to consider the case without variable drift. 
We first suppose that t ^ a 2 . By (5.3) and (2.7), we get 



E x [Y t ] sc e^ x -°^ [ e^ t Pt {x,y)dy^Ce^ x - a hm{TTx/a) f - sin(vry/a) dy. 
Jo Jo a 



The last integral is independent of a. Summing over x yields (5.10) as well as (5.9) in the case 
t > a 2 . Now, if t < a 2 , by the Many-to-one lemma and Girsanov's theorem, we have 

W[Y t ] = e 0omt w x\ e n(x t -a)^ H AH a <t\= e n2t/{2a2) W x [H A H a < t]e^ x ~^ . 
Summing over x yields (|5.9|) 



In order to prove (5.8), we have by Lemma 3.5 



E x f [Z 2 ]=B x f [ J] w z (X u (t)f 

ue.yY ^{t) 



/3 m 2 



f fp f s (x,y)(Bf> s) [Z t ]) 2 dsdy. (5.12) 



By Lemma 5.1 and the fact that Zt is a martingale with respect to the law P x , this yields 

E *[Z t 2 H Ce-^/ ffl2 ) + °( E "(M) (e*[ V w z (X u (t)) 2 ]+ n t p s (x,y)w z (y) 2 dsdy). (5.13) 

V u e ^ a{ t) J Jo Jo J 

Now we have wz{x) 2 = (asin(Trx/a)e~^ a ~ x )) 2 < 7r 2 (a — x ) 2 e~ 2 ^ a ~ x ^ ^ Cwy{x) for x 6 (0,a), 
because \x 5= 1/2 by hypothesis. This yields 



Si := E*[ 2 w z {X u {t)f\ sj CE x [Y t ] sc Cw Y (x), 
ue<Ao,a(t) 



(5.14) 



by (5.9). Now, by (5.3) and (2.11), we have 



S 2 := f f Ps (x,y)w z (y) 2 dsdy = ae^ x -^ f ae^^ sin 2 fry /a) J a (x, y, t) dy. 
Jo Jo Jo 



Lemma 2.1 



now gives 



S 2 sJ Cae^ x - a) f i 
Jo 



e w sin (iry/a) ( £ sin(7rx/a) sin(7ry/a) + ay ) dy 



< 



t 1 



Cae^- a )(sin(7nc/a)4r + 



e-^y 3 dy, 



(5.15) 



the last line following again from the change of variables y *—> a — y and the inequality sinx ^ x. 



Using again the fact that [i ^ 1/2, equations (5.13), (5.14) and (5.15) now imply 



E^[Z 2 ] ^ Ce-fW^+^^^wzix) + wy(x) 



(5.16) 



If we write the positions of the initial particles as x\, . . . ,x n , then by the independence of their 
contributions to Z t , 

\ a v f (Z t ) = £ VaijW) < 2 E / 1 ^ 2 ]- ( 5 - 17 ) 



Equations (5.16) and (5.17) now prove (5.8). Equation (5.11) is proven similarly. 



□ 
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5.3 The number of particles 

In this subsection, we establish precise first and second moment estimates for the number of particles 
alive at a time t. These estimates extend those of [7], which are effective only when t » a 2 . For 
r e [0, a] and t ^ 0, we denote by Nt(r) the number of particles in [r, a] at time t. 

Proposition 5.3. Suppose u ^ 1/2. Let t > ; x, r e [0, a) and suppose that x ^ (r + a/20)l( t ^ o a) . 
T/ien 



E*[7V t (r)] = e 



/J,(x— r) 



2(1 + fir)e n */( 2a i^/i t 
vd' a 



„ sin(7rx/a) v lr+<„2i 
, ; , O ( (1 + r 3 ) V 7 ^ (4 ^ a ) 



(5.18) 



In particular, if a/20 + rl(^ a 2) ^ x < a, £/ten 



E x [JV t (r)] 



2tt(1 + /jr)e^ a - r ) 



wz(x) 



' (a-x) 2 + (l + r 2 ) / wg(g) 

I i V Wy (2)1(^2) 



Moreover, for every x e (0, a) and t > 0, we have 

E*[JV«(r)] = [ 3 

a- 3 

wi/i | error | sj E t / a 2(l + 0(l/a 2 )). 



wz(x) ( 1 + error) , 



(5.19) 



(5.20) 



Proof Let £ ^ and ^ x < 1. By (5.3) and Brownian scaling, we have 

B ax \N t 2 (r)l = e ^ ax +' w2t l 2 I ,-r"'rJi 



J 

Jr/a 



z Pt (»i *) dz. 



(5.21) 



Note that (d/dz)p\(x 1 z)| z= o = —29'(x,t), by (2.4). Taylor's formula then implies that there exists 

£ e [0, z], such that 

<9 2 
pKx,*) - -2z6'(x,t) + z 2 1 ^p 1 t (x,z)\ (5.22) 



dz 2 



We first note that for all x e (0, 1) and t ^ 1, or for all x £ [^g, 1] and t < 1, one has by (2.1 ) and 

|6>'(x,t)| ^ C(sin(7rx) v l( t< a))- 



(2.2), 



(5.23) 



Furthermore, if t ^ 1, we have 

CI 2 



= it 2 V e 7r2n2 ' //2 n 2 sin(7rnx) sin(7rn£) 

n=l 

oo 

< Asin(Trx) 2 ^e"" 2 ™ 2 * 72 < Czsin(7rx)e-" 24/2 , 



(5.24) 



n=l 



by the inequality | sinnx| =% nsinx, x £ [0,7r]. Equations (5.22), (5.23) and (5.24) now give 



J ' 

Jr/a 



,-(iaz 1 



p t (x, z) dz 



2 2 

[i A a L 



l + He^lO (l + r 3 )^ 



. Sin 7TX 



3 -ir J */2 



(5.25) 



Now suppose that t ^ 1, x, r 6 [0, a) and x ^ r + a/20. By (2.4) and the mean value theorem, there 

exists £' e [0, z], such that (<5 2 /dz 2 )pJ(x, z) = 2z9'"(x + £',£). With ( |2.2[ ), one then easily sees that 

for z < x — a/40, one has 

(9 2 

^(x,z) < Cze^ c felll /4 . (5.26) 



dz 21 
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Equations ( |5.22[ ), ( |5.23[ ) and ( |5.26[ ) now give, for x' = x — a/40, 

29'{x,t) 



rx' /a 
Jr/a 



D —[iazl 



r/a 
Furthermore, we have 



p t (x, z) dz 



fi 2 a 2 



[1 + ^)6-^ + 



l + r 3) e -A»r 



I e-» az p\ {x, z) dz sc e-^' ( p] [x, z) dz 

Jx'/a Jx'/a 



< e -M* 



(5.27) 



(5.28) 



Equations (5.21), (5.25), (5.27), (5.28) and the hypothesis on /i now imply (5.18). 



In order to prove (5.19), let x e [1/20,1]. By Taylor's formula and the fact that 9(x,t) is a 



solution to the heat equation with diffusion coefficient 1/2, there exists £ e [x, 1], such that 

- e'(x,t) = 2(1 - x )~e(i, t) - (i - x) 2 9'"(a,t). 

With (|2.1[) and (|2.2[), one now easily sees that \8'"(£,t)\ s$ C(e _7r <//2 sin(7rx) v l(t<i)). Equations 



(5.29) 



(|5.2|), (|5.18|) and (|5.29|) now readily imply (|5.19|). For the last equation, by (|5.3|) and ( |2.7[ ), we have 

(5.30) 



for every x 6 (0, a) and £ > 0, 

2 e A"J r a 
B x [N(ta 2 )] = (1 + err)u; z (x)^- e^sin(7ry/a) dy, 

a Jo 



with |err| ^ Et- Evaluating this integral and using (5.2) and the hypothesis on [i yields (5.20). □ 



We will often use the following handier upper bounds on E^fiV^r)]: 
Lemma 5.4. Suppose fi ^ 1/2. Let t ^ 0, x,r e [0,a] and suppose that x ^ (r + a/20)l( t ^ a 2\. 



e n{a-r) ( 1 + r 2 



E x [N t (r)] ^ C(l + r) e —^(w z (x) + -^L w (x)l ( ^ a2) ) . 



a" v a 

Furthermore, we have for all x, r e [0, a] and t ^ 0, E^fiV^r)] ^ Ce^ x ~ r > . 



(5.31) 



Proo/. One sees from (2.1) and (2.2) that 9"{x,t) «S C for (x,t) e [0,1] x [l,oo) u [1/20,1] x [0,1]. 



Equation (5.31 ) now follows from (5.18) and the mean value theorem. The second statement on the 

□ 



other hand follows from the inequality Nt(r) ^ e^ a r >Yt together with ( |5.9[ ). 

Lemma 5.5. Suppose [i ^ 1/2 and r =% 9a/10. For every t ^ and x e [0, a], we have for large a, 



B X f [N t (r) 2 ] Sj C e -/(*/« 2 )+0(Err(/,t)) 



Z1 fia 



a" 



„ / 1 + r 4 
I + r 2 1 + =— ir 



2 ^(u; z (x)^+w(x)) (5.32) 



Proof. As in the proof of (5.8), we have by Lemmas 3.5 and 5.1 

"a rt 



E x [N t (r) 2 ] sc e- / ( t /° 2 ' +0 ( Err (W)(E s: [JV t (r)]+ J 8om2 [ f p s (x,z)(W[N t _ s (r)]) 2 dsdz). (5.33) 



By (5.3), Lemma 5.4 and the hypotheses on x and r, we have after a change of variables z — > a — z 
in the integral, 



f 

Jo 



p s (x,z)(E z [N t _ s (r)]) 2 dz^C 



x i(l + r 



Jo 



,/«( 



fi(x-a-2r)+TT 2 s/(2a 2 ) 



a- x,z)e ^ z ( z 2 + 



1+r 4 



dz + a 6 ( p a s (a-x,z)e- flz dz). (5.34) 

Ja/20 ' 
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Integrating ( J5.34 ) over s from to t and splitting the interval at t a a 2 , we have by ( |2.8[ ) (for the 
first piece) and (2.7) (for the second), 



ft 

Jo Jo 



Ps (x,z)(B z [N t _ s (r)]) 2 dsdz^C 



,/"< 



„ / 1 + r 4 
x 1 + r 2 ) 1 + s— |( 



>(l + oV" fl / 20 )( Wz (4+wW) (5.35) 



The lemma follows from (5.33) and (5.35), together with the second statement of Lemma (5.4) and 
the hypotheses on fi and r. □ 

5.4 The particles hitting the right border 

In this section we recall some formulae from [7J about the number of particles hitting the right 



border of the interval. We reprove these formulae here for completeness and because Lemma 2.1 
makes their proofs straightforward. For most formulae we will assume that / = 0, i.e. that we are 



working under the measure P. Only Lemma 5.7 contains an upper bound on the expected number 
of particles for general /, which will be useful in Section [7j 

For a measurable subset S c R, define Rs to be the number of particles killed at the right 
border during the (time) interval S, i.e. 

R s = #{(«,*) : u 6 Jf(t) and H (X U ) > H a {X u ) =teS}. 

Lemma 5.6. For any initial configuration v and any ^ s < t, we have 



\B[R M ] - ^ 3 S) Z Q \ sj %3gl(Yo a E s/a2 (l a (t - s)/a 3 )Z ), 



(5.36) 



where E s is defined in (2.6). Furthermore, if /x > 1/2 and ^ t ^ a 6 , then for each x e (0,a), 



B x [Ri]^G^ m [^w z (x) + w Y (x)), 



(5.37) 



Proof. Recall that for a real- valued process, Hq and H a denote the hitting time functionals of and 
a and H = Hq a H a . Note that W*AH < go) = 1 for all x e [0, a]. By Lemma 3.2, we then have 



E X [Rs] = ^ X [ J] 1 (X u (t)=a,teS) 
(u,t)eJ? H 



w x 



B mH a -i , 
& *-(H >H a eS) 



^ x - a h a (x,S), (5.38) 



where the last equality follows from Girsanov's transform and (2.10). Equation (p. 36) now follows 



Lemma 2.1 The second-moment estimate is obtained from Lemma 3.3 and ( 5.38[ ) by a calculation 
similar to the one in the proof of Proposition 5.2 (see [4"Tl Lemma 5.8] for details). □ 

Lemma 5.7. Let f be a function as in Section \5. 1\ Then for every x e (0, a), we have 

B}[R S ] < B X [R S ]. 



Proof. As in the proof of Lemma 5.6, we have 



E x f [R s ] = Wl 



i-'t 



e ±(H >H a eS) 



w a n 



-tf a HtdB t -tf a ^/2dt+H a /2 1 



(H >H a eS) 



by Girsanov's theorem and the definition of f$Q. On the event {Hq > H a }, we now have by integration 

by parts, 

rH a rH a 

fit dB t = /x(o -x) + a{fi Ha - fi) - B t dfi t > K a ~ x )> 

Jo Jo 
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since B t 6 [0, a] for t e [0, H a ]. This gives 



by the proof of Lemma 5.6 □ 



We finish this section with a lemma which links BBM with absorption at a critical line to our 
BBM with selection model. The proof is simple, see \47\ Lemma 2.5.10] for details. 



Lemma 5.8. Let £ > 1, y > 1, /j, ^ 1/2 and f be a barrier function (defined in Section 5.1). 
Suppose that ^ t /a ^ y + ( and \\f\\ ^ -y/a. Let (xi,ti)f =l be a collection of space-time points with 

Xi = a- y + (1 - p)U - f(s/a 2 ), i=l,...,N, 

and U ^ C / or o" *■ Define Z = ^ u>z(xj), 1" = £^ w(^«) °^ W 7 ?/ = ye~ y N. Then, 

2-^,(1 + 0(1)) „„, y _ l^(i + o(l)). 

In particular, for large a, we have 

Y ^ Z/y. 

5.5 Penalizing the particles hitting the right border 

In this section, let (U u ) ue u be iid random variables, uniformly distributed on (0,1), independent 
of the branching Brownian motion. Furthermore, let p : R + e (0, 1] be measurable and such that 
p(t) = for large enough t. Recall that H = Hq a H a . We define the event 

E = {$(u,t) e ££ H : X u {t) = a and U u ^ p(t)}. 

In order to describe the law P^ = P x r(-\E), set 

h{x,t) = pf' t \E), Q{x,t) = J] qifyhfat)*- 1 , q(x,t,k) = q(k)h(x,t) k - l /Q(x,t) (5.39) 

The law P^ can be seen as a Doob transform of P^ by the space-time harmonic function h(x,t). 

Under P?, the BBM stopped at J£h is again a branching Markov process where (see |47| Section 3.4] 
and the errata in the arXiv version) 

• particles move according to the particles move according to the law obtained from Brownian 
motion with drift —fit (stopped at and a) through a change of measure by the martingale 

tAH 



h(X tAH ,t a H) exp (- f " /3 (1 - Q(X„ s)) dsj . 



• a particle located at the point x e (0, a) at time t branches at rate /3qQ(x, t)l xe fo,a)j throwing 

k offspring with probability q(x, t,k). 
Together with Lemma |3.2| this immediately gives the following useful Many-to-one lemma for 
the conditioned process stopped at the stopping line Jzfj = «5fjj A t : Define the function 

e(x,t) =PoJ]( k + !)(! ~ h&tf-^qik) < /3 (m 2 + m)(l - h(x,t)). (5.40) 
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Lemma 5.9. For any measurable function g : [0, a] — ► R + , x e (0, a) and t > 0, we have 



& f [ 2 5 (X u (i))' 
veSft 



w° 



l-'l 



h(x,0) 



ds 



(5.41) 



In particular, if we denote by p((x,y) the density of the P x r-BBM, then 



Pt(x,y) = 1 $4pt^y)w t x ' t ' y 



h(x,oy 



taboo 



,-ll e ( x s,s)ds 



<(/i(x,0)) l pt(x,y) 



and for general f , 



Pt(x,y)^(h(x,0)) l p{(x,y). 
This lemma immediately gives an upper bound for the quantities we are interested in: 



(5.42) 
(5.43) 



Corollary 5.10. Let x e (0, a), t ^ and g : [0, a] — ► R + be measurable with g(0) = g{a) = 0. 
Define S t = Xm Gj% g{X u {t)). Then, 






(5.44) 
(5.45) 



Proof. Equation (5.44) immediately follow from (5.42). In order to prove the second-moment esti- 



mates, we note that by Lemma |3.5| and the description of the conditioned process, 



W f \ 



[S?] = Ej\ V g(X u (t)f] + f ( a p s (x,y)m 2 (y,s)(3 Q(y,s) (eJ^S*]) 2 dyds, 



where fh~2(x,t) = 2fe>ok(^ ~ ^)<l( x 't,k). By (5.39), we have rri2{x,t)Q(x,t) ^ h(x,t)rri2- Equation 
flBT45] ) then follows from ([5I2J), ([542]) and ( |5.44| >. D 



The following lemma gives a good lower bound on the first-moment estimates in the case where 
/ = 0. 



Lemma 5.11. Suppose fi ^ 1/2, t < a 3 and p(s) = /or a// s ^ a 3 , Lei £* fre as in Corollary 5.10. 
We have 

B x [S t ] > B x [S t ](l - CfczJplD- (5.46) 



This follows from the following estimate on h(x, 0), which will be sharpened in Lemma |6.3 
Lemma 5.12. Suppose p(s) = for all s ^ a 3 . Then for all x 6 (0, a), we /iawe 

1 - h{a - x, 0) sc C||p||ao(x + l)e~^ x . 
Proof. By Markov's inequality, we have 

1 - h(x, 0) sj E^(#{(n, s) e ^ a z : X tt (a) = a, [/„ < \\pU}) sj \\pUW f {R a z). 



The lemma now follows from Lemmas 5.6 and 5.7 and the inequality svnx ^ x, x 6 [0, 7r]. 



□ 



Proof of Lemma 5.11\ By ( 5.42 ) , 

W[S t ]>W[S t -\ inf (MV,*)WS£ 



J/ 

taboo 



-J* e(X s ,s)ds 



(5.47) 
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By Lemma 5.12 and the hypothesis on //, we have for every y £ (0, a), 

h{a -y,t)> h{a -y,0)>l- C|H|ooe^ /3 . 
This gives 



^ taboo [ e 



-J*e(X s ,s)ds 



> i - wsa, 



f e(X S)S ) 

Jo 



d.s 



^ 1-C||p||oo, 



(5.48) 



(5.49) 



by (5.40), (5.48), Lemma 2. 2 1 the inequality e x ^ 1 — x for x > and the fact that the law of 



the Brownian taboo process is preserved under the map y h-* a — y. The lemma now follows from 
dOTI), KM and KM. D 



Finally, we study the law of Rt under the new probability. 
Lemma 5.13. We have for every x e [0,a], 

E x f [R t ] - b||ooE^[« t 2 ] ^ E x f [R t ] ^ (hix^r^Rt], 



(5.50) 



and if there is a p e [0, 1], such that p(s) = p for s < t, then we even have 

B][R t ]^E f [R t ]. (5.51) 

Proof. Define the stopping line §£% = {(u, s) 6 ^£n a ■ s ^ t}. By the definition of the law P, 



m[Rt\ 



E-' 



E / 



n ( «, a )ea t (l-P(-X-«W)) 



(5.52) 



Now the denominator is /i(x, 0) by (5.39), which yields the right-hand side of (5.50). The left-hand 
side follows by noticing that 



E 



Rt Y\ (l-p(X u (s)))] > EJ[R t (l - \\ P U Rt ] > W f [R t ] - b||ooE^LR t 2 ]. 

(u,s)sMt 



For (5.51), we note that if p(s) = p for s < t, then by ( |5.52 ), 

W f [R t (l - p) R t] 



E%R t ] 



E?[(l-p)**] ' 



Since (1 — p) is decreasing in k, this yields (5.51). 



D 



6 BBM with absorption before a breakout 

In this section, we study branching Brownian motion with drift —fi and absorption at until a 



breakout occurs, an event which will be defined in Section 6.1 and which corresponds to a particle 



going far to the right and spawning a large number of descendants. In (6.20), we decompose the 



system into a particle conditioned to break out at a specific time T (this particle will be called 
the fugitive) and the remaining particles, which are conditioned not to break out before time T. 



These two parts will be studied separately, the former in Section 6.4 and the latter in Section 
6.3| Before that, in Section |6.2[ we study the law of the time of the first breakout, showing that 
it is approximately exponentially distributed. First of all, however, we start with the necessary 
definitions: 
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6.1 Definitions 

We will introduce several parameters which will be used during the rest of the paper. The two most 
important parameters are a and A, which are both large positive constants. The meaning of a is as 
in the previous sections: It is the right border of an interval in which the particles are staying most 
of the time. The parameter A has a more subtle meaning and controls the number of particles of 
the system and with it the intensity at which particles hit the point a. In Section [7j we will indeed 
choose the initial conditions such that Zq * e . 

When we study the system for large A and a, we first let a go to infinity, then A. Thus, the 
statement "For large A and a we have..." means: "There exist Aq and a function ao(A), both 
depending on the reproduction law q only, such that for A ^ Aq and a > ao(A) we have...". 
Likewise, the statement "As A and a go to infinity..." means "For all A there exists ciq(A), 
depending on the reproduction law q only, such that as A goes to infinity and a ^ ao(A). . . " . These 
phrases will become so common that in Sections [7] to [9] they will often be used implicitly, although 
they will always be explicitly stated in the theorems, propositions, lemmas etc. We further introduce 
the notation o(l), which stands for a (non-random) term that only depends on the reproduction 
law q and the parameters A, a, e, rj, y and £ and which goes to as A and a go to infinity. 

The remaining parameters we introduce are all going to depend on A, but not on a. First of all, 
there is the small parameter e, which controls the intensity of the breakouts. Indeed, when Zq m e , 
the mean time one has to wait for a breakout will be approximately proportional in e. Morally, one 
could choose e such that e ' 2 « e « .A -1 , but for technical reasons we will require that 

e < CJOA -17 , and (6.1) 

e > 0^- A,Q - (6.2) 

Another protagonist is r), which we will choose as small as we need and which will be used to bound 
the probability of very improbable events, as well as the contribution of the variable Y. It will be 
enough to require that 

V < e- 2A , (6.3) 



which, by (6.2), implies 

7] s; Ce 12 . (6.4) 



The last parameters are y and £, which are defined as in Lemma 4.2, with r] there being the r\ 
defined above. Note that the parameters r], y and £ appeared already in [7] and had the same 
meaning there. 



We can now proceed to the definition of the process. Recall the definition of fi in (5.1). As in 



Section 5.1, we denote by P x and E x the law and expectation of branching Brownian motion with 



drift —fi starting from a particle at the point x 6 R; we extend this definition to general initial 



distributions of particles according to a counting measure v. Recall from Section 3.1 that JV(t) 
denotes the set of individuals alive at time t. We want to absorb the particles at and do this 
formally by setting 

^V {t) = {ueyy{t):H {X u )>t}, 

where Hq is again the hitting time functional of 0. 

Instead of absorbing particles at a, we are now going to classify them into tiers as described 
in the introduction. Let u 6 U, t > 0. We define two sequences of random times (r n (u)) n ^-i and 
(cr n (u)) n ^o by T-i(-u) = 0, (Tq(u) = and for n $s 0: 

T n {u) = inf{s > a n {u) : X u {s) = a}, 
crn+x{u) = inf{s Ss T n (u) : X u (s) = a - y + (1 - fi)(s - r n (u))}, 
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Figure 2: Left picture: A graphical view of the tiers. The trajectory of one particle u is singled out (thick 
line) and the times ti(u) and ai(u) are shown. The tier particles are drawn with straight lines, the tier 1 
particles with dashed and the tier 2 particles with dotted lines. Right picture: The stopping line jV t (encircled 
particles). 



where we set inf ; 
stopping lines 



oo. See Figure [2] for a graphical description. We now define for t > the 



,(0 



'A 



(0 



{(u,s) e U x [0, t] : s = t\{u) and u £ J/${s)}, I > — 1, and 
{(u,s) £ U X [0,i] : s = <Ji{u) and u £ r yV Q (s)}, I > 0. 



(6.6) 
(6.7) 



That means, 

S) 



M 



contains the particles of tier I at the moment at which they touch the right 



barrier and ^ contains the particles of tier I at the moment at which they come back to the 
critical line. Note that the sets M\ and ^ are increasing in t and a/ - } ' - y '•' : 

1^0. We also set 



f 1} < y t {l) < M\ l> for every 



R' 



(i) 



m 



(i) 



In order to extend the definitions of the variables Zt and Yj to the current setup, we could simply 



replace -4^, a (t) by j9q(V) in their definition (see Section 5.1). However, it will be more useful to 



take special care of the individuals u for which ti{u) < t < ai + i(u) for some I ^ 0, since these are 
in some kind of "intermediary" state which is difficult to analyse. We therefore define the stopping 
lines 



^ ' = {(«, s) £ U x R + : u e jVq{s), r;_i(u) ^ s < ti(u), 



and either t < ai(u) = s or ai(u) < t = s}, I ^ 0, (6. 1 



and 



^t = {]A 



(0 



(6.9) 



feo 



,(/) 



In other words, the stopping line jY^ contains the particles of tier I that have already come back 
to the critical line at time t, as well as the descendants of those that haven't, at the moment at 



30 



which they hit the critical line. We then define for I > (recall the definitions of wz and wy from 
Section [5]), 

zf = 2 wz(X u (t)), Y®= 2 w Y (X u (t)), 



ue,_AC 



(0 



ue,yK 



m 



Furthermore, for any symbol S and ^ k ^ I < 00, we write, 

/ 



q(k;l) _ V 1 c(i) £('+) = gC; 00 ) 5 = 5(0+) 



i=k 



?<J) 



For a particle (u, s) e ^ , we now define the stopping line 

y{u,s) = ^ v ^ £ jj x R . y e ^( r ) ; ( u ,t) < (v,r) and r = a l+1 (y)}. 



This stopping line yields a collection (X„(r), r — s)r v 



(v,r )ey( u < s '> 



of space-time points and we denote by 



Z( u,s >, y( u,s ) and Wy' s the quantities from Lemma 5.8 corresponding to this collection of points 



r(u,s) 



(in particular, Wy U ' a) = ye y #.y^ u,s >). Of course, we have chosen the stopping line in such a way 



T (u,s) 



that the variable W y ' follows the same law as the variable W y defined in (4.2). We also define 
Tmax = m - ax ( V; r)ey( u ^i( r ~ s ) ■ We then define the "good" event 

Q(U,S) = ^y{u,s) ^ Q n ^{U,S) ( - U x ^J} 

n { sup #{w e ^(r) : («, s) < (v, r) < ^("> s )} < C} (6.10) 



and the event of a breakout 



B {u,s) = | Z (u,«) > £e A| (J ^ G (« 



,s)\c 



(6.11) 



(the inclusion of the "bad" event (G( u,s >) c is for technical reasons). If the event B^ u,s > occurs, the 
particle u is then also called the fugitive. We set 



PB = P a (B(0>°)), (6.12) 

and define the law of BBM started at a with the first particle conditioned not to break out: 

Q a (-) = P a (-|5 c ) = — ^ S 

I-PB 

where we set B = _B(0>°). We further set Z = Z^' ' and W y = W y ' and note that by Lemmas 



5.8 



and 4.2, we have for large a, 



P a (\Z - ttW\ > 2r?) + P a ((G (0 ' o) ) c ) < ry, 



where VF is defined as in (4.2). Hence, by ( |6.3[ ), (LL5J) and Lemma 4.2, we get 

1 



PB 



7T + 1 



ee' 



(6.13) 



(6.14) 



which goes to as A and a go to infinity, by (6.2). Furthermore, (6.13) yields for large A and a, 

Q a [Z] = (B[7rWl {7rW ^ e A {1+oW)+om ] + 0( V ee A ))(l + 0(p B )) 
= ir(A + loge + qeJ5\+ °( 1 ))' 
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by (4.6), (6.1), (6.3) and (6.14). In particular, we have for A^\ and large o, 



Q a [Z] ^ CA. 



Moreover, by (4.5), (6.13) and (6.3), we have for A^\ and large a, 

Q a [Z 2 ] ^ Cee A . 



(6.16) 



(6.17) 



Finally, note that by Lemma 5.8. we have Y ^ rjZ, Q a -almost surely, a fact that will often be used 
without further reference. 

We now define for every I e N the time of the first breakout of a particle of tier I, 



T® (w) = inf {t H:we J £ (u ' s) } , 



(0 



and set 



T (0;I) 



min T W , 



(u,s)e£%l 



with T = T(° ;0 °) = minT«. 

3>0 



(6.18) 



(6.19) 



Now fix t > and IeNu{oo}. We want to describe the system conditioned on T^ 0;l ' = t. 
For this, suppose that at time the particles are distributed according to a counting measure 
v = Xir=i <W We denote by % the fugitive of the breakout that happened at time T^ ' 1 ' and define 
Pi = P u (i < % | T^ 0; " = t). This yields a law (pi)f =1 on the initial particles, depending on v and 
t. Since the variable T^ 0; ", the time of the first breakout, is the minimum of the variables T^ , 
i = 1, . . . ,n, the times of the first breakout of the BBM descending from the particle i, we can 
decompose the process into 



Ylduj {i) T^ =t)=jftx P*'(du; w | T(° ; ') = t) x J^[P^(da; (i) | T^ > i). (6.20) 



i=l 



j=l 



J^« 



That is, we first choose according to the law (pt)iLi the initial particle that is going to cause the 
breakout. This particle spawns a BBM conditioned to break out at time t. The remaining particles 
spawn independent BBM conditioned not to break out before time t. 

Remark 6.1. Note that many results in this section can be done for BBM with varying drift given 
by a barrier function /. For example, with Lemma 5.7 one gets immediately that Py(T'°) > t) ^ 
P(T(°) > t) for all t > and it will be clear from the next section that in fact P/(T > t) ^ P(T > t) 
as well. However, in order to simplify notation and because we will not need the results often in 
this generality, we state them here only for BBM with fixed drift. 



6.2 The time of the first breakout 

For large A and a, define 



fcTL 



i „3 



20.4 



a . 



(6.21) 



We want to prove that the random variable T defined in the previous section is approximately 
exponentially distributed with parameter ps^Zo/a 3 , which is the statement of the following propo- 
sition: 



Proposition 6.2. Let < t < fe.211 - For A and a large enough, we have 

P(T > t) = exp ( - ttpbZq^ (l + 0(At/a 3 + p B )) + O(p B Y ) 



(6.22) 
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The proof proceeds by a sequence of lemmas. Lemma 6.3 gives a estimate on P(T"(°) > t). This 



is used in Lemma 6.4 in order to obtain an estimate on Q a (T > t), using a recursive argument. 



Finally, Proposition 6.2 is proven by combining Lemmas |6.3| and |6.4| 
Lemma 6.3. Let < t ^ a 3 . Suppose that p B ^ 1/2. Then, 

P(T(°) > t) = exp ( - 7rp B Z ^(l + 0(p B )) + O(p B Y )y (6.23) 

Proof. Let xi,...,x n be the positions of the initial particles. Since the initial particles spawn 
independent branching Brownian motions, we have 



P(T<°) >t) = Y[P X *(TW >t). 



(6.24) 



We have for every x e (0, a), 

px^ r (0) > f ) = E x 



[ n 



Lg(u,s) 



(M,s)G.« t (0) 



E J 



:i-pb 



i-^t 



(») 



(6.25) 



since by the stron g br anching property, the random variables Z^ u ' s > are independent conditioned on 



?l . By Lemma 



5.6 



and the assumption t ^ a , we have 



? (o) 



lE*^] - ™ z (x)t/a d | sj Cw Y {x), 
V x [(Ri 0) ) 2 ] ^ C(w z (x)t/a 3 + w Y {x)). 



By Jensen's inequality, (6.26) and the inequality | log(l — z)\ < z + z for z < 1/2, 



e*[(i-pb: 



a 



(o) 



^ exp ( - wp B w z (x)^(l + Cp B ) - Cp B w Y (x] 



(6.26) 
(6.27) 



(6.28) 



Furthermore, the inequality (1 — p) n < 1 — np + n(n — l)p /2 and Equations (6.26) and (6.27) give 



E : ' 



vI-Pb) 



j(0) 



^ 1 - 7rp B -u;^(x)^7(l - Cp B ) + Cp B w Y {x). 



(6.29) 



The lemma now follows from (6.25), (6.28) and (6.29) together with the inequality 1 + z =$ e z for 
zeR. □ 

In the following lemma, note that according to the definition of the tiers, a particle starting at 
a starts immediately in tier 1. 

Lemma 6.4. Let < t ^ ^ |6.2il - Then, for large A and a, 

Q a (T>t)^exp^-Cp B A(^+r,)). 
Proof. We have 



(6.30) 



Q a (T>t) = Q a \ [7 P Xu(s) (T>t-s)j ^Q a [P u {T> t)], 
where v = £\ x o>(i) $x u (s)- Since T > t implies T( ) > t, we have 



(6.31) 



P U (T >t) = P U (T > t\T<® > t)P u (T^ > t). 



(6.32) 
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Let Z = Z(0'°> and Y = y(0>°), such that Y < rj Z, Q a -almost surely, by Lemma 
definition of the "bad" event B^™' '. By Lemma 6.3, we have for large A, 



5.8 



P"(T(°> > t) > exp ( - Cp B Z (4j + r?) 
Furthermore, with the notation from Section [5.5[ with p(s) = pb, 



and the 



(6.33) 



\T> t\T® >t) =P U (T> t) =P U \ j [ Q a (T> t-s)] ^P u \Q a (T> t) R t 0) 

I.. „^,®(o) 



{u,s)eM\ 



By Jensen's inequality and Lemmas 5.13 and 5.6 this implies 



(«) n 



V V {T > t\T<® >t)> Q a (T > t) EV l R < ] 55 Q a (T > £)^('/" 3 + 



o(i)) 



(6.34) 



Equations (6.31), (6.32), ( |6.33[ ) and (6.34), together with Jensen's inequality and (6.15), now yield 
for large A and a, 



Q a (T >t)> Q a (T > t) 



w 2 A(t/a s +0(r))) 



x exp -Cp B A(A + ri) 



(6.35) 



By the hypothesis on t and (6.3), the exponent of Q a (T > t) in ( 6.35[ ) is smaller than 1/2 for large 
A and a. This yields the statement. □ 



Proof of Proposition 6.2, The upper bound follows from Lemma 6.3 and the trivial inequality T ^ 
T^> . For the lower bound, we note that as in the proof of Lemma 6.4 , we have by Jensen's inequality 
and Lemma |5.13[ 



?C°)i 



p(T > t) = P(r > 1 1 r<°) > t)p(r(°) >t)> Q a (r > tfW J ]p(r(°) > t). (6.36) 

By Lemmas |6.4| and |5.6[ we have 



A°h 



Qa{T > t) m\ ] > exp - Cp B A(£ + V )(^Z + Y 



(6.37) 

The lower bound in ( J6.22 ) now follows from (6.36), (6.37) and Lemma 6.3, together with the 
hypothesis on t, (6.3) and (6.14). □ 

Lemma 6.5. Define 7 = (irpsZo) -1 . Suppose that Yq ^ C and let a ^ and n e N. Then, for 
large A, for every I e N u {00}, 



E[(TW) /o 3 + Qri(TMsaS)HC 2_»^T)»" 



fc=0 



k\ 



Furthermore, if > 5 = A 1 o(l), i/ien /or large A and a, 

E[(T/a 3 )l (T ^ a3) ] = 7 (1 + 0(A5 + p B )) + 0((S + 7 )e"°^) 
Proof. We first note that we have, for n ^ 1 and p > 0, 



(t + Q) n - 1 e-" 4 dt= V 

J fc=0 



n 



l)!a fc 



fe!p 



,n—k 



(6.38) 



(6.39) 



(6.40) 



Now, we have 



E[(T (0;0 /a 3 + a )«i (T(0) ^ 3) ] = [V + a)" P(T(°^/a 3 e di) 

< n f (t + a)™" 1 P(T(° ; ') > to 3 ) dt 
Jo 



a 
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The inequality (6.38) now follows from Lemma 6.3 and (6.40) and the hypothesis on Yq. For the 



second part, we note that 

r-6 

E[(T/a 3 )l (T ^ a3) ] = P(T > to 3 ) dt - 8P(T > 5a 3 ), 
Jo 

and by Proposition |6.2| and the hypothesis on Yq, we have for t ^ 6 and large A and a, 
P(T > to 3 ) = (1 + 0(p B )) exp(- 7 - 1 t(l + 0{A5 + PB ))). 



Equation (6.39) now follows from the last two equations. 



□ 



We now show how we can couple the variable T with an exponentially distributed variable: 



Lemma 6.6. Suppose that e~ Zq = 1 + 0{e 3 ' 2 ) and that Yq =5 t]Zq. Then there exists a coupling 
(T,V), such that V is a random variable which is exponentially distributed with parameter TrpBe A , 
T is a (V) -measurable and P(|T/a 3 — V\ > e 3 ' 2 ) ^ Ce 2 for large A and a. 

Proof. For brevity, set p := 7rp B e . Let F be the tail distribution function of T, i.e. F(t) : = P(T > 
t). The number of individuals in a BBM tree being at most countable, T has no atoms except oo. 
We can therefore define a random variable U which is uniformly distributed on (0, 1) by setting 

U = F(T)l (r<00) + U'F(cc)l {T=Q0 ), 

where V is a uniformly distributed random variable on (0,1), independent of T. Now define 
V = — p~ 1 logU. Then V is exponentially distributed with parameter p and T = F^ 1 (e^ p ), 
where F^ 1 denotes the generalised right-continuous inverse of F. Hence, T is cr(V)-measurable. On 
{T < oo}, we have by Proposition 6.2 for a large enough, 

V = -p-\Kp B e A e- A Z Q T/a\l + 0(AT/a 3 + p B ))) + O(p B Y ) 

= T/a 3 {l + 0{e 3 / 2 + AT/a 3 +p B )) + 0{p B e. A r]), 



(6.41) 



by the hypotheses on Zq and Yq. Hence, by (6.2), (6.4) and (6.14), we have for a large enough, 

\T/a 3 -V\ = 0(e 3 / 2 T/a 3 + A(T/a 3 ) 2 ) + 0(e 2 ). 



But now we have by Lemma |6.3| for large A and a, 

P(T/a 3 > J/ 8 ) ^ P(T(°)/a 3 > e 7/8 ) ^ Ce- 0( -^ ^ e 2 /2. 
The statement now follows from (6.41[) and (J6.42 ) together with (6.1) 



6.3 The particles that do not participate in the breakout 



(6.42) 

□ 



In this section, we fix r < l\6.2i[ We define the law of BBM conditioned not to break out before r 
in the tiers 0, . . . , I by 

P,(-) =P(.|T(° ;Z ) >t), With P = Poo. 

Expectation w.r.t. P; i s den oted by Ej. Under the law P;, the process stopped at Jz?# then follows 
the law P from Section 5.5 



with 
pit) =p B l(t^ T) 



(l-p B )Q-(T^^r-t), 



(6.43) 



such that by Lemma 6.4, (6.3), (6.16) and the trivial inequality T^ 0; > ^ T, we have for large A and 
a, 

blloo < Cp B . (6.44) 



In particular, by Lemma 5.12, (5.2) and (6.44), we have for large A and a, 



{h{x,0))- l ^l + Cp B . 



(6.45) 



Lemmas 5.9 and 5.11, together with (6.44) and (6.45) now immediately imply the following: 
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Corollary 6.7. Let xe(0,o), O^Kt^ fe.2il and g : [0,a] — > R + be measurable with g(0) 
g(a) = 0. Define S t \ — 2 ^-(o) g(X u (i)) or S\ = R t . Then, with E; as above, 



Ef K (0) ] = (1 + O(p B ))E-K (0) ] and Ef[(5 t (0) ) 2 ] < (1 + 0(p B ))E*[(4°>) 2 ]. 
Moreover, as in the proof of Lemma |5.13[ one can show that 

Qf[Z] = (1 + 0( PB ))Q a [Z] and Qf[Z 2 ] < (1 + 0(p B ))Q a [Z 2 ]. 
We define two nitrations (^)/>o and {J^Oi^q by 



(6.46) 



such that Sf/ c j£f c ^+i for every /. Now define for measurable la [0, r], 



Z 



(J) 



V 7(«»*) vW 



S 



y(u,t) 



(u,t)e3?} l 'n(UxI) 



(M,i)e,y; (0 n((7x/) 



with Z 



(0 



Z 



(0 

0,[O,t] 



and Y, 



(0 



r (J) 



Y~'r Q i . Furthermore, recall from Section 



ZL and ^0 and the corresponding quantities for Y. 

Lemma 6.8. We have for all =% k < I, x e (0, a) and large j4 and a, 



6.1 



the definition of Z 



(k;l) 



%[z^ i ^ ^ (tt + cp B )Q a ^](^4 fc) + cad**; 



In particular, for large A and a, 

%[Z% +1) ] ^ CA(^Z + CY ) (vr 2 ^ + q^gF?))". 
In £/ie case fc = 0, we also have for large A and a, 

%[Z$] > (tt - C PB )Q a [Z](^o - C%). 
Proof. We have for ^ fc ^ I, 



Vi[Z, 



(fc+i) 



|J^] = E,[ 2 Z< u - 



t) 






Qt-dz]Ri k \ 



(6.47) 
(6.48) 
(6.49) 
(6.50) 



since conditioned on J^., the random variables Z^ u,t > are iid under P/ of the same law as Z under 
Qf_ k and independent of Mk, by the strong branching property. Now, we have 

E Z [R^ | Sf fc ] = 2 EX " (t) [4-t I T > ^ " *] . 

(u,i)e t y T (fe) 

such that by Lemma |5.6| and Corollary |6.7[ 

E t [RW | gy < (i + ( PB )) (vr Jzg> + ^teSSrl } ) • (6-51) 



Equations (6.50) and (6.51|) together with (6.46) give (6.47). Equation (6.48) follows easily from 



<k) 



,(k) 



(6.47) by (6.16) and the fact that YX, =% r\Z\J , P;-almost surely for 1 < k < I. Now, in the case 



k = 0, we have % = ^o by definition. Denote the positions of the initial particles by xi, ... ,x n . 
By Lemmas pT6l|5.13| and ( |6.44[ ), 

n n 

e,[4°)] = J]e*[4 0) ] ^ 2 EXl [^ 0) ] - bllooE Xi [(4 0) ) 2 ] 



1=1 



1=1 
r 



5s vt^Zq - Qo#0 - llpllooQOTltia^O + >o) , 



which yields (6.49). 



D 
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Lemma 6.9. Suppose that t < r. Then for large A and a, we have 

v2 



B[Z t ] = Z (l + TrQ a [Z]i. + 0( PB + (A%) 

If moreover t ^ 2a 2 , then 

C 



O{AY ), E[Y t ] ^C(Y + nA^Z ) . 



B[Y t l {R[ta2tr0) ] sc _(Z + AY ) and ?{R [t ^ A ± 0) < C V (Y + Z Q ). 



Proof. First note that we have mf xe r 0a -ih(x,0) ^ 1/2 for large A, by (6.45). Furthermore, by 
the hypothesis r < 46.211 and (6.3), we have 7r 2 A( J ^ + Crj) < 1/2 for large ^4 and a. The first 
two inequalities now follow from Proposition 5.2 Corollary 6.7 Lemma 6.8, (6.46) and (6.3) by 
summing over k. In particular, if t ^ 2a 2 , then for large ^4 and a, E[Z t _ a 2] < C (Zq + Alo) and 
E[y t _ a2 ] sS Crj(Z + Y ). Together with Proposition 
the other two inequalities. 



5.2 



Lemma 



5.6 



and Corollary 6.7, this proves 

□ 



In order to estimate second moments, we will make use of the following extension to the Many- 
to-two lemma (Lemm a |3.5[ ). For x, z £ (0, a) and < t < r, we define fn 2 (x,t) to be the quantity 
m 2 (x,t) from Section 5.5| corresponding to the penalisation from (6.43) and p t (x,z) to be the 
density of tier I particles at position z and time t under the law P x , not counting the particles u 
with t < o~i{u). Then set Pt = pj • 

Lemma 6.10. Let w : [0, a] — > R + be measurable and define St = Yi( u t)eJi w (X u (t))- Then 

E*[S 2 ]=E X [ V w(X u (s)) 2 ] + /3 n a m 2 (z, s )p s (x,z)E^[S t ] 2 dzds 

,.. ~Z „ Jo Jo 



(u,s)e,yY t 



J] E^ x nM^)[S t ]E^2M^[S t ] . (6.52) 



Proof For (vi, si), (v 2 , s 2 ) £ ^ we write (vi, si) a (t; 2 , S2) for their most recent common ancestor. 
Then define for I ^ 1, 

•<*i = {((n,si), (v 2 , s 2 )) £ ^ 2 : Vi + v 2 and if (vi,s{) a (v 2 , s 2 ) = (v , s ), 

then rj_i(vo) < so < ^(v )}, 
and for / 5= 0, 

^2 = {((«i, s l)> ( u 2, s 2 )) £ ^ 2 :»i^2 and if (vi,s{) a (v 2 , s 2 ) = (v , s ), 

then erj(v ) < s < n(vo)}, 
We then have 



E J 






((«i,«i),(«2,sa))e^i 



^ J] E( x "i( Sl )' Sl )[5 t ]E^2( S2 )' S2 )[5i], (6.53) 



and by Lemma |3.5| (see also Remark 3.4), 



E J 



2 



w(X 1 , 1 (si))w(X„ 2 (s 2 )) 



((«l,»l),(«2,S2))e^5 



(0 



A 



'0 f f m 2 (2,s)p«(x,z)E(^[S t ] 2 dzds. (6.54) 
Jo Jo 



Equations (6.53) and (6.54) and the tower property of conditional expectation yield the lemma. □ 
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Remark 6.11. An analogous result holds for Rt. 
Lemma 6.12. We have for every t < t, 

t 



Var(Z t ) < Cee A (-^Z + Y ) and Var(i2j) < Cee A ^Z + Y ) . 



Proof. By Lemma 6.9 and the hypothesis t < r ^ l\6.2i\ , we have for every x e (0, a) and s ^ t, 



B^ x ^[Z t ] ^CA 






^CAe-'^ a - x) . 



(6.55) 



Now, as in the proof of (5.16), we have by (6.55) and Corollary 6.7 



which yields 



f [\f ) {x,z)% { - z > s) [Z t \ 2 dzds^CA 2 (—wz{x) + w Y {x)), (6.56) 

Jo Jo ^ a ' 



f \ m 2 (z, s)p s (x, z)^ z ' s \Z t f dz ds 
Jo Jo 

C&\ ^ f \ a p^s (Xu(so),z)B^[Z t fdzds 



sc 



< 



(u,s )ey t ' 
t 



(6.57) 



CA^l- E &[z 0m ] + &\y 0m ] 



^CA 2 (^wz{x)+w Y {x)), 



by Lemma 6.8 and the hypothesis t < r < ^ J6.2il - Furthermore, by (6.55), (6.3) and Lemma 5.8 
have 

E*|" V V % l {Xv 1 {si),s 1 )Y Zt ^ i {x V2 {s 2 ),s 2 )Yz t ^ 

(u,s)eMt (v 1 ,s 1 ),(v2,s 2 )e,y( u - a \v 1 ^v2 



we 



^CQ a [Z 2 ]B x [R t ]^Cee A (^wz{x)+w Y (x)), (6.58) 



as Lemma 



6.9 



by (6.17) and Lemmas 5.6 and Corollary 6.7 Lemma 6.10 together with (6.57) and (6.58) as well 



and the inequality w z < Cwy gives 



Var (Z t ) ^ E x [(Z t ) 2 ] sc Cee A (^w z {x) + w Y (x) 



(6.59) 



Summing over the initial particles yields the inequality for Var(Zt). The proof of the second in- 
equality is analogous, relying on Lemma |5.6| and Remark 6,11| D 



We finish the section by a corollary which will be useful in the next section. 

Corollary 6.13. Suppose that t ^ r and x € (0, a). Then for large A and a, we have E z [Zt] ^ 
CAe -(a-x)/2 ) &[ Z f] sc Cee A e- ( - a - x ^ 2 and W[Y t ] s= Ce^'^ 2 . 



Proof. Immediate from Lemmas 6.9 and 6.12 and (5.2). 



□ 
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6.4 The fugitive and its family 



We now describe the BBM conditioned to break out at a given time. Recall that 9/ denotes the 
fugitive. For simplicity, we write 77 and 07 for 77 (^) and ai(fy), respectively, I e N. On the event 
T = r"', we define %j to be the ancestor of ^ alive at the time o~j, j = 0, . . . , I. By the strong 
branching property, we have the following decomposition: 

Lemma 6.14. Ie«fceN,IeNu {00} with l> k andt^O. Conditioned on & 0) T( 0;Z ) = T^ = t, 
To and %), the BBM admits the following recursive decomposition: 

1. The initial particles u ¥= &o spawn independent BBM conditioned on T^ 0; > > t, 

2. independently, the particle ^b spawns BBM conditioned on a particle (call it %q) hitting a 

for the first time at the time tq, all the children of which born before To being conditioned on 
T (0;0 >t 

(a) If k = 0, this particle %q spawns BBM conditioned on the event B of a breakout. 

(b) If k > 0, it spawns BBM starting at the space-time point (a, tq) conditioned on T' 0; ^ = 
T^ ' = t. In particular, if we write 5? = ^^ o ,T °' , then conditioned on &y, the particles 
in 5^ spawn BBM starting from the collection of space-time points S? , conditioned on 

T (0;/-l) = T (k-1) =t 

Note that in the case 2b above, the subtree spawned by (^q, to) follows the law Q conditioned on 
T(°' 1 ) = jW = t, hence the law of ^W' r °) is not the same as under Q. In particular, E[Z(' ?/ o> T o)] <£ 
CA. Indeed, conditioning on one of its descendants breaking out at a later time corresponds to a 
kind of size-bias on the number of particles. However, it is still true that Z^o^o) < ee , by the 
definition o f a br eakout. 

Lemma 6.14 gives a decomposition of the BBM conditioned on T^ k > into k + 1 pieces. In order 
to describe what happens in a single piece, define &<& to be the a-field generated by the family of 
events {{T^ = T^}) keN and by &, {X^{s)) a ^ s<Tj , (d u , k u ) u ^ and «y(^), for j = 0, . . . , k, 
on the event {T^ ^ = T^}. Note that in particular, on the event {T( 0; ') = T( fc )}, tj, {^j,o-j) are 
&<% -measurable, j = 0, ...,k. It is plain that conditioned on J^V, the subtrees spawned by the 
children of the ancestors of % during the intervals [o~j,Tj) are independent BBM conditioned not 
to break out before t. 

It remains to describe the trajectory of the fugitive and its reproduction. We could do by a 
decomposition of the process at the first branching, but it is faster to use the Many-to-one lemma 
instead, which is the method of proof of the following lemma, which we state for general penalisations 
p(s). This result is essentially \24\ Theorem 1]. 

Lemma 6.15. Let t ^ 0, x 6 (0, a) and p : R + — > [0, 1] be measu rable with p(s) = for s large 



enough. Denote by P the law associated to p(s) as in Section 5.5, Recall the definition of e(x,t 



from (5.40). Then, given a family of J^co) -measurable non-negative random variables {Y u ) ue jj, we 

"*t 

have 



W 



Y, 



3u e U : (u,t) 6 



?(0) 



E* 



Y f 



•J e(6,s)ds 



E* 



6 e 



H (0 > H a (£) 
#o(0 > H a (0 = 



t 



where under ~p*' x } the spine moves according to standard Brownian motion and spawns particles 
with rate fh~x(x, t)/3oQ(x, s) according to the reproduction law ((rni(x, t)) kq(x, t, k))k^o, which start 
independent P-BBM. In particular, conditioned on (X u (s))o^ s ^t o-nd (d v , k v ) v < u , the children of the 
ancestors of u follow independent P-BBM. 

Proof. We have 



B x 



Y u 



Sue U : (u,i) e 



M 



E d 



^(u,s)e,4 0) 1 (H a (X u )edt)Yu 



E-< 



2( U , S ) 6 «(0) \H a {X u )edt) 



(6.60) 
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and by Lemma |3.1[ 

E»[ 2 

(«,s)e.* t (0) 



1 (H a (X u )edt) Y u 



E* 



Ye e S" i (? 3 ' s )/ 3o( 3(€ s ' s ) d ' s l 



(H (0>H a (Oedt) 



(6.61) 



According to the description of the conditioned process in Section 5.5 and the description of the 
spine in Section |3.3| the particles on the spine follow the Doob transform of Brownian motion with 
drift — \x by the harmonic function h(x, t) and spawn particles with rate rni(x, t)/3oQ(x, s) according 



to the reproduction law ((mi(x,t)) kq(x,t,k))k^o- With Girsanov's transform, (6.61) yields, 



E- 



1 2 

{u,s)eM 



l (H a (X u )edt) 



Y u 



,-Ai(o-!E) MM) 

h(x,0) 



E 



*,x 



Y it€ jo l {Ho{i : )>Ha 



(Oedt) 



(6.62) 



Equations (6.60) and (6.62) yield the first statement. The second statement follows by taking 
appropriate test functionals (Y u ) ue jj. □ 



Corollary 6.16. In addition to the assumptions in Lemma 6.15, suppose t e [0,o ] and \p\co ^ c p , 

(6.63) 



where c p is some universal constant implicitly defined below. Then, 



E x \Y u hue U : {u,t) e 



,(o) 



^CE* 



Yf, 



H (0 > H a (£) e dt 



Proof. By Lemma |2.2[ the inequality e x ^ 1 — x and the hypothesis on t, we have 



E* 



,-^0 e (is,s)ds 



Ho(0 > H a (0 =t]= W t x ^ [e-U^,s)d s] >{l _ CM 



The statement now follows from (6.64) and Lemma 6.15 



(6.64) 

□ 



We come back to the BBM conditioned to break out at a given time and set up the important 
definitions. Recall that % denotes the fugitive. We will denote by a bar the quantities referring 
to the particles spawned between the times cri(^/) and ti(&) for some I 5= 0, and by a check those 
referring to the particles spawned between the times ti(%) and gi+\(%). Formally, we set 



1={UE 



Jlft : (u,t) A(W,T)eUx \J[ai(^),n(^))h 



feo 



{(u,s) e U< } : (u,s) A(W,T)eUx {J[at(W) , n (W)))} 



2>0 



feO 



and 



We then define 



jY t = {ueAi: («,t) a {<&,T) eUx (J[^_i(«r), ai(W))}. 
Z t = ^ w z {X u (t)), Y t = J] w Y (X u (t)), Rt = #M t , 



ue._At 



ue,yVt 



and 



Zt = J] w z (X u (t)), Y t = J] «*(*«(*)). 

Note that on the event T = T^- ' , we have ^¥{T) = by definition. 



By Corollary 6.13, the following J^-measurable functional will be of use in the study of the 
bar-quantities. 

<% = S 2 l( ff «<du<ro(*« ~ l)e-(-^^))/ 2 , on {T = r»}. (6.65) 

where we recall that d u and fc u are respectively the time of death and the number of children of the 
individual u. 
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Lemma 6.17. For large A and a, we have for every I e N, 



Proof. By Lemma 6.14 and Corollary 6.16 (which can be applied because of (6.44)), we have for 
every Z e N, 

i 



E[<% | T = T® ^ fed < C S E[Wi[ r ^ e^ 8 "* 

4 =o L LJ o 



<)/ 2 di 



7 1 = X"(0 < 



where W% is the law of a Brownian bridge of length n — o"j from X^(o~j) to X^-(tj). The statement 
now follows readily from Lemma |2.2| □ 



We can now study the probability that the fugitive stems from a given tier. 

Lemma 6.18. Suppose that C\e < Zq < C^e and Y$ =% ??^o- 27ien /or large A and a, we have 
for le {1,2}, 



p( T G+) <r (0iJ-i)) <^C(eA) 1 , 
Proof. By (6.2), (6.14) and ( |6.46[ ), we have for large A and a, 

P(T(°) > fell) ^ exp(-C,4/e) < exp(-C/Vi). 



(6.66) 



Now, for the rest of the proof, let t < fe,2il and let z^ = 2X=i <W We have by the decomposition 
(6.20[) of the process conditioned on T' 0; ' _1 ) = t, 



p ^ T 0+) > t | T (0;J-i) = f ) = £ w p"-*» < (t( j+ ) > j|t(0;1-i) > t) x p**(T^ > 1 1 3-(o^-i) = t ). 



i=l 



(6.67) 

Define q = 7r 2 ^4(t/a 3 + C J6.48K ?)i which is less than 1/2 for large A and a by the hypothesis on t. By 
Lemma 6.8 and the hypothesis, we have for every j < Z, 



-t ! 



and moreover for every (x,s) e (0, a) x [0, £], 

\(x,s) vy{l-jh 






For every k < Z — 1, (6.68) applied to the particles in ^ , j = 1, . . . , k yields 
E[Z§ | T« = T^'" 1 ) = t] < C £ e A 2 4" j ^ Cee^" 



l-k 



i=i 



Moreover, by (6.69), Corollary 6.13 and Lemma 6.17, we have 



e[z<? | r^ = t^-v = t] < c^4- fe . 



In total, we get by (6.68), (6.69) and (6.70), the hypothesis on Zq and (6.2), 



E[zf + Z<J + Zg | T« = r^'" 1 ) = t] < Ce A f4 + ec\- k 



Jensen's inequality and Proposition 6.2 together with (6.72) and the inequality 1 — e x < x give 
p( T (J+) < 1 1 T (k) = r (o ; /-i) = t j ^ Ccfe'^ + cj- fc ). 



Summing (6.73) over k, the law of total expectation gives 

p( T (i+) < t | r (0;I-l) = t ) ^ Cq^- 1 ^ + C t ). 



(6.68) 
(6.69) 

(6.70) 

(6.71) 

(6.72) 

^ive 
(6.73) 

(6.74) 



The lemma now follows by integrating (6.74) over t from to l\6.2i\ and using Lemma 6.5 and 
([6~66]). □ 
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Remark 6.19. One may wonder whether one can simply calculate P(T"' e dt | T^ 1 ^ 1 ' > t) for every 
I ^ 1 and t ^ 46.2il > using only the tools from Section |6.3[ This would require fine estimates on the 
density of the point process formed by the particles from tier I — 1 hitting a just before t. These 
estimates can be most easily obtained if one stops descendants of the particles hitting a at a (large) 



fixed time C, instead of the line from Lemma 5.8, with which the results in this paper would hold 



as well. However, in order not to lose generality we stick to Lemma 6.18 which is enough for our 
purposes. 

7 The B-BBM 

We will now define properly the BBM with the moving barrier, also called the B-BBM (the "B" 
stands for "barrier"), which will be used in the subsequent sections to approximate the iV-BBM. We 



will still use all the definitions from Section 6.1 with one notational change: Recall that by (6.20), 
we can decompose the process into two parts; the first part consisting of the particles spawned by 
the ancestors of the fugitive and the second part consisting of the remaining particles. As in Section 
6.4 the quantities which refer to the particles of the first part will be denoted by a bar (e.g. Z) or 



check (e.g. Z). The quantities of the se cond part will be denoted with a hat in this section (e.g. Z), 



in reference to the law P from Section 6.3 Moreover, the quantities related to the descendants of 



the fugitive will be denoted by the superscript "fug" (e.g. Z fug ). 

7.1 Definition of the model 

Suppose that we are given a family (f x )x^o of non-decreasing functions f x e ^ 2 (R, R+), such that 
for each x ^ 0, f x (t) = for t ^ 0, f x (+cc) = x and for each 5 > small enough there exist 
%o = xo(5), to = to(5), such that 

1. xq{5) — > oo as 5 —*■ 0, 

2. \\fx\\ sSfT 1 for all x e [0,x ], 

3. f x (t) >(1- S)x for all t > t Q and 

4. t < S- 1 . 



where ||/|| is defined in (5.4). A family of functions with the above properties exists, for example 



the following, which we will choose in Sections [8] and [9j 

f x (t) = log (1 + (e*-l)9(t))- (7-1) 

For every A and a, let v$ = u ' be a (possibly random) finite counting measure. We now define 
the B-BBM with initial configuration of particles uq. Starting from BBM with constant drift — jjl 
with initial configuration uq, we define for each n 6 N a stopping time n and for each n 6 N* the 
barrier process (X\ ■ )t G [e n _ li e n ] as follows: 

1. We set e = 0. 

2. Denote by T = T\ th e tim e of the first breakout of the BBM absorbed at and by % the 
fugitive, as in Section 6.1 We set X\ = for t e [0,T]. 



3. Define T~ = (T - e A aF)/\ and 

A = log (e- A {Z T - + Z 0vT - + Z^^)) v 0. (7.2) 

where Z 0vT - = Z T _ ]£T~ ^ tq(^) and = Z)J otherwise. The above quantities are defined 



in Section 6.4, after Corollary 6.16[ 



Define T+ - T x + = T + t^L^ ] and 6i = (T + e A a 2 ) v T+ . Note that T+ and therefore also 
01 is a stopping time for the BBM. Now define 

xl 1] = f A (^=Jly te[T+j@ll 

42 



We then give to the particles an additional drift — {d/dt)X\ for t e [T + , Gi], in the meaning 
of Section 15.11 

5. We have now defined T±, T^ , Gi and X^K We further define v\ to be the measure formed by 
the particles at time Gi, which have never hit 0. To define T2, Tg", Q2 and XL 2 !, we repeat 
the above steps with the process formed by the BBM started from those particles, with the 
definitions changed such that the barrier process starts at X^ = Xk, time starts at Oi etc. 

6. We now construct the barrier process X\ from the pieces by X\ = X\ , if t e [0 n _i, n ]. 
Define the event Go = {suppz^o c (0, a), \Zq — e \ =% e 3 ' 2 e , Yq ^ n}. Recall the definition of 



the phrase "As A and a go to infinity" from Section 6.1 Define the predicate 
(HB) (The law of) vq is such that P(Go) ->1 as A and a go to infinity. 
(HBo) vq is deterministic and such that Go is verified. 
Furthermore, if for some n > 0, {t ■ ,a )?=i e R+ for all ^4 and a, then define the predicate 

(Ht) There exists < t\ < ■ ■ ■ < t n , such that for all j e {1, . . . , n}, t • ,a — > tj as A and a go to 

infinity, with t 1 ,a = if £1 = 0. 

Theorem 7.1. Suppose (HB). Let n > 1 and (£• ,a )?? =1 6 R™ suc/t tfoai (77£,) is verified. Define the 
process 

(X t ) t>0 = (X [ $-K 2 At) t>0 . 

Then, as A and a go to infinity, the law of the vector (Ia,«)L converges to the law of (Lt-) r - =1 , 



where (Lt)t^o is the Levy process from Theorem\l.l\ 



A stronger convergence than convergence in the sense of finite-dimensional distributions is con- 
vergence in law with respect to Skorokhod's (Ji-)topology (see J28, Chapter 3]). Obviously, the 



convergence in Theorem 7.1 does not hold in this stronger sense, because the barrier is continuous 
but the Levy process is not and the set of continuous functions is closed in Skorokhod's Ji-topolog}Q 
However, if we create artificial jumps, we can rectify this: 

Theorem 7.2. Suppose (HB). Define Jt = Xq , if t e [6„,0 n +i), for n e N. Then as A and a 
go to infinity, the process (X[)t^o = (J a 3 t ~ ^At)^^ converges in law with respect to Skorokhod's 



topology to the Levy process defined in the statement of Theorem 7.1 



For each n ^ 1, we define the event G n to be the intersection of G n _i with the following events 



(see Section 6.1 for the definition of JVt): 

• supply c (0,a), 

• jY 6n c(Jx {G n } and n > T+ (for n > 0), 

• \e- A Z 0n - 1| < e 3 / 2 and Y @n sj n. 
Note that G n e ^e„ for every n 5= 0. 



The core of the proof of Theorems 7. 1 and |7.2| will be the following proposition: 

Proposition 7.3. Fix A 6 R and define 70 = (71p.Be). There exists a numerical constant b > 0, 
such that for large A and a, we have P(G n | J^o)1g ^ ne 1+b for every n ^ and 



logE 



e e « 



Lg = (n 7 o(«(A) + iAvr 2 ^ + o(l) + 0(e fe ))) 1 G(J . (7.3) 



where k(X) is the cumulant from (1.1) and where o(l) and O(e) may depend on A. 



One could prove convergence in the Mi-topology (cf. [58, Section 3.3]), but the use of this less standard topology 
would lead us away from our real goal: the proof of Theorem 11.11 
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7.2 Proof of Proposition |7.3 



First note that conditioning on #0j we can and will assume without loss of generality that the 
hypothesis (HBo) holds and will only state it again in the lemmas. In the same vein, we will always 
assume that A and a are large and will omit stating this fact, except again in the statements of the 
lemmas. 

Furthermore, it is only necessary to treat the case n = 1, since the rest can be obtained by 
induction: Suppose that we have shown the result for n = 1. Since the process starts afresh at the 
stopping time n , we then have P(G n \G n+ i | J^e„) ^ e 1+b lG n for every n ^ 0, which implies the 
statement on the probability of G n . Similarly, if we set k'(X) = 70 (k(A) + i\ir 2 A + o(l)), where 
o(l) is the term from (7.3), then 



E[e 4A( < ] +1 -< ]) l Gn+1 I ^ J = (e^ + °^ b ) + Q (P(G„\G n+1 | JT J) ) l Gn 



l + 0(e 1+b ))e K '^ + °^ b h Gn , 



because 70 = e(l + o(l)) by (6.14). It follows that 



E 



ixxy 



"Wl 1 r < 

-"-Crn + l 



e K'(X) + 0(e b ) E 



e e 



[00] 



1 + &N 



1 + 0{e^°) l Gn 



e K'(\)+0(e b ) E 



ixx^r' -, 
e e„ ± Gn 



By induction, and the fact that P(G,^ +1 ) ^ ne 1+b ^ e°' e ', we then obtain (7.3). 



Throughout the proof, we will work on several "good sets" , which we all define here for easy ref- 
erence, since they will be reused in the following sections. For this reason, some of them (for e xam ple 



7.3 



G) are actually more restrictive than what would be necessary for the proof of Proposition 

• G<% = {e A a 2 sj T sj ^a 3 , T = T^ 1 ), <% sc e A ' 3 }. 

• G fug = {Z^?) sc e A /e, Y^^ sc zW/y) n G^^ (the event G^> T ) was defined in 

(6T0l))- 

• G = {\Z T - - e A \ sc e l ' A e A , Y T - sj a - l e 3A ' 2 }. 

• G = {Z 0vT - sc e 1/4 e A , Y 0vT - ^ e~ A/2 } (we define Y 0vT - analogously to Z 0vT -). 

• Ga = G>% n Gf ug n G n G. 

• G n bab ("nbab" stands for "no breakout after breakout") is the event that no bar-particle 
breaks out between T and ©i + e A a 2 and that no hat-, fug- or check-particle (fug-particles = 
descendants of ,y^> T )) hits a between T~ a nd 61 + e A a 2 . 



Recall the definition of &<% from Section 6.4 and set J^a := -^ a u v &t- v ^t-i such that 



Ga e =^a and the random variable A is measurable with respect to #a- Define further Z& 
Z T - + Z^' T ^ + Z 0vT - and 1a analogously. Finally, define P n bab = P(- 1 G nba b)- 



The probability of Ga- 

Lemma 7.4. Suppose (HBo). F° r large A and a, we have P(Ga) ^ 1 — Ge 5 ' 4 . 

Proof. By Proposition |6.2| and Lemma |6.18 we have 



P(e A a 2 s$ T s; 



<ea 



T = T^ l h>l-CA 2 e 2 . 



(7.4) 



Since \fe =% fe.2i| for large A by (6.1 ), Equation (7.4) and Markov's inequality applied to Lemma 
yield 



6.17 



2_2 



P(G^) > l-CA z e 



Furthermore, by (6.13) and the last statement of Lemma |5.8[ we have, 

P(GL g ) < Pb 1 ( P ° (* W > ?> A le - 2 V ) + n) ^ Ce 2 , 



(7.5) 
(7.6) 
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by d45l), ((6^1), Kty and (|6~14]). 



In order to estimate the probability of G, we will calculate first- and second moments of the 
quantities in the definition of G. These estimates will be needed again later on, when we calculate 
the Fourier transform of the variable A. 



By Lemma 6.9 and the hypotheses, together with (6.14), (6.2) and ( |6.3[ ), 
E[e- A Z T - | J%]1 G ^ = 1 + 7rQ°[Z]£ + O (e 3 / 2 + (A^f 



which together with Lemma |6 . 1 2| gives 

E[(e- A Z T - - l) 2 | &v]l Gv ^ Var(e- A Z T - | &«)l Gv + C {(A^) 2 + e 3 ) 



(7.7) 



(7.8) 



by the hypotheses on the initial configuration. Lemma 6.5 and (|7.7|), together with (6.1), (6.4) and 

(7.9) 



(6.15), now give 



B[e- A Z T - | G v ] = 1 + 7T 2 7o(^ + loge + qui o(l)). 



Similarly, (7.7) and (7.8) and Lemma |6.5| give 

E[(e- A Z T - - l) 2 | G v ] = 0(A 2 e 2 ). 



(7.10) 



Furthermore, by Lemma 6.9 and Markov's inequality, together with (6.3) and the hypotheses, we 
have 

P[Y T - > e 3A/2 /a | G*] sj Ce- A/2 . (7.11) 

By the same arguments, and since Z\J =% ee A l/ T=T (\y\ by the definition of a breakout, we have 

(7.13) 



E[e- A Z 0vT - | &v]l Gv < Cel {T=TW) 
E[(e- A Z 0vT -) 2 | & v ]l Gv ^ Ge 2 l (r=T(1)) 



E[y 0vr -|^*]lG« <e^<e- A , 



(7.14) 



be (6.3). Equations ( 7.11[ ) and ( JT.14 ) and Chebychev's inequality applied to (7.10) and (7.13) 
together with (6.1) now give 

P(GuG|G^) < l-e 5/4 . (7.15) 

D 



The lemma now follows from (7.5), (7.6) and (7.15). 



The probability of G n b a b- 

Lemma 7.5. We have P((G n b a b) c I ^~a)1g a ^ Ce 2 for large A and a. 



Proof. Define R\ 
mas 



Rv 



?fug 



[T-,ei+e^ a 2] + R [T,ei+e A a,2] + ^p 1 - 6i+e A a 2 ] ■ % Markov's inequality, Lem- 



5.6, and 5.7 and Corollary |6.7| we have 

P(Ri > | ^a)1g a < C(Z A e A /a + Y A ) sc Ce 2 , 



by the definition of Ga, (6.2) and (6.3). Now, by Corollary 6.13 



E[Z T | #a]1g a < CA&n «S CAe A/3 and E[Y T | ^a]1g a < C<% < C eA/3 > 



?(i) 



6.2 



and the 



where we used the fact that Zlt' < ee by the definition of a breakout. By Proposition 
inequality 1 — e -a: ^ x, the probability conditioned on J^a and on Ga that a check- or bar-particle 
breaks out between T and 0i is now bounded by pse ' 3 =% Ge 2 by (6.14) and (6.2). This yields 
the lemma. □ 
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The particles at the time Oi. The probability of Gi. Recall that from the time T + 



T + t, 



(W,T) 



on, we move the barrier according to the function f A , which is equivalent to having 



(7.16) 



the variable drift —[i t = —fi — fA(t/a 2 )/a 2 . On Gf ng , we have T + ^ T + ( and Q\ = T + e a 2 . 
Furthermore, < A ^ Clog(l/e) on G A and by the hypotheses on the functions (f x ) and (6.2) it 
follows, 

on G A : U/aII < Va and A - ig = °( e ~ A/2 )> 
where || • || is defined in (5.4). 

Lemma 7.6. Suppose (HBq). Then P(Gi) > 1 — Ce 5//4 for large A and a. 

Proof. By Proposition 5.2 Lemmas 5.1 and 7.5 Corollaries 6.7 and |6.13 and ( JT.16 ), 

E nbab [e- A Z ei | ^a]1 Ga = e- A - A (Z A + M«){1 + 0{e- A ' 2 )) = 1 + 0{e- A l 2 ) 
Var nbab (e- A Z ei | ^a)1g a ^ Cee- A - A (Z A e A /a + Y A + &*) sj Cee' 2 ^, 



such that by Chebychev's inequality and (6.2), 



Pnbab(|e- A Z ei - 1| > e 6 ' 2 ) < Ce 



3/2 N 



Furthermore, by Lemma |6.9[ Corollary 6.13 and Markov's inequality, we have 

Pnbab(suppi/ ei c (0, a), y Gl < rj | ^a)1g a < e A /(a??) + Cr}e A sc Ce 2 , 



(7.17) 



(7.18) 



by (6.3) and (6.2). The lemma now follows from (7.17) and (7.18) together with Lemmas 7.4 and 

El □ 



7.6 



obviously still holds if one replaces G\ by G\ n {|e ^^0! — 1| < le 3 ' 2 }. 



Remark 7.7. Lemma 

This will be needed in Sections [8] and 

The Fourier transform of the barrier process. Define Adrift = e ~ A (Z T - + Z^ vT -) — 1 and 

Ajump = e Z> ' ', which are independent random variables. Recall that A > on G A . By (7.12) 
and Lemmas 6.18 and |7.4[ we have 

B[e- A Z 0vT - | G v ] sj CeP(T = T« | G*) ^ Cvle 2 , 



such that with (7.9), we get, 

E [Adrift | G v ] = vr 2 7o (A + loge + c + o(l)). 
Furthermore, (7.10) and (7.13) and the inequality (x + y) 2 < 2(x 2 + y 2 ) yield 

E [(Adnft) 2 I Gw]=0(e 2 A 2 ). 



Equation (7.20) and (6.1) imply 

P(|Adrift| > s 1 / 3 , G« ?/ ) = O^A 2 ) = 0{e 7 h- 



(7.19) 



(7.20) 



(7.21) 



Fix A 6 R. Since log(l + x + y) = log(l + x) + log((l + y)/(l + x)), by (7.16), (7.21) and 



Lemma 7.4 



iXX, 



[00] 



E[r"^"L ] - E[e^ A l (|Adrift|<£l/3) l G J + 0{e^ 
E[e ,Aio g (i + A dr , t ) 1(|Adrift|<ei/3)lGAe 
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iA1 °g(i+ifer) 



(7.22) 



■££>] + 0( e 9/8) > 



for any A £ R. We will first study the term concerning Aj ump . Write Z = Z^ )T > and let p be a 
real- valued constant with |p| < e 1 ' 3 . Then, 



E[e *Alog(l+ -=)■] = E[el Alog(l + ^/) | z > £eA] + O( pa ((G (0,O) )C)) 



r 



g(x) P(e~ A Z e dx \ Z > ee A ) + 0(e 2 ), 



(7.23) 



where 



g(x) = exp ^Alog^l + — 
By definition of ps and rj, we have 



f xP(e- A Z e dx | Z > ee A ) = (p B + 0(r ? ))- 1 E[ e - A Zl (ee A <z ^ ) ] = 7r 2 7o (-log£ + o(l)), 



by (4.6), (6.2), (6.3), (6.14) and (6.13). It follows that 



f 



g(x)P( e - A Zedx\Z>£e A ) 



rCO 

l + avr 2 7o(l + / 9)" 1 (- lo ge + °( 1 ))+ I h(x)P( e - A Z e dx\ Z > ee A ), (7.24) 



where we define h(x) = g(x) — 1 — iX(l + p) xlfx^i) for x > 0. Denote by /i (x) its left-hand 
derivative. Note that |/i(x)| ^ C(l a x 2 ) and |/t~(x)| < C(x _1 a x 2 ) for x ^ 0. By integration by 



parts, (4.5) and (6.13), we have 

rCO 

I h(x)P(e- A Zedx\Z> ee A ) 

= h(e) +P B 1 (l +o(l))( h~(x)P(Z > xe A ) dx + {h{l+) - h{l))P{Z > e J 
h-{x)-dx + (/»(1+) - /i(l)) + o(l)) 

r 00 i 

= vr 2 7o(l + o(l)) h(x)^dx. 
Jo x 

Now, one readily sees that 



(7.25) 



f 00 1 

Jo KX) ^ 



dx = i\{c + o(l) + 



o(p)) + r 

Jo 



iXx 



1 — iAxl^^!) A(dx), 



(7.26) 



where A(dx) is as in the statement of Proposition 7.3 and d is a numerical constant. Equations 



(7.23), (7.24), (7.25) and (7.26) and the Taylor expansion of e x at x = now yield 



E[e4 Alo g (l + 4 f HE) ] 

= exp 



r-00 

.2 I JXx 



7o(iAvr 2 (-loge + c / + o(l) + vr 2 I e lAx - 1 - iXxl {x<1) A(dx)) + 0(e\ \oge\p) . (7.27) 



Coming back to (7.22), we have by the Taylor expansion of (1 + x) at x = 



E r e JA(log(l+A dritt )+0(e|loge|A d] 



ift))- 



(|A drift |< e i/3) 



lG A ] 



E[(l + iAA drift + 0(A 2 rift ))l ( ^^ a3) ] + 0(e 9 / 8 ) 

1 + iA7r 2 7o(A + loge + c + o(l)) + P(T > Vea 3 ) + 0(e 

expiAvr 2 7o(vl + loge + c + o(l) + 0(e 1/8 )), 



by Lemma 7.4 and (7.21) 



by (7.19) and (7.20) 
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where the last equation follows from Lemma 6.3 and the Taylor expansion of e x at x = 0. This 
equation, together with (7.22) and (7.27) and the fact that Aj ump is independent from Adrift, T and 



Gbuik) yields (7.3) in the case n = 1. 



7.3 Proof of Theorems EH and O 

We define the process (Xi') t>0 by X'J = X^ ] - ir 2 At. 

Proposition 7.8. Suppose (HB). Then, as A and a go to infinity, the process (X")t^o converges 
in law (with respect to Skorokhod's topology) to the Levy process {L^jt^o defined in Theorem \l.l\ 

Proof. Conditioning on #o, we can assume without loss of generality that the initial configuration 
v$ is deterministic and that Go is verified. Denote by (J^/')t>o the natural filtration of the process 

J^T - • In order to show convergence of the finite-dimensional 

L*7 J 



X'J and note that M' = &" , _ x 

distributions, it is enough to show (see Proposition 3.1 in [UJ or Lemma 8.1 in [28], p. 225), that 
for every A e R and t, s > 0, 



E 



E[e iXX " + s 



*n 



e i\X'> e SK{\) 



0. 



(7.28) 



as A and a go to infinity. Now, define n := [t^f J and m := [(t + s)7 J, such that (m — n)7o 
s + A~ l o{l), by (|6.1[) and (6.14). Then we have by Proposition 7.3 



E[e ^' +s -^)|^ Ji^ 



In total, we get for A and a large enough, 



-iXAs 



E [ e ^(4 



00] 



v LCC -h , „ 



^eJl Gn 



cxp 



»(k(A) + o(l) + 0(e b )))l 



(7.29) 



L G„ 



E 



E[ e iA(x "+ s - x t') | J?/'] - e SK(A) 1 sj e SK ( A )E[|e s(o(1)+ ° (eb)) - 1|] + P(G 



By Proposition 7. 3| this goes to as A and a go to infinity, which proves (7.28). 



In order to show tightness in Skorokhod's topology, we use Aldous' famous criterion [Tj (see also 
|10j . Theorem 16.10): If for every M > 0, every family of (J^/')-stopping times r = r(A,a) taking 
only finitely many values, all of which in [0, M] and every h = h(A, a) ^ with h(A, a) —* as A 
and a go to infinity, we have 

X" +h — X" — * 0, in probability as A and a go to infinity, (7.30) 

then tightness follows for the processes X' t ' (note that the second point in the criterion, namely 
tightness of X" for every fixed t, follows from the convergence in finite-dimensional distributions 
proved above). Now let r be such a stopping time and let V T be the (finite) set of values it takes. 
We first note that since G n 3 G n+ \ for every n e N, we have for every t 6 V T and every A and a 
large enough, 

P(Gf t7 _ lj )^P(Gf M7 _ lj ) = 0(M £ b ). (7.31) 

by Proposition 



7.3 



Moreover, since &'J a ^r _i f° r every t ^ 0, we have for every A > 0, 

L*7 J 



t) 



teVr 



E[ e ^" +h -*")] = 2 E[e^ x "^- x "h {T 

Mx>{ +h -x») | ^ 

Xk(A) + o(1)+0(£ 6 )) 



V E\E[e tA( ^-^)\3? T ]l(r=t)lG 



0(Me 6 



by (7.31) 



1 



by (7.29) 



0(Me 6 )) + 0(M £ (, ) 1 

which converges to 1 as A and a go to infinity. This implies ( 7.30| ) and therefore proves tightness in 
Skorokhod's topology, since M was arbitrary. Together with the convergence in finite-dimensional 
distributions proved above, the lemma follows. □ 
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A coupling with a Poisson process. Let (V n ) n ^o be a sequence of independent exponentially 
distributed random variable with mean 70. In order to prove convergence of the processes X[ and 
Xf, we are going to couple the BBM with the sequence (V^) in the following way: Suppose we have 



constructed the BBM until time n — i- Now, on the event G n — 1, by Lemma 6.6, the strong Markov 



property of BBM and the transfer theorem (|38j. Theorem 5.10), we can construct the BBM up to 
time 6 n such that P(|(T„ - 6 n _i)/a 3 — V n \ > e 3//2 ) = 0{e 2 ) (recall that T n denotes the time of 
the first breakout after n -i). On the event G^_i, we simply let the BBM evolve independently 
of {Vj)fen- Setting G' n = G n n {ij ^ n : \(Tj - @j-i)/a 3 - Vj\ sj e 3/2 }, there is by Lemma 



6.6 



and 



Proposition |7.3| a b > 0, such that for large A and a, 

P(G' n )>P(G )-nO(e 1+b ) 
Furthermore, we have n = T n + a 5 ' 2 on G' n , whence for large A and a, 

on G' n : |(0 n - e n _!)/a 3 - F n | < 2e 3 / 2 . 
This construction now permits us to do the 



(7.32) 



(7.33) 



Proof of Theorem \7.S\ Let d denote the Skorokhod metric on the space of cadlag functions D([0, 00)) 
(see [2S], Section 3.5). Let <& be the space of strictly increasing, continuous, maps of [0, 00) onto 
itself. Let x,x±,X2, ■ ■ ■ be elements of Z?([0, 00)). Then (|28j, Proposition 3.5.3), d(x n ,x) — > as 
n —* 00 if and only if for every M > there exists a sequence (<p n ) i n ^j such that 



sup |Vn(i) — *| - *• 0, and 
te[0,M] 



SUp \x n ((p n (t)) ~x(t)\ 
te[0,M] 



0. 



(7.34) 



If (x' n ) ne N is another sequence of functions in D([0,ao)), with d(x' n ,x) — > 0, then by the triangle 
inequality and the fact that <3? is stable under the operations of inverse and convolution, we have 



(7.35) 



d(x n , x) — > if and only if there exists a sequence (<p n ) i n 3>j such that the first inequality in (7.34) 
holds and 

SUp \x n ((p n (t)) ~ x' n (t)\ -* 0. 
te[0,M] 

For every A and a, we define the (random) map (fA,a e ^ by 

¥A,aho( n + r)) = ((1 - r)G n + r9 n+1 )a" 3 , with n e N, r e [0, 1]. 
Let M > and define um = [-^7ol- Then we have by the monotonicity of {X\ 



t^o, 



SUp \(pA,a(t) ~ t\ ^ max | a ®n — 70 n 

te[0,M] ne{0,...,n M } 

sup |Xf-X; (t) |< max A|a- 3 6. 

te[0,M] ^ ' V ^ ns{0,...,n M } 



70™ 



(7.36) 

(7.37) 



By Doob's L inequality and the fact that 70 = i£[Vi], we get 



P j max 

■ ne{0,...,n M } 



Yu Vi ~ n7 ° 



i=i 



> e 1/3 ) sj Ce- 2/3 n M V™(Vi) = 0(e 1/3 ). 



Furthermore, on the set G' n , we have for every n ^ um, |@n — Ya=o ^*l ^ 0( n M£ 3 ' 2 ) = 0{e 1 / 2 ). 
In total, we get with (7.36) and (7.37), as A and a go to infinity, 



VM > : sup |yu, a (£) - t\ v |X t " - X^ 

te[0,M] 



VA,a(*)l 



0, in probability, 



(7.38) 
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which is equivalent to 



M=l 



£ 2- M 1 a sup \<p A , a (t) - t\ v \X» - X ' {t) 

L V tG[0,M] ,y ' 



0, in probability. 



(7.39) 



Now, suppose that A and a go to infinity along a sequence (A n ,a n ) ne ^ and denote by X' A a , X" A a 
and (fA„,a n the processes corresponding to these parameters. By Proposition 7.8 and Skorokhod's 



representation theorem (|10j. Theorem 6.7), there exists a probability space, on which the sequence 
(X A ) converges almost surely as n — > oo to the limiting Levy process L = (Lt)t^o stated in the 
theorem. Applying again the representation theorem as well as the transfer theorem, we can transfer 
the processes X' A a and y>A n ,a„ to this probability space in such a way that the convergence in 
(7.39) holds almost surely, which implies that the convergence in (7.38) holds almost surely as well. 



By the remarks at the beginning of the proof, it follows that on this new probability space, 

d(X' An , an ,L) < d(X' Antan ,X' A ^ a J + d(X' AnAn ,L) - 0, 
almost surely, as n — > oo. This proves the theorem. 



□ 



,A,a\ 



Proof of Theorem 7.1. Let (t, ' )™ =1 be as in the statement of the theorem. By the virtue of Theo- 
rem [7T2l it suffices to show that 



'(vi:X{£ 



7 3 t- ' a- 5 



1. 



(7.40) 



Suppose for simplicity that t^ ' a = ti for all i; the proof works exactly the same in the general case. 
Let n := \2(tk + 2)/7o], such that n = 0(e _1 ), by (6.14). By Chebychev's inequality, we then have 

n n 

P(2 Vi < t k + 2) ^ P(^(Vi - 70) ^ -%) = 0{n Var(V 4 )) = 0(e). (7.41) 



i=l 



i=l 



Furthermore, define the intervals 1% = t% + [—2ne 3 < 2 — e /a, 2ne 3 ' 2 ], i = 1, . . . , k and denote by & 
the point process on the real line with points at the positions Vi, V\ + V2, V\ + V2 + V3, . . .. Then 
& is a Poisson process with intensity 7^ = 0(e _1 ) and thus, 



'(fn\J h*0) = 0(e l l 2 ). 



(7.42) 



We now have 



P(Vt:X!.„, 



Lcoj 



Ji 



Ua 3 



>p(t(i,j):t i a 3 e[e j -T j _ 1 ,& j ]j 

n k 

>p(G' n , YjV^tk + 2, ^n [jli 



i=\ 



j=l 



by definition 

by definition of G' n 



> P(G ) - 0(e b ) 



by (7.32), (7.41), (7.42) 



Letting A and a go to infinity and using the hypothesis (HB) yields (7.40) and thus proves the 
theorem. □ 



8 The B b -BBM 

In this section, we define and study the B -BBM, which is obtained from the B-BBM by killing 



some of its particles. It will be used in Section 10 to bound the iV-BBM from below. 
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8.1 Definition of the model 

We will use throughout the notation from Section Yn Furthermore, we fix 5 e (0, 1/2) and K > 1 
such that Ek ^ S/1Q, where Ek was defined in (2.6). Define N = \2ire + a~ 3 e fm ], where fi = 
■\J\ — it 2 /a 2 . We will then use interchangeably the phrases "as A and a go to infinity" and "as 
N go to infinity". Furthermore, in the definition of these phrases (see Section 6.1), Aq and the 



function ao(A) may now also depend on 5. The symbols C$ and Cg a have the same meaning as C 



(see Section 5.1), except that they may depend on 5 or 5 and a as well, a being defined later. 
The B -BBM is then defined as follows: Given a possibly random initial configuration uq of 



particles in (0,a), we let particles evolve according to B-BBM with barrier function given by (7.1), 
where, in addition, we colour the particles white and red as follows: Initially, all particles are 
coloured white. As soon as a white particle has N or more white particles to its right, it is coloured 
recQ Children inherit the colour of their parent. At each time n , n ^ 1, all the red particles 
are killed immediately and the process goes on with the remaining particles. See Figure [3] for a 
graphical description. 




Figure 3: The B^-BBM with parameter N = 6 (no breakout is shown). White and red particles are drawn 
with solid and dotted lines, respectively. The blobs indicate when a white particle is colored red, the circles 
show the six white particles living at that time (thick circles correspond to two particles at the same position) . 



e =5f^ d r if and only if the 



For an interval / cz R + , we define the stopping line -SfS*^ by (u,t) c ~^0 j 
particle u gets coloured red at the time t and has been white up to the time t, with t e I. We 
then set Z^j and KS fj by summing respectively wz and wy over the particles of this stopping line. 
Furthermore, we define jV t and ^ whlte to be the subsets of jVt formed by the red and white 
particles, respectively, and define Zf ed , Y^ ed , Zf^ lte and Y^ whlte accordingly. 



6 This can be ambiguous if there are more than one of the particles at the same position, for example when the 
left-most particle branches. In order to eliminate this ambiguity, induce for every t^Oa total order on the particles 
in ,yV(t) by u < v iff X u (t) < X v (t) or X u (t) = X v (t) and u precedes v in the lexicographical order on U. Whenever 
there are more than A^ white particles, we then colour the particles in this order, which is well defined. 
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Let v\ be the configuration of white particles at the time t and abuse notation by setting 
v n = Vq ■ We set G_ 1 = il and for each n e N, we define the event G n to be the intersection of 
G n _ l with the following events (we omit the braces). 

• supple (0,a), 

• -^ hite c U x {6 n } and 9 n > T+ (for n > 0), 

• |e~^Z2 hite - II < e 3/2 and Y£ hitc sc r>. 

• P ( z d [e„,e, l+ ^a2] + Y d [e n ,e n +i<a2] < a_V2 ^e») > 1 -e 2 . 

The last event is of course uniquely defined up to a set of probability zero. Note that G n e J^Q n for 

each n e N. Furthermore, we define the predicates 
(HB ) (The law of) uq is such that P(Gq) — » 1 as A and a go to infinity. 
(HBq) ^o is deterministic and such that G holds. 

(HBj_) uq is obtained from [e iVj particles distributed independently according to the density pro- 
portional to sin(7rx/a)e _M:r l(o ia )(x). 

We now state the important results on the B-BBM. 

Lemma 8.1. (HB* L ) implies (Hit ) for large A and a. 



Proposition 8.2. P ropo sition 7.3 still holds for the It -BBM, with G n replaced by G n . The same 
is true for Theorems 7.1 and 7.2, with (HB) replaced by (HB'). 

Recall the definition of qu^ and x a from the introduction and of (Ht) from Section ul 

Proposition 8.3. Suppose (HBj_). Let (t- '") satisfy (Ht) and let a 6 (0, e~ ). As A and a go to 
infinity, 

P(Vj : qu^a^.) > x^) - 1. 

8.2 Preparatory lemmas 

In this section, we will establish some fine estimates for the number of particles of the process, which 
will be used later to bound the number of creations of red particles. If a particle u 6 «/^(i) satisfies 
ti(u) ^ t < oi + i(u) for some I ^ 1, we call it an "in between" particle (because it is in some kind 
of intermediary stage) and otherwise a "regular" particle. We then define 

Nt(r) = J] lpr„(*)>r)- 
(u,s)ejVt 

The main lemma in this section is the following: 

Lemma 8.4. Suppose (HBq). Let Ka? < t ^ 46.211 - Then, for ^ r < 9a/10 and every a > there 
exists Cs, a , such that for large A and a, 

P(JV t (r) > e -^ 2+r ^N | T > t) < C 5 , a A 2 s^ + ^e^l"^ 

Furthermore, conditioned on #a? for t < ©i + e a 2 , for large A and a, 

Pnbab(iV i (r) > e~W 2+r MN\ ^a)1 Ga < C^AVe"^^. 

The following lemma about BBM conditioned not to break out will be used many times in the 
proof: 
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Lemma 8.5. Let f be a barrier function, < t < to ^ 46.211 an d suppose that Err(/, t) < C. For 

x e (0, a), define P^ = Pi(- 1 T > to)- Then, for all r ^ (9/10)a and /or /arge ^4 and a, 

B x f [N { t 1+) (r)] sC CAe- A N(l + /xr)e^ r (l + {r)r 2 /a)) (4^(*) + Wy(xj\ (8.1) 



E^(iV t (r))^] sc C(l + r 4 )e^M^e- A A^ ( ^w z (x) + w Y {x) ) . 



i.2 



Remark 6.1). This gives, 



Proof. By Lemma 5.7, the upper bound of Lemma 6.8 is still valid with varying drift (see also 
E/^0%] < CA{±wz(x) + w Y (x)) , and E^ 1 ^] ^ ^[zg^]. (8.3) 



where the second inequality follows from Lemma 5.8 Equation (8.1) then follows from Lemma 5.4 
and (8.3), noting that a — y — f((/a 2 ) ^ r + a/20 for large a, by the hypotheses on r and /. 



For the second equation, suppose for simplicity that / = 0. By Lemma 6.10, we have for all 
x e (0, a), 



E x [(N t (r)) 2 ]^E x [N t (r)]+(3 m 2 f fp{' (0+) (x, z)E< 2 ' fl >[JV t (r)] 2 dzds 

Jo Jo 



E4 ^ J] E^i( Sl )' Sl )[iV t (r)]E( x "2( S2 )' S2 )[iV t (r)]j. (8.4) 



Now, for large A we have by (5.43) and Lemma 5.1 

■*t ra 



pf/ 0+ \x,z)E^[N t (r)] 2 dzds^CE\ V \ p s ^ (X v (s ),z)E^[N t (r)] 2 dzds 

J JO L /„ . „x c - Jso Jo 



(«,ao)e^i 



• 5) 



With Lemmas 5.4 and 6.8 and (8.1), we have as in the proof of Lemma 5.5, for any so ^ t, 

-t ra 



J ^p s - S0 (x,z)E^[N t (r)] 2 dzds ^ CA 2 e- 2A N 2 {l + r i )e- 2 ^[^w z {x)+w Y {x)). (8.6) 



Equations (8.5) and (8.6) together with Lemma 6A 

"t ra 



now give 



f rpf/ 0+ Xx,z)E^[N t (r)] 2 dzds sc CA 2 e- 2A N 2 (l + r i )e-^ r (-, 
Jo Jo Vfl,J 



w z (x)+Awy(x)). (8.7) 



Moreover, we have 



2n 



E x \ J] 2 E x -i^[N t {r)]E x ^M[N t {r)] 

(u,s)eMt (!)i,si),(t)2,s 2 )ey( , '> s ),si<s2 

^E x [ J] ( J] E x ^[N t (r)] 

^C(l + r 4 )e-^M 2 e- 2A iV 2 E x -[ £ ( Z ^) 2 

{u,s)e3e t 

sc C(l + r 4 )e- 2 » r ee- A A 2 N 2 (^wz(x)+w Y (x)), 



where the last inequality follows from (6.17) and Lemma 5.6, Equation (8.2) now follows from (6.2), 
dO, |IO), dO| and (Ob. D 
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Corollary 8.6. With the hypotheses from Lemma 8.5, we have for large A and a, 

% X f[N t {r)} < CAe- A N{l + r 2 )e^ r e- (a ^ )/2 . 



Proof. This follows immediately from Lemmas 5.1[ 5.4 and 8.5 and Corollary 6.7 



□ 



Lemma 8.7. Suppose (HBq). Let a > 0. For large A and a, we have for every < r < 9a/10 and 
t ^ 6i + e A a 2 , 



Pnbab(iV t (r) > e- A / 6 - r / 3 N\^ A )l GA ^ C a A 2 ee- 2A / 3 e-d- a >. 



Proof. Conditioned on J^a> define P = P(- 1 T > Q\ + e a ). Then 



Var nbab (iV t (r)|^' A )l GA ^c2 £ l (ff . <da<Ti) (k u - 1)E 



(Xu(du),d u ) 

xm 



[Nt(r) 2 ], 



i-0 u<9/ 



and Lemma 8.5 now implies 



Var nbab (N t (r) | ^a)1 Ga ^ C(l + r*) e - 2r A 2 ee- A N 2 i 
Furthermore, by Corollary |8.6[ 

E nbab pV t (r) | ^a]1g a < CAe- A N(l + r 2 )e~^S % . 



(8.9) 



5.10) 



J.11 



5.12) 



Equation (8.9) now follows from (8.11) and (8.12) together with the conditional Chebychev inequal- 
ity and the fact that S^ ^ e ' 3 on G<% a Ga- D 



Proof of Lemma 8.J^, Assume Ka < t < te2T} By the hypothesis on Zq, the definition of K, 



Proposition 5.3 and Corollary |6.7[ we have for large A and a, 

E[iV t (0) (r) \T > t] ^ (1 - 75/8) (1 + /xr)e^ r iV, 



5.13) 



Moreover, we have by Lemma 8.5 



r(i+) 



E[iV t li+j (r) | T > t] < C(l + /ir)e-^ r iV(l + rjr)A e - A (^Z + Y ) sc Cy/eA{\ + r 2 )e^ r iV, (8.14) 



by the hypotheses on Zq and Yq. Equations (8.13) and (8.14) now give for large A and a, for all 
re [0,o], 

E[iV t (r) | r > t] < (e^ 2 - <5/8)e~ r/3 iV, (8.15) 



Furthermore, by Lemma |8.5| and the hypotheses on Zq and Yo> we have 

Var(JV t (r) | T > t) < C^ 2 e(l + r 4 )e~ 2r (t/a 3 + r>), 



.16) 



Chebychev's inequality, (8.15) and ( J8.16 ) yield the first equation of the lemma. 

Conditioned on J^a> let T < t < ©i + e a 2 . Define 7V^ (r) to be the number of hat- and 
check-particles to the right of r at time t that have not hit a between T~ and t and likewise 

~r rii 

N t us (r) the number of fug-particles with the same properties. Set M% = e * = (1 + A^) -1 , where 
A t = 9({t - T+)/a 2 )A. Note that we have on G A : |A - e- A Z (a// ^\ < e x /8. 



Define a r = (1 + //r)e ^ r . By Lemma 5.1 and Proposition 5.3, we have for large A and a, 
E[iV t bulk (r) | J^a]1g a ^ a r e o(1) iVM t (l + E^e"^- + Z 0>t -)1 Ga < a r NM t (l - 35/4), (8.17) 
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by the definition of G and G. Furthermore, we have by Proposition 5.3 and Lemma 
E[iVf s (r) | ^ a ]1 Ga ^ a r e°^NM t (9 ((* - T)/a 2 ) + O ((y + A + r) 2 r,/a)) e - A Z^> T h GA 
sc a r NM t (A t + 0{e l ' & + a^r 2 )) 1 Ga , 



5.1 for large a, 
.18) 



Equations (8.17) and (8.18) now give, 



E[AT bulk (r) + Nl us (r) \ ^a]1 Ga ^ a r N (l - 5/2 + 0{ a - l r 2 )) . 



Similarly, one has by Lemma |5. 5 



Var« ulk (r) + A^ ug (r) | ^a)1g a < Ce- A {\ + r A )e 



-2fir 



3.19) 



(8.20) 



Equations (8.19) and (8.20) and the conditional Chebychev inequality now yield for large A and a, 
P(iV bulk (r) + ivf g (r) 5* e ~^/5-r/3 N | ^^x Ga sc C s (l + r 4 ) e - A e~^ r . 



This, together with Lemma |8,7| and the fact that the hat-, fug- and check-particles do not hit a on 
G n hah finishes the proof of the lemma. □ 



We finish this section with a result which fills the gap in Lemma 8.4 under the hypothesis (HBjJ. 

Lemma 8.8. Suppose (H$ L ). Then P(\e~ A Z - 1| sc e 3 / 2 , Y < n) < o(l). Moreover, let a > 0. 
Then for large A and a, we have for every t < Ka 2 and < r < 9a/10, 



P(iV t (r) > e-W 2+r ®N, R Ka 2 = 0) < C a , 5 e- 



-(S-a)r 



/a, and P(-R*r a 2 = 0) ^ 1 - C s e A /a. 



Proof. Recall that initially, there are [e A^J particles distributed independently according to the 
density (p(x) proportional to sm(irx/a)e~ flx l(o :a -)(x). An elementary calculation yields that 



E[Z ] = e A (l + o(l)), Var(Zo) sj Ce A /a 6 , E[Y ] sc Ce A /a. 



3.21) 



This immediately yields the first statement, by Chebychev's and Markov's inequalities. Moreover, 



.21) with Lemmas 5.6 and 5.7 yields 



/K \ Ke A 

P(R Ka2 > 1) < E[R Ka2 ] sc C(-E[Z ] + E[Y ] < C 1 , 

V a / a 



(8.22) 



which gives the last statement. Note that on the event {Rko 2 = 0}> the B -BBM equals BBM with 
absorption at and a. Since the density <f>(x) is stationary w.r.t. this process, a quick calculation 
shows that 



ELV t (0) (r)] = E[N^°\r)] = [e~ 5 N\ f° <f>{x) dx 



< e ^(1 + /w)e~ 



A"' 



Furthermore, by the independence of the initial particles and the law of total variance, 
Var(AT t (0) (r)) = [e~ s N\( Var(E x LY t (0) (r)]) + E[Var x (AT t (0) (r))] 



(8.23) 



124) 



where AT is a random variable distributed according to the density (j) and the outer variance and 
expectation are with respect to X. By Lemma 5.5 we have for every x e [0, a], 

Var x (AT 4 (0) (r)) < ^[^(r) 2 ] sj Ce- 2A N 2 (1 + r 4 )e- 2 ^ r ((t/a 3 )w z (x) + w Y (x)), 



and a simple calculation then yields for t =% Ka 2 , 

E[Var x (A^ (0) (r))] = f° V&r x (N { t °\r))cb(x) dx sc C{K/a)e- A N{l + r 4 )e~ 2 '"'- 
Jo 
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3.25) 



Moreover, by Lemma 5.4 and the hypothesis on r, we have for every x e [0, a], 



E*[iV t (0) (r)] < C(l + r 2 )e-^(e^l (x ^ ) + e- A N(wz(x)+w Y (x))l ix> i^a ) 



which yields 

E 



19a 
20 



E A [A; m !/)))- > | ^C(l + r 4 )e- 2 ^( I "" ae^ dx + e' 2A N 2 f 
^Ca- 3 e- A N(l + r 4 )e- 2tMr . 



(a — x + 



1 N3 e Kx-2a) dx 



J ; i - j 

20 



Equations (8.24), (8.25) and (8.26) now yield for large a, 

Var(iV f (0) (r)) sc C{K/a) e - A N(l + r 4 )e _2/xr . 



The lemma now follows from (8.23) and (8.27), together with Chebychev's inequality. 



(8.26) 



127) 



□ 



8.3 The probability of G\ 

Most of the work in this section will be devoted to bounding the number of creations of red particles. 
For this, we will discretize time: We set t n = Ka 2 + n5/A and I n = \t n ,t n+ \) for all n 6 N. 
Furthermore, define the &<% -measurable random variables n\ := max{n : t n < T} and ri2 ■ = 
max{n : t n < ©i + e A /a 2 }. 

For an interval / c: R + and r e [0, a], we now define 

Uj{r) = #{(«, t)eJ?g d j:X u (t)>r}, 

the number of red particles created to the right of r during the time interval / and set U n {r) = Uj n (r). 
Furthermore, we denote by iV t te (r) the number of white particles with positions > r at time t, 
including the "in between particles" . 

Lemma 8.9. Suppose (HB$). Let I = [£j,i r ], with d := t r — t\ =% (5/4. Then for every a > 0, 
for sufficiently small 5, there exists C$^ a , such that for large A and a, for every r e [0,9a/10], if 
Ka 2 < tj < 



E[Uj(r) |T > t,] < C 5 , a A 2 e( ^ + r/je-tf-^iV, 
P(3t e / : iV t white (r) > JV | T > *,) s; C 8 ^A 2 e(\ + r/)e- ( t" a)r , 



Ep7/(r)l (J j ttj _ Ciir]=0) | T > t,] < Q,a^ 2 e(^ + r^e-d^iV. 



Likewise, conditioned on j^\, we /iai>e /or T < i; < Oi + e a 2 , 

E nbab [[//(r) |^ a ]1 Ga < C (5)a A 2 £ 2 e-(f- Q >iV, 
P nbab (3t e / : A7 hitc (r) 5* N | ^a)1 Ga ^ C^Ve" (§-">. 

Finally, suppose (HEr^). Then 



E[C/ [0 ^ a2] (r)l (ilKa2=0) ] ^ C s , a ae-^-^N. 



(8.28) 
(8.29) 
(8.30) 



.31) 
.32) 



.33) 
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Proof. Morally, (8.28) follows from the fact that particles get created with rate N/2 and the prob- 



ability that a particle at the point r turns red is given by the probability that there are more 



than N particles to the right of r, an upper bound of which is given by Lemma 8.4 If there were 
no in-between particles, we could make this argument precise and give a short proof for the ex- 
pected amount of red particles created in an infinitesimal time interval. However, the presence of 
in-between particles, which do not follow Markovian dynamics due to the conditioning, makes the 
proof more complicated and forced us to discretise time. 

Define the stopping time r = inf{£ e I : N™ e (r) > N}, with inf = +00. Define further the 
law P = P(- \T > ti) and the event E = {N tl (r') < e" 5 / 2 "^/V}, where r' = r - «fr v 0. Then 
e -4r /3 «^ C a e~ r ( 4 ' 3 ~ a > for every a > 0, such that by Lemma 



8.4 



P(r < 00) ^ P(r < 00, E) + P(E C ) < P(r < 00 | E) + C 5 , a A 2 e(^ + rj e"^"^. (8.34) 



we have 



*-(*L 4- X-(f-«)' 



For t e /, we can bound A^ t whlte (r) by the sum of 

1. Nf : the number of descendants at time t of {the particles in JV^ which are to the right of 
r'}, 

2. N' 2 ^: the number of descendants of {the particles in JVti which are to the left of r'} and which 
reach the point r before the time t r and 

3. iV t : the number of "in between" particles at time t. 

Conditioned on ^ ; and the event E, Nj; is stochastically bounded by the number of individuals 
at time t — ti in a Galton- Watson process with reproduction law q and branching rate /3o, starting 
from [e~ <5//2 ~ r / 3 iV| individuals. Let (GWf)^o be such a process, then e"* mt GWt = e~ t//2 GWi is a 
martingale and by Doob's L 2 -inequality, we get 

E[ sup {e- t/2 GW t - GW ) 2 ] sj 4Var( e - d / 2 GW d ) < CdN e - r ' /3 . (8.35) 

te[0,d] 



By (8.35), Chebychev's inequality and the hypothesis on d, we get, 



P(supiV t [1] > Ne~ s/4 ) < P( sup GW t > Ne' 5 ^) 
tel ' te[0,d] 



sj P( sup e- t/2 GWi - GW > Ne- 3S / 8 (1 - e^ 8 )) ^ C 5 d e - r ' /3 /N. (8.36) 
te[0,d] 

Now, if r < 1, we have N* 2 ) = 0, so suppose that r ^ 1. Let the random variable T T denote the 
number of particles in BBM which reach r before the time d. Since there are at most N white 



particles to the left of r' at time ti, we have by Lemma 3.1 



ELyM I & tl \ sc N sup B x [T r ] sj CNW( sup B t > y/f) < CNd- 1/2 e- r/( - 2d) , (8.37) 

xe[0,r'] te[0,d] 

where (Bt) follows standard Brownian motion under W. F or each (u,t) e JVti-, let Nu be the 



contribution of the descendants of u to N<- 2 K By Lemma 3.5 

VarfiVP] I & tl ] < 2 E [( A7 'i 2] ) 2 ] < N SU P EX [ r r] 

(u,t)eJH, xe[0,r'] 

d a ( 8 - 38 ) 

<N( sup E^rj + cf I Pt (x,z)E( 2 ' t )[r r ] 2 dzdt). 
^xe[0,r'] Jo Jo ' 

Trivially, E x [r r ] ^ e d < 2 and $£ Pt(x,z) dz ^ e d / 2 for any x e R, t ^ d Together with ( |8.37| ) and 



(8.38), this yields, 



Var[JVM I J*tJ sJ CiV sup E z [r r ] sc CNd- 1/2 e- r/(2d) . (8.39) 

xe[0,r'] 
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Chebychev's inequality and the inequality t 1 ' 2 e 1 /' 4i ) ^ Ct now yields 



UO) 



As for iVM, we have E[12?] = Var(i? t2 ) + B[R tl ] 2 < Ce 2A by Lemma 



6.12 



6.8 



as well as Lemmas 5.6 
and Corollary [6/7] such that by the definition of a breakout event, E[sup te j(iV t ) 2 ] ^ Ce 2 ( 2 



Markov's inequality then gives 



P(supiV t 
tel 



[3] 



> 



( e -5/4 _ e - 5 / 8 )iV) <c C 5 e 2A C 2 /A r2 



.41) 



Equations (8.34), (8.36), (8.40) and (8.41) now yield (8.29). As for (8.28) and (8.30), we first 



note that the expected number of individuals created until time t in a Galton- Watson process with 
reproduction law (<?(fc))fe>o is bounded by the number of individuals at time t of a Galton-Watson 
process with reproduction law (1 — q(0)) q(k)k^i, whose expectation is bounded by Ct. This gives, 

Ep7 (T)tr] | r < oo] < C(dN + (e A ) . (8.42) 

because there are at most N white particles at the time ti and the total number of in-between 



particles is bounded by £E[-RtJ < CQe by Lemmas 



5.6 



and 



6.8 



and Corollary 



6.7 



As for the 



overshoot, i.e. the particles that have been coloured red at the time r, we first note that 
E[[/ {t} 1 (t<oo) ] ^ B[U {t} 1 {t<OD:U{t] ^ n) ] + NP{t < oo). 



Now, 



U{ T }l( T <'X,u, T} >N) < max{/c u : ue U,d u e I, k u > N, X u {d u -) > r} 



(8.43) 



.44) 



By the definition of a breakout event, the set on the right-hand side of (8.44) contains no "in 



between" particles for N > (. Furthermore, the expected number of branching events of "regular" 
particles to the right of r is bounded by 



i 



2\„-r 



C\ B[N t (r)]dt^CdN(l + r 2 )e 



where the inequality follows from (8.13). This yields, 

E[tf M l( T <ao,tf M >tf)] < CdN{\ + r 2 )e- r E[Ll {L ^ N) ] sj Cd(l 



r )e 



.45) 
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where the last inequality follows from the Cauchy-Schwarz and Chebychev inqualities. By ( |8.29| ), 
(18.421), (18.431) and KM we finally get for large A and a, 



B[U I (r)]=E[(U {T} (r) + U { 



.r. tr ](r-))l(r<x)] < C s , a NA 2 e^ + v y-(i- a > 



S.47) 



For the proof of (8.30), we note that on the set Rn,— ct r ] = there are no in-between particles 



during /. Taking away the corresponding terms in the above proof yields (8.30). 



The proof of the other three equations is very similar. The proof of (8.31) and (8.32) uses the 



second half of Lemma 8.4 instead of the first and (8.12) in addition to (8.13) for the proof of (8.45). 



The proof of the last equation draws on Lemma 8.8 instead of Lemma 8.4 and requires covering the 
interval [0, Ka 2 ] by pieces of length d. □ 

Corollary 8.10. Suppose (HBq). For an interval I a R +; let G a /2j be the event that no particle 
is coloured red during the time interval I and to the right of a/2. Then, 



?(G v \G a/2t[Ka 2 tTi ) < C s o(l) 



and 



Pnbab ( (G a / 2 ,[T,ei+e A a 2 ]) C I <^A ) 1<3 A ^ C S o{l). 
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Proof. Follows immediately from (8.29) and (8.32) by summing over the intervals I n (note that 
T ^ ^/ea i on G f ). 



□ 



Lemma 8.11. Suppose (HB$). Set T" = T a yfea and define the intervals I = [Ka ,T'] and 
J = [T, ©i + e a 2 ] For large A and a, 

E[(Z^j + Y^)l Ga/2I ] < C s A 2 e\ A , and E nbab [(^ cd J + 1^)1^ I ^a]1g a < C s e 2A /a. 

Now suppose (HE>J. Then E[(Z^ Ka2] + Y^^l^^} ^ C 5 e A /a 2 . 



Proof. By Lemma 8.9 we have for every r 6 [0, a/2] and n ^ no := [(y/ea — Ka )/(ii — to)\, 



«() 



»■(> 



E[C//(r)] < 2 E [^( r ) I T > ^ P ( T > *») < C^ 2 ^-^^ £ R + V)P(T > t n ) 



n=0 



sj C 5 A 2 ee^ r e A E 



n=0 

T a v^a 3 



f] 



< C^Ve - ?*", (8.48) 



by Lemma 6.5 and (6.3). Now, by integration by parts, we have for every (possibly random) interval 
J, 



Z T S d i = f «4(r)^/(r) dr and *£* = ^(0) + f u^ 
Jo Jo 



w Y (r)Uj(r) dr. 



.49) 



The first equation in the statement of the lemma now follows from (8.48) and (8.49), since 



w' z (x) =% (7(1 + x)e x a and w' Y {x) =% Ce x a for x e [0, a]. The remaining equations follow similarly. 

□ 



The following lemma will only be needed in the proof of Proposition 8.3 



Lemma 8.12. Suppose (Hb\). Let r e [0, 9a/10] and (K + l)a? ^ t ^ fen - ITien, /or a// a > 0, 
for large A and a, 

P[N r t ed ' {0) (r) >aN,T>t]^ (^(cT^l + ^r) e -^ r A 2 e^ + rj\ + e< 

Proof. For an interval 7, denote by N t e (r, /) the number of red particles to the right of r at 
time t which have turned red during /. Define the event 

zred i T^red 



G = G a /2,[o,t] n {R[t- a *-c,t] = 0} n {^0 i[OiA - a2] + Y0 t[OtKa 2 ] < a 



-l/2i 



and note that P(G C , T > t) ^ Ce 2 for large A and a, by Corollary 8.10, Lemma 6.9 and the 
hypothesis. 

Write E = E[- 1 T > t]. KE>' s )[iV t (0) (r)] < /(», «) for some function / with /(0, s) = f(a, s) = 0, 
then by integration by parts, 



E[iY t rcd '(°)(r,/)l G ] ^ ^^ 2 E[U ds (x)l {R[t _ it ^ 0) ] 



dx 



f(x, s) dxds. 



Define /i = [0,t-a 2 ] and I 2 = [t-a 2 ,t]. By (5.20), (5.21) and Corollary 6.7 we have for x e [0,o/2], 



£(*.«) [JV t (0) (r)] ^ 



C(l + /ir)e _Mr A r e" A u' Z (x), ifaeii, 
Ce^ S^ e~ fiz p^ s (x,z), if seJ 2 , 
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Similarly to the proof of Lemma 8.11 we now have by the inequality w' z {x) < C(l + x)e x a , for 
large A and a, 

E[JV t red ' (0) (r, h)l G ] sj Ce- A iV(l + /xr)e^ r E[C/ /l (x)l G ](l + x)e x - a dx 

Jo 

sj C s A 2 e((t/a 3 ) + r?)(t/a 3 )(l + fi^e^N, 



by Lemma 8.9 and the definition of G. As for the particles created during I 2 , we have by Lemma 8.9 
the definition of G and the hypothesis on t, for all s e I2, 

B[U ds (x)l G ] sj CsA 2 e(X,+ri]Ne-l x ds. 



Now, with (2.5), one easily sees that for large a, {d/dx){e^ x p a s {x, z)) > for all x < a/2 and s ^ a 2 , 
whence, by (2.8) and Fubini's theorem, we get 



E 



■a 2 ra/2 



[N r t cdM (r,I 2 )l G ] ^C\ \ B[U ds (x)l G ] 

JO Jo Jr 



d 



fiz "• / „fJ,x a 



JO 
■a/2 



d.r 



(e^rffoz)) dzdxds 



< 



r a/z -. r 00 d 

C E[U ds (x)l G ] e-^ — ie^a^ix a z)(a - x v z)) dzdx 

Jo J r dx y 



^C s A 2 e{-,+r])N 



J -00 ra/z 
I e^ (:c/4+z) (l + x)dxdz 
r Jo 



sc C s A 2 e{{t/a 3 ) + r^e'^N. 



In total, we have 



E[N T t edM (r)l G ] < C 5 A 2 e(-^ + 77) (1 + /zr)^ r iV. 



An application of Markov's inequality finishes the proof. 



□ 



Lemma 8.13. Suppose (HB* ). Let Gf red be the event that the fugitive does not get coloured red. 
Then, for large A and a, 

P(G& n G a / 2) [Ka 2 ,T]\Gired) < C&E . 

Proof. Set G a / 2 = G a / 2y [Ka 2 ,T]- Let Po be the law of BBM with absorption at zero from Section [6| 
(i.e. we do not move the barrier when a breakout occurs). Let T red be the first breakout of a red 
particle. Then 



P(G* n G a/2 \G fred ) = Po(T = T red , % n G a/2 ) < P (T red sc V^ 3 , G a/2 ). 
Define T" = T a y^a 3 . Then, since -5^f^/ is a stopping line, we have by Proposition 



Corollary 8.10 



6.2 



Po(T rcd < Vea 3 I ^ rcd/ )l Ga/2 < Ce-^e-^Z^, + Y£%,). 



.50) 



and 



.51) 



Recall that P(Z v ^ r 



Y; 



red 



0,T' 



_ 21 < a" 



J 0,[Q,Ka 2 ] ^ * 0,\O,K a 2 ] 

conditional expectation and (8.51), this gives 



1 ' 2 ) > 1 — e 2 by hypothesis. By the tower property of 



Po(T rcd ^ V^ 3 , Ga/2) < E[P (T red sc y/ea s | ^ rod/ )l Ga/2 ] + e 2 + o(l) 



< 



CE[Ce 1/2 e A (^0 d [Ka 2 jT , ] + Y^ d [Ka 2 T/] )l Ga/2 ] + Ce 2 



The lemma now follows from (8.50) together with Lemma 8.11, (6.1) and the hypothesis 



□ 
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Lemma 8.14. Suppose (HBq). There exists a numerical constant b > 0, such that for large A and 
a, 

P(G\) ^l-e 1+b . 

Proof. Recall the event G\ from Section [7] and change its definition slightly by requiring that 
\e~ Zq 1 — 1| ^ e 3 > 2 /2, call the new event G^. Define the random variable 



zrcd 



,[ei,6i+Xa 2 ] 



^red 

,[Gi,Gi+ftTa 2 ] 



- 1/2 |^ 01 ). 



By Proposition 7.3 and Remark 7.7, it suffices to show that 

P (G[\ te d sj e 3 / 2 /2, X > 1 - e 2 }) < e 1+b , 



for some numerical constant b. 



Let T' , I and J as in Lemma 8.11 We then have by Markov's inequality, Lemma 8.11 and 

(8.52) 



Corollary 8.10 



E nbab [l-X|^ A ]l GA ^C 5 o(l). 

Markov's inequality applied to ( 8.52[ ) then yields 

P(X > 1 - e 2 , G A n G nbab ) sj C s o{l). 



(8.53) 



Now define G red = {T < ^/ea 3 , Z r ^ Ka2] + Y^ Ka2] «c a l l 2 } n G / 2) j n G frcd . Conditioned on 
T and J^red , on the set G re d, the particles from the stopping line ^q% then all spawn BBM 
conditioned not to break out before T (because neither the fugitive nor any in-between particles 



are on the stopping line). By Lemma 8.11, Lemma 6.9 and the tower property of conditional 
expectation, we then have 

E[^ d l Grcd ] ^ c(Aa- 1 ' 2 + E[(Zg d 7 + AY£ d j)l Ga/2I ]) < C s {Aefe A 

E[^ ed l Grcd ] < G(Aa- 1 ' 2 + E[(Zg d + Y^j)l Ga/2I ]) < C s A 2 e 3 e A . 

By Lemma |6.9[ we now have 

E nbab [^ d | ^ A v ^ t ]1 Ga ^ Cs(z? d + Z£j + A(Y± cd + Y^j)), 



and Lemma 8.11 and the tower property of conditional expectation give 

E[Zgfl GAnGnbabnGrcd ] ^ C 5 (Aefe A . 
Furthermore, by the hypothesis, Corollary |8. 10| and Lemmas |7.4[ |7.5] and 8.13 we have 

P(G A n G nbab n G rcd ) > 1 - e 1+b . 



154) 



5.55) 



The lemma now follows from (8.53), (8.55) and Markov's inequality applied to (8.54), together with 
(|6lb. D 



8.4 Proofs of the main results 



Proof of Lemma 8.1, By Lemma 8.8, it remains to estimate the probability of the last event in 



the definition of G . Define the random variable X 



"p ( vred 
tr ^ Zj 0,[O,Ka^ 



^1,^1 ^«- 1/2 W- 



By Lemma |8.8| and Markov's inequality applied to the last equation of Lemma 8.11, we have 
E[(l — X)lf R 2= o)] ^ Cse A /a. By Markov's inequality, we then have 

P(X < 1 - e 2 , R Ko? = 0) < e- 2 E[(l - A-)l (fl o)] < C s o(l). 



Together with (8.22), this finishes the proof. 



□ 
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Proof of Proposition 8.2 . Let n > 0. Conditioned on ^@ n , the barrier process until n+ i is by defi- 
nition the same in B-BBM and B-BBM. Furthermore, G n a G n and G n e ^e„. The first statement 



then follows by induction from Lemma 8.14, as in the beginning of the proof of Proposition |7. 3 



As for the second statement, inspection of the proofs of Theorems |7.1| and 7.2 shows that they 



only rely on Proposition 7.3 and on the existence of the coupling with a Poisson process constructed 



in Section 7.3 But this construction only relied on the law of T n+ \ conditioned on J^e n an d G n 
and thus readily transfers to the B-BBM. □ 



Proof of Proposition 8.3, Define p 
Section 



7.3 



Jb r^olS 



, such that (1 + up)e~ w ;> ae 2d (l + o(l)). Recall from 
the coupling of the process (6 n /a 3 ) n ^o with a Poisson process (V n ) n ^o of intensity 
7Q -1 = pBe A 7r. As indicated above, this coupling can be constructed for the B-BBM as well. In 
particular, setting & n = G* n {Vj sj n : \(Tj - Qj-i)/a 3 - Vj\ ^ e 3/2 }, we have as in ( |7.32[ ), for 
every n > 0, 

P(Gi) > P(Gl) - nO(e 2 ) > P(G ) - nO(e 1+b ), (8.56) 

by Proposition 
suppose that t 



.2 



A,a\ 



(where b > is a numerical constant). Now let (t i ' ) satisfy (Ht). For simplicity, 
= ti for all i, the proof of the general case is exactly the same. We want to show 

hite/ 



that P(V1 ^ i ^ k : N^ e (p) 5s aN) -> 1 as A and a go to infinity. By ( |7.42[ ), the probability that 
there exists (i, j), such that ti and tj+i are in the same interval [@j, Qj+i] is bounded by 0(y/e) for 
large A and a. We can therefore suppose that k = 1, the general case is a straightforward extension. 
It thus remains to show that for every t > 0, P(./V^3 lte (p) > aiV) -> 1 as i and a go to 
infinity (the case t = is a straightforward calculation). Let £ > and no := \jq (t + 1)]. Then 
E [^n ] > * + 1 and Var(F no ) ^ 2-/ (t + 1), such that P(K„ < t + 1/2) < C70 = O(e) (for this proof, 
we allow the constants C and the expression O(-) depend on t). Defining the event E n {s) for n ^ 
by 

E n (s) = {Q n >8}n{ije{0,...,n}:se [Tj-i, @j + e A a 2 ]} n G b „, 

(define T_i = 0) we then have by ( |7i42| ) and J8^56| |, P(P„ ) ^ x ~ Ofc 6 )- 

We now prove by induction that there exists a numerical constant b > 0, such that for every 
1 < n < no, 

S„ = sup sup P u (NJ hitc {p) < aN, E n {s)) sj C s , a rie 1+b , 
ueG\ se[0,ta 3 ] 

as A and a go to infinity (we accept the abuse of notation v e G ). By definition, Sq = 0. Now, 
suppose the statement is true for some n ^ 0. We then have for every v e G and s e [0, £a 3 ], 



p , (iV wmtc (p) < aiVj £ n+1 ( s )) sj P^(iV s willtc (p) < «JV, En+iCa), a > 6i) 

+ P-(AT s white (p) < aAT, E n+1 (s), s < 9i) =: Pi + P 2 . 

We first have, by the definition of E n+ i(s) and the fact that the process starts afresh at the stopping 
time 0i, 

A < E*TP(iV7 hite (p) < aN, E n+1 (s) | ^ejl^e^}] ^ S n . 

Second, we have again by the definition of E n+ i(s), 

P 2 ^ P^iV7 hitc (p) <aN,T> s)l { ^ eAa2) ^ P»(N: hitc +°\p) <aN,T> s)l {s ^A a ,y (8.57) 



Write P v = P u (- 1 T > s). For e a < s ^ ^ |6.2ik . we now have by Proposition 5.3 and Corollary 



6.7 



for the first inequality and Lemma |5.5| and Corollary 6.7 for the second, 

B u [N^(p)] > (1 - 0(e 3/2 ))e 25 ae- 5 N > (1 - 0(e 3/2 ))(l + <J)a, 
Var"^^)) ^ C Q a 2 e- 2j4 iV 2 , 
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such that 

P y (JVJ hite >(°)(p) < aN, T > s) < P u (N(°Xp) < (1 + 5/2)aN) + P U (N T S C d +°\p) ^ (5/2)aN, T > s) 

^Cs, a {e- 2A + Ah(s/a s + V )). 



by Chebychev's inequality applied to the first term and Lemma 8.12 to the second. Altogether, this 
gives, 

P2 ^ C S JP(T > Via 3 ) + e- 2A + A 2 e 3 / 2 ) < C 8 , a A 2 e 3 / 2 , 



for large A and a. Finally, we get 5 n+ i < S n + Cs, a A 2 e 3 ' 2 for all ^ n < no, which yields 
S n ^ Cs^nz 5 for large A and a, by (6.1). This gives 



P(^3 ite (p) < aN) < S no + P(K < t + 1/2) + (1 - P(E no (s))) - 0, 
as A and a go to infinity. The statement follows. □ 

9 The Btt-BBM 



In this section, we define and study the B"-BBM, a model which will be used in Section 10 to bound 
the iV-BBM from above. 

9.1 Definition of the model 

As in the previous section, we will use throughout the notation from Section [7J fix 5 6 (0, 1/2) and 
define K to be the smallest number, such that K < 1 and Ek ^ 5/10. We now define however 
iV = \2ire a^ 3 e fia \. Again, we will use interchangeably the phrases "as A and a go to infinity" 
and "as ./V go to infinity" and Aq and the function ao(A) in the definition of these phrases may now 
also depend on 5. The symbols C$ and Cs, a have the same meaning as in the last section. 

The B"-BBM is then defined as follows: Given a possibly random initial configuration uq of 



particles in (0,a), we let particles evolve according to B-BBM with barrier function given by (7.1) 
and with the following changes: Define t n = n(K + 3)a 2 and I n = [t n ,t n +i) (note that these 
definitions differ from those of Section [8]). Colour all initial particles white. When a white particle 
hits during the time interval I n and has at least iV particles to its right, it is killed immediately. 
If less than iV particles are to its right, it is coloured blue and survives until the time t n+ 2 a Oi, 
where all of its descendants to the left of are killed and the remaining survive and are coloured 
white again. At the time 0i, the process starts afresh. See Figure [4] for a graphical description. 

For bookkeeping, we add a shade of grey to the white particles which have hit at least once 
(and call them hence the grey particles). We then add the superscripts "nw", "gr", "blue" or "tot" 
to the quantities referring respectively to the non- white, grey, blue or all the particles. Quantities 
without this superscript refer to the white particles. 

In particular, we define Bj and -Bj ot to be the number of white, respectively, white and grey 
particles touching the left barrier during the time interval I with less than iV particles to their right 
(i.e. those which are coloured blue). We set B n = Bj n and B^ = B^. 

Let v\ be the configuration of (all) the particles at the time t and abuse notation by setting 
vli = v q . We set G_i = £1 and for each n e N, we define the event G n to be the intersection of 
G n _ l with the following events (we omit the braces). 

• suppi^ c (0,a), 

• ^ef c U x {9 n } and G n > T+ (for n > 0), 

• |e-^Z* ot - II < e 3 / 2 and Y* ot ^ ri, 
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< N particles! 



< N particles! 



Figure 4: The B^-BBM with parameter N = 3 (no breakout is shown). White, blue and grey particles are 
drawn with solid, dotted and dashed lines, respectively. A cross indicates that a blue particle is killed. 



N(O n ) ^ N and for all j 5= with 6 n _i + tj < 9 n : iY(9 n _i + tj) > N, 



* I: ^ B [B n ,0 n +u\ 



j < e- A e fia /a 



■F< 



e„ 



^1-e 2 



The last event is of course uniquely denned up to a set of probability zero. Note that G n e && n for 
each n e N. Furthermore, we define the predicates 
(HB") (The law of) uq is such that P(Gq) — > 1 as A and a go to infinity. 
(HBq) uq is deterministic and such that G holds. 

(HB^) uq is obtained from \e N] particles distributed independently according to the density pro- 
portional to sin(-7rx/a)e~ Ma; l(o ! a)(a ; )- 
We now state the important results on the B"-BBM. 

Lemma 9.1. (HBj_) implies (HIt) for large A and a. 



Proposition 9.2. Proposition 1.3 still holds for the B*-BBM, with G n replaced by G n . The 



same 



is true for Theorems 1.1 and 1.2, with (HB) replaced by (HB*). 

Recall the definition of qu^ and x a from the introduction and of (Ht) from Section pi 
Proposition 9.3. Suppose (HEr ± ). Let (t ■ ,a ) satisfy (Ht). Let a e (0, 1). As A and a go to infinity, 

P(Vj:qu^(zA , .,) < .»•,.. ,..0 - L. 






Lemma 9.4. Define a variant called (J-BBM of the &-BBM by killing blue particles only if there 
are at least N particles to their right. Then Propositions 9.2 and 9.3 hold for the G^-BBM as well. 



9.2 Preparatory lemmas 

We first derive upper bounds on the probability that the number of white particles is less than N 
at a given time t. We do not impose any particular condition on the initial configuration. 
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Lemma 9.5. Let P = P(- \T > s) with s =% / jg^il and let Ka 2 < t ^ s. We have for small 5 and 
for A and a large enough, 



r(0) 



P(Nr <N\^ )l (ZoHl _ s/2)eA) ^ C 5 e 



-At„-\ , „-A 



e- A Ya). 



Proof. By (5.20), Corollary |6.7| and the definition of K, we have for A and a large enough, 
E[iV t (0) | & ] > (1 - Cp B )E[iV t (0) I J?o] > (1 + V^e^Zo. 



By the conditional Chebychev inequality, we then have 



r(0) 



Var(iV i (0) | ^ ) 



p(ivr < n i ^0)1(^(1-^)^ < c ,/ m 



The lemma now follows from Lemma 5.5 and Corollary 6.7 



□ 



Corollary 9.6. Under the conditions of Lemma 9.5, suppose furthermore thatt ^ (K + l)a 2 . Then, 
for small 5 and for A and a large enough, 



r(0) 



,(o) 



-A, 



,-A. 



P(iVr < N, Z^' Ka2 > (1 - <5/2)e A | ^ ) < Qe" A (l + e - A Z )/a. 



Proof. First condition on ^t-Ka 2 an d apply Lemma 9.5 Then condition on J^o an d a Pply (5.10) 
and Corollary |6.7| □ 



For s ^ 0, define the event 



Gz,. 



sup 



A ~(0) 



e- A Z 



<<y/4k 



and set G^.n = G\z t„ . In using the last two results, the following lemma will be crucial: 

Lemma 9.7. Suppose |e" A Z - 1| < V 8 - Let s < fem a™d se£ P = P(- |T > s). T/ien P(G^ ) < 
C^e /or Zarge A and a. 



Proof. Let h(x,t) = P<>>*)(T > s) and note that h(x,t) >1- Cp B by (|6.45[). Define 



z 4 fe = 2 (&(*«(*)> *)rW*«(*))> 



U gJ* 0) 



which is a supermartingale under P by Lemma 5.9, with E[Z^] = (1 + 0(pb))Zq by Corollary 
For large A and a, we then have by Doob's L 2 inequality, 

P(G|J sj P( sup \Z y t l - V[ZH ^ e A Vl6) «S C s e- 2A Var(Z%). 

0sCt:£s 



6.7 



The lemma now follows from Proposition 5.2 and Corollary 6.7 



□ 



The following lemma is the analogue of Lemma 9.5 for the system after the breakout. For an 
individual u, we define N^ u and _/V t ug '^" to be the number of hat-, respectively, fug-particles which 
are not descendants of u. 

Lemma 9.8. Let t e [T, 0i] and let (u,t u ) e ,yV T - u ^C^ > T ) . Then, for large A and a, 

Pnbab(A7 u + N^ u < N | J^a)1 Ga ^ Csve- A /e. 
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Proof. As in the proof of Lemma 



8.4 



set Mi 



-x 



[i] 



(1 + A*)" 1 , where A t = 6{{t - T+)/a 2 )A. 



u 



Recall that on G A - | A - e- A Z^^\ ^ e 1 / 8 . For a particle (u, t u ) e ,jV t - u y(®> T ), we define N^ u '° 
and N t s '^ u ' to be the number of hat-, respectively, fug-particles which are not descendants of 
and which have not hit a after the time T~, respectively, T. 

Now, by Lemma |5.1| and Proposition |5.3[ we have for large A and a, 

E[N^ U '° | J^a]1 Ga > Ne- A M t (Z T - - |u>z|| x )(l " */8)1g a ^ tfM t (l + 5/2) 1 Ga , and 

E[Ar fug,o,^ | ^ a j 1ga = AT e -^M t Z^' T )(^((t - T)/a 2 ) + 0(y/a))l GA 
>NM t (A t (l + S)-S/8)l GA . 

In total, this gives for large A and a, 

E[JV7 U '° + Af "° I ^ a ]i Ga ^ N (l + S/4)1 Ga . 

Moreover, by Lemma |5.5[ we have for large A and a, 



(9.1) 



Var(JV, 



-u,0 



Nt 



fug,— 'U,0 



^A) 



< 



C(iVe- A ) 2 ((Z T - + Z^> T ))(i - T-)/a 3 + (Y±°J + Y^^)) ^ CN 2 r,e~ A /s, (9.2) 



by the definition of G A - Now note that under P nbab , Nf u = N^ u '° and N^ u = jv^'^ 



almost surely. Equations (9.1) and (9.2), together with Lemma 7.5 and the conditional Chebychev 
inequality now yield the lemma. □ 



For a barrier function /, let P° f denote the law of BBM with drift — /i s defined in (5.5) starting 



./' 



from a single particle at (without absorption) . Let Rt be the number of particles hitting a before 
the time t and define Zt, Yt by summing wz, respectively wy over the particles at time t. 

Lemma 9.9. Let t ^ and f be a barrier function. Then, 



E° f [Z t ] s; ae*?* - "*, E° f [Y t ] sc e*?* - " and E° f [R t ] sc e^^. 
Proof. By Lemma |3.2| and Girsanov's theorem, we have 



EjM 



e /3 mt e -/xa^y0 



^l(^o) 



«S e 



Pomt^fia^O 



^(X t -f(t/a?)) 



=% e2a : 



t—fia 



which implies the first two inequalities, since Zt ^ aYt- As for the third one, by Lemma 3.1 and 
again Girsanov's theorem, 



E°[R t ] = W°, 



BomHa - 



l (H a ^t) 



«S e 



-^w 



e^ Ha l 



(H a !it) 



=5 e2a : 



rt—^ia 



Furthermore, we have E^fi^] ^ E°[i?i] as in Lemma 5.7 This finishes the proof of the lemma. □ 



9.3 The probability of G\ 

For an interval I, we define the random variable Lj = £]( u s )e^ ^-(sei) counting the number of 
particles hitting the origin during the time interval / and set L n = Lj n . 

Lemma 9.10. Let n ^ 0, s =% 46.2il > I — [ti,t r ] c [0, s] and x e [0, a]. Furthermore, let f be a 
barrier function with Err(/, t r ) ^ C '. Then, for large A and a, we have 

™ rr .„ jCe^a- 3 (t r -t l )(Aw Y (x)+w z (x)), if t t > a 2 or x > a/2 

i Ce tm {wY{x) + a 6 {t r ~ti)wz{x)), for any ti and x . 
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(0) 



Proof. Write E[-] = E[-|T > s]. Define L} = £(«,«)&?*. l(H a (x u )> s6 r>- As in the proof of 



Lemma 5.1 we have for every interval /, x e [0, a] and t ^ s, 

E x f [Lf } ] < e °( E "(/- 4 ))E x [4 0) ] ^ Ce^ x I a {a - x,I), 



(9.3) 



by (2.10) and Corollary 6.7 Lemma 2.1 and (|9.3[) now give, 



(x„(«),«wo)-| 



E^[L / ] = E^4 0) ]+E^ J] E / 

(«, S )e^ ( r 1+) 

^ Ce"" («*(*)/> - *, J) + ^E^zg+J^]) . (9.4) 
now gives for large a, I a {a - x,I) ^ C (a~ ci (t r - £;) sin(7ra;/a) + l( tl<a 2 and x<a / 2 )) , and 



Lemma 



2.1 



by (5.38), Lemmas 5.7 and 6.8, Corollary |6.7| and ( |6.15 ), we have 

& f [z£+l u] ] sj CAw Y (x) (l + ^sin(7rx/a; 



(9.5) 



The lemma now follows from (9.4) and (9.5), together with the hypothesis t r ^ s < 46.211 ^ Ca /A. 



D 

The following lemma is crucial. It will permit to estimate the number of particles turning blue 
upon hitting the origin. 

Lemma 9.11. Suppose that uq is deterministic with Zq > (1 — <5/4)e . Let s ^ fe.2il and I = 
[ti,t r ] c [Ka 2 , s]. Write P = P(- \T > s). Then, for large A and a, 



E 



[ S 



(u,t)eJ? HQ 



(teI,N^ 0) <N) 



^C (5 e- A (a- i + e- A yo)E[L / ]. 



Proof. For an individual u e U and i ^ 0, define iV t = ^, , „(o) l(„£ Uo ), where uo is the 

ancestor of u at the time 0. By the trivial inequality iV t ^ iV 4 for every u and by the 

independence of the initial particles, 



" (u,t)GJSf Ho 



(teI,N { t 0) <N) 



< £ [ E l( te /)P(^ (0) -" < TV) 
(«,t)ei?H 



The lemma now follows from Lemma 9.5, since for every it e ^V(0), we have Zo — tu^(X„(0)) > 
Z - C ^ (1 - <5/2)e A for large A, by hypothesis. □ 



Before the breakout. For every n > 0, define the (sub-probability) measure P n = P(-, T > 
t n+ i). The following lemma gives an estimate on the number of white particles (not counting the 
grey ones) turning blue during the interval I n . 

Lemma 9.12. For every n > 1 with t n +i ^ h§.2i\ , we have for large A and a, 

e fia 

~E n [B n l Gzn ] ^Cs—y- 
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if = t - (K + l)a 2 



Figure 5: The key step in the proof of Lemma 9.12 Conditioned on J^V , the descendants of u are independent 
of the descendants of the other particles by the strong branching property. Furthermore, conditioned on 
&v, the probability that at time t there are less than N particles not descending from u is bounded by 
Cse~ A (a~ 1 + e~ A Yp ''^ u ). These facts imply that the expected proportion of white particles turning blue 
upon hitting a is bounded by C$e~ /a. 



Proof. Let B n be the number of white particles turning blue during I n and which have not hit a 



before the time if = t n — {K + l)a 2 = £ n _i + 2a 2 and let B^ 
now yield, 



B n — Bn . Lemmas 



9.10 



and 



E„[BWl G _ \&t,]^C s e 



-A 



1 



-A 






^C 5 e^Y^ + e- A {Y^f)l Gztn 
Proposition |5.2| now gives 

B n [B^l Gz J ^ Cse^(a- 2 + e- 2A a-^[Y^_ i+a2 l Gztn J) < C s e^/a 2 . 
As for the remaining particles, by the strong branching property, 

B n [B n 1+ h Gzn | J? HaAt ,] = J] [ E^)[L trir+dp) ]P n (7V r W < N, G z , n | Jf«) 

n_^(01 Jin 



(u,s)e; 



(0), 



^C s (e- A /a)R^'(e^/a)Q a [Z] 



by Corollary 9.6 and Lemma 9.10 Lemma 5.6 and ( 6.46| ) then give 

B n [B n 1+ h Gz J ^ C 5 (e» a /a 2 )At n+1 /a 3 ^ Cge^/a 2 , 



by the hypothesis on t n+ \. The lemma now follows from (|9.6|) and (9.7) 



9.11 



(9.6) 



(9.7) 

D 



Up to now, we have only considered the white particles turning blue and will now turn to B n ot , 
n > 0. For n > and < k < n — 2, we define Bk >n to be the number of particles that turn blue at 
a time t 6 I n , have an ancestor that turned blue at a time t' e Ij., have none that has hit between 
tk+2 and t n , have not hit a between if and tk+2 and did not break out between tk+2 and t. 

Let Gb,n the event that no blue particle hits a before t n and G g ^ n the event that no particle that 
turned from white to blue before t n -\ has a grey descendant that breaks out before ^.2il) - Then set 

G™ = G z ,n n G Wl n G g , n n {5 ^ eV a /a 2 , F < v}- 
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Lemma 9.13. For every n > with t n+ \ < fa.2il fl^rf k ^ n — 2, we have for large A and a, 



E n [-B fcn l G tot] < Cs~E k [B k l G tot] 

Proof. Let £$k be the stopping line consisting of the particles that turn blue during If., at the 
moment at which they turn blue (hence, -E>Jjl ot = j^SS^). Note that S$k, a jV+ 

,(0) 



(0) 



fc - -*s 



such that the descendants of 3$^ and of ^ _ are independent, given their past, by the strong 
branching property. Note also that Gj(, ot e & (o) , since by definition this a-field contains all 

*n— 1 



the information about the descendants of the particles in B§j, j < k. By Corollary 9.6, we then 
have 



En[-Bfc,nlGt ot I &t 



,aJC 



(0) ] 



{u,s)c3ti k 



S 



L(H a {X v )>t k+2 ,X v (t k+2 )E(0,a),tEl n ) 



■ (v,t)eJf ak+2 



(9.1 



9.9 



and 



9.10 



each summand on 



where here H^ +2 (X) = inf{£ 5= ife+2 : X t = 0}. By Lemmas 

the right-hand side of the above inequality is bounded by Cge t ^W} l [Y tk+2 ^ s ] ^ C$. The lemma 

follows. D 

Lemma 9.14. For all n ^ 1 with t n+ \ ^ fe.2il , we have E n LB^ ot l G tot] ^ Cse^/a 2 for large A and 
a. 



Proof. By Lemma 9.12 the statement is true for n = 1, because B\ ot = B\ by definition. Now we 

"and the fact that B^ = B < e A e^ a /a 2 on G%\ 



have for every n, by Lemmas 
n-2 



9.12 



and 



9.13 



-A ™-2 



E n LB^ ot l G tot] = N^ E n [i? fcin l G tot] + E n LB n l G tot] < CM 2j Efet-^fc iG* *] 

fc=o a fe=l 

The lemma now follows easily by induction over n, since n ^ a by hypothesis. 



3 /ia 



(9.9) 



a 



Lemma 9.15. Suppose (HBq). For large A and a, we have P n ((G^ ot ) c ) ^ C$e 2A ^ 3 n/a + e 2 for 
every n > 0. 

Proof. Let n > 1. By definition, the event G| 1 „ 1 \Gb,n implies that a descendant of a particle in 
SSn-i hits a before £ n . Markov's inequality and Lemma 9.9 then imply, 



P n (G^ 1 \G'b,n) ^ E n [£C 2 l G totJ x supE"^] < CsaT\ 



t^t 2 



(9.10) 



by Lemma |9. 14 Furthermore, by Proposition |6.2| and Lemmas |9.9| and |9.14 we have 



n {{GT-i n G b , n )\G &n ) sj E n [S^ 2 l G totJ x supp^fe^ + Y t ] sj C 5 e 



-2A/3 a -l 



t^t 2 



(9.11) 



by (6.14) and (6.2). By ( |9.10 ) and (9.11), we have for large A and o, 

P„(Gz,n\0 < P«(Gz,n\Cl) + Qe-^^a" 1 < P n _!(G Z ,n-l\G^l) + C^^oTK 
By induction over n and the hypothesis, this gives for every n, 

Pn(G z ,n\Gn 0t ) < Cge-^^n/a) + e 2 . 



Together with Lemma 9.7, the statement follows. 



D 
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Define the random variable no := L-^AiJ ana - n °te that t no -\ > T for large A. Define the 
events G nw = {Zf^ < e A '\ Y* < e A ' 2 /a, B^_ x < ^/a} and G™ = U n (GT n {n = n}). 

Lemma 9.16. For large A and a, P((G*£ n G^)\G nw ) s= C 5 e- A l 2 . 



Proof. By definition, #^ a &%„ . Lemmas 



»o" 



6.9 


9.9 



and 



9.14 



then give 



E 



00 

[Z£ o l G tot nG J ^ CsupE°[Z t + Ay t ]( J] E[E n [S^lGto. | ^sfejle* 



ts=<2 



n=0 



< 



[Via 3 /*lJ 
n=0 



(9.12) 



Similarly, we have 

P( J Rf; ji > 0, G™, G v ) ^ Csv, and Ep^l^ _^ 0) l G tot nG J ^ G 5/ 



(9.13) 



Finally, we have as in (9.12), by Lemmas 9.13 and 9.14 

no -2 
E [ B fo-l 1 G^- 1 nC«] = E [ 2 E [- B *'.»>0-1 1 G^_ 1 I ^Hq] 1 ^ 



fc=0 



|V^ 3 /*iJ 

^C 5 (e- A /a) 2 Efc^rictot] < C s e- A+ ^/a 2 . (9.14) 

fc=0 



The lemma now follows from (9.12), (9.13) and (9.14), together with Markov's inequality 



□ 



After the breakout. We study now the system after the breakout. Recall that uq = \T/t\\ 
by definition and define ri\ := [0i/ti + lj. We define ^ A = ^a n ^jyv A ag n _ x and the event 

G A = Ga n G nw £ ^j\- We denote by the superscript "gr<" the quantities relative to particles the 
descending from those that were grey before or at time t no . Then let G^ bab be the intersection of 
Gnbab with the event that none of these particles hits a between t no and @i + e a? before hitting 
0. Define then the (sub-probability) measure P^ = P(-, G nhab ). 

Lemma 9.17. For large A and a, P(G A ) Ss 1 - Ce 5 / 4 and P((Gj ibab ) c | ^ A )l G a ^ Ce 2 . 



Proof. On Ga, we have t no < ^g 21I for large A and a, in particular, no ^ a. By Lemmas 9.15 and 



9.16 and (6.2), we then have for large A and a, 

V(G A \G A ) < P(G A \G t „° t ) + P((G A n G^)\G nw ) sj Ce 2 . 



With Lemma 7.4 this yields the first statement. For the second statement, we have by Lemmas 5.6 

s\ Tl C\ 1 1 

P«" n x a 21 > I J1)1hi < C{{e A /a)Zf + Yf )l ni sj Ce A l\ 
Together with Lemma |7,5[ this implies the lemma. 



□ 

We change the notation of B n a bit: it is defined to be the number of particles hitting during 
I n for the first time after t no -\. In particular, we also count the descendants of the grey particles 
at that time. 

Lemma 9.18. For large A and a and n ^ no — 1, we have E^ bab LB n | ^ A ]l r $ < Cse~ 5A ' 3 e tia /a. 
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Proof. Define Jzf := .jVj-- u ^v^' 7 "). We have for t e J n , as in Lemma 9.11 



ELb^n + ^l-^i]^ <C 2 E^«>[L B ]P(JVT u + iV?*" ni <iV|^A)l 



G\ 



sc 



(u,t)6jSf 

(u,t)eif 



G» 



by Lemma 9.8. By Lemma 
a, 



9.10 



(6.1) and the definition of G and Gw, this gives for large A and 



E 



nbabL 



[B n + B** | J^]l„„ sj C (S e^r ? (ee A a)- 1 (Z T - + AY T _ + Z^> T ))1 



fi 3 . 



(9.15) 



sj C s e^ a rje- 2 /a. 
Further more , by the independence of the check- and bar particles from the others, we have by 



9.8 



Lemmas 



large enough, 



and 9.10 and by the inequality wy(x) ^ a 1 e ' a x ')' 2 , valid for every x ^ a/2 and a 



K b ,dBn + B n \^ A ]l 



< C 5e ^r ? ( £ e^)- 1 a- 1 (%vT- + AY 0vT - + AS* + ZF + AYF)\^ < Qe^Va, (9.16) 



by (6.1) and the definition of G\. Similarly, we have for n > no, 



nbabL 



(9.17) 



and B^ _ l < rje^/a on the event G A . The lemma now follows from (9.15), (9.16) and (9.17), 



together with (6.3) and (6.2). 



□ 



For n ^ no — 1, define G^ ot to be the event that no descendant of a particle which has been 
coloured blue between £ no _i and t n hits a before 0i + e A a 2 . 

Lemma 9.19. We have for every n ^ no and for large A and a, 



E 



nbab^n 0t l^l^I]l G i^^^ 5A/3! 



/ia 



Proof. The proof is similar to the proof of Lemma 9.14 Substituting Lemma 9.5 by Lemma 9.8 
one first shows that for every no — 1 ^ k ^ n — 2, one has 



E 



nbab 



[B Kn i &ot | ^ji G » ^ Csve- A K^d B r Igto. I ^£]1 G » , 



which, by a recurrence similar to (9.9), yields the lemma 



(9.18) 

□ 



Lemma 9.20. For large A and a, P nbab ((G t „° 1 t ) c | &l)l G » ^ 1 



a. 



Proof. Similarly to the proof of Lemma 9.15, but using Lemmas 5.6 and 5.7 instead of Proposi- 
tion RT21 we have for large A and a, 



nbab*. 



{d^AG™ I ^i)l G » < C s Bl hah [B^ 2 l &oK | ^j> ]1_ B e A -^/a < C s a~ 2 , 



where the last inequality follows from Lemma 9.19 and (6.3). The lemma now follows by induction 
over n. □ 
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Lemma 9.21. Suppose (HEr Q ). Then ~P(G\) > 1 — e 4//3 for large A and a. 
Proof. By Lemma |9.5[ Corollary |9.6| and the union bound we have 

P(3n sc n : N(t n ) < N, G z ,n , G<*) sj C & e~ A . 



Furthermore, by Lemma 9.8, we have 

P(3n < n < m : N(t n ) < N or JV(9i) < N, G A ) < C s r)/e sj C 5 e- A , 
by (6.3) and (|6.2[). Now, by Proposition 5.2 and Lemmas 9.9 and|9.17 



(9.19) 



(9.20) 



ni- 1 



"l A \ u n / 

n=no — 1 

^Ce^ + C^Ce^ 2 , 



G\ 



(9.21) 



by Lemma 9.19 Similarly, we get 
Moreover, we have by Lemma |9. 19 



(9.22) 



Setting X = P(sg 1)ei+tl] ^ e-^e^/a & @1 ) , we then have E[(l-X)l G ^ ng ^ nG ^] sj C^- 2 ^ 3 
by Markov's inequality. Applying the Markov inequality once more to X as in the proof of 



Lemma 8.1 yields 



P(X < 1 - e 2 , G* bab n G£ n G» A ) < C^V 2 ^ 3 < C 5 e 2 , 



(9.23) 



by (6.2). The lemma now follows from (9.19), (9.20) and ( |9.23 ), Lemmas 9.17 and 9.20 and Markov's 
inequality applied to (9.21) and (9.22). □ 



9.4 Proofs of the main results 



Proof of Lemma 9.1, The property about Zq and Yq follows as in the proof of Lemma 8.8 from 



21). Moreover, again as in the proof of Lemma 8.8 we have by (8.21), (8.23) and (8.25), for 
every i ^ \e N], 

P(iV t 7 < N, R tl = 0) < C 5 e- A /a and P{R tl > 0) < C 5 e A /a, 

where we denote the initial particles by 1, ... , \e °N~\. By the independence of the particles, we now 
have, 

\e s N] 

E[B l {Rti=0) ] ^ ^ V(Nt7<N,R tl =0)E Xiio) [Lol {RH=0) ] 



i=l 



^ C 5 (e- A /a)e^E[y ] < Cge^/a 



fj.a/2 



by J83i| >. Setting X = P(B < e _A e^ /a J^o), we then have E[l - JTJ < C,se A /a by Markov's 
inequality. Applying Markov's inequality once more to X as in the proof of Lemma |8.1| yields the 
lemma. □ 
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Proof of Proposition 9.2. Exactly the same reasoning as in the proof of Proposition 8.2 but using 
Lemma 19.211 instead of Lemma 18.141 □ 

O. Define p = x ae -2s, such that (1 + pp)e~ fip ^ ae~ 2S (l + o(l)). As in the 



Proof of Proposition 
proof of Proposition 



3] (see (8.57)), it is enough to show that 

sup sup P u (N$°\p) >aN,T>t)^ C 5 e 



5/4 



(9.24) 



for large A and a. Let v e G and t ^ e a 2 . If t > y^a 3 , then the above probability is bounded 
by Ce 2 by Proposition 6.2 Suppose therefore that e A a 2 ^ t ^ \feo? and write P = P(- 1 T > t). 



By definition, we have (1 + x a )e~ Xa = a. Let n be the largest integer such that t n+ \ < t; note that 
n > 1 for large A. By Lemmas |6.9| and 6.12 together with Chebychev's and Markov's inequalities, 
(6.2) and (6.1), we have for large A and a, 

P(\e- A Z tn - 1| < e 1 '\ Y tn sc e 4A / 3 /a) > 1 - C s e^. (9.25) 



And as in the proof of Lemma 9.16, the same holds for Zf r and Yf r as well, which yields 

P(|e-X g " II < ^ C < e 4A / 3 /a) > 1 - C^l\ 
Lemmas 



(9.26) 



5.6 



and 



9.9 



and Corollary 6.7 then show that P(-Rf°* f] > 0) < o(l). 



Let iV t (p) be the white and grey particles to the right of p at time t which descend from the 



white and grey particles at time t n and which have not hit a between t n and t. By Proposition 5.3 
and Corollary |6.7[ we have for large A and a, 

(9.27) 
(9.28) 



mr S (p)\^ n ]^e- s / 2 a Ne- A Z^, 
Var(Af«(p) | J^J < C 5 , a a 2 N 2 e- 2A {a- l Z^ + Y?). 
Chebychev's inequality and ( 9.26| ) then give for large A and a, 

P(iV t wg (p) > e^aN) sc Ce 5 / 4 . 



9.11 



(9.29) 

iave 
and 



Furthermore, denote by N^ (p) the number of blue particles to the right of p at time t which have 
turned blue after t n + a 2 and which have not hit a between t n and t. Then by Lemma 
(9.3), we have 

B[Nt luc (p)\^t n ]^C s e- A (a- 1 + 



sj C & e- A (a- 1 



e- A Y. w 



YZ S ) 2 I" E[L [r , r+dT] (n, S )]E°[iV t _ T (p)] 

[ I a (a-X u (s),dT)E°[N t _ tn _ T (p)]. 

Ja 2 

(9.30) 



\ ' i Z r , 

s « 

(n,s)G^ wg 



fj,X u (s) 



Set d = t — t n . Note that [K + 3)a 2 < d < 2(if + 3)a 2 . By Lemma 3.2 and Girsanov's theorem, we 
have now for every r ^ 0, 

L (X T 



ir 2 it 2 p° 

E°LV T (p)] = e^ mr W\{X T >p)= e^ T W°[e-^l {XT>p) ] = e^ T <T» z g d - r {z) dz, (9.31) 

Jp 



where g T (x) = (2ttt) 1 ' 2 e x '^ 2t ' is the Gaussian density with variance r. If r ^ a 2 , then 
sup^ g T (z) ^ C/a, such that for every x e [0, a] and z > p, 



J-d— a 
J a (a - x, dr) 5d _ T (2) < ^"^"(a - x, [a 2 , d - a 2 ]) sj CKaT x sin(7rx/a) 
a 2 



(9.32) 
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t ^ a 2 , such that 



by Lemma 2.1 Moreover, by Lemma 2.1 we have I a (a — x,dr) < Ca sin(7rx/a) dr for every 



rd ra* 

I I a (a - x, dT)g d - T (z) < CoT 2 sin(7rx/a) r~ 1/2 dr sj Ca -1 sin(7rx/a). 
Jd-a 2 Jo 



(9.33) 



Equations (|9_30|), (|9_3lJ), (|9_32j) and fl9.33[ ) now yield 

E[iV t bluc (p) | J^J < Ce-V^a-^a- 1 + e"^)^ < C 5 ae- 2A iV(l + oe"^)^. 

(9.34) 

Furthermore, if Nl° st (p) denotes the particles to the right of p at time t descending from those 
turning blue between t n -\ and t n + a 2 , then by ( 9.31[ ) and the supremum bound on g T (z), 



E[Nl esi (p)Gi ot ] < Csa-h-WEiiB? + BZUlc**] < Cfeae"^, 



(9.35) 



by Lemma 9.14 Markov's inequality applied to (9.35) and (9.34) together with (9.26), Lemmas 9.15 
and |9.7| give 



P(N t 



rest 



A^ 



blue 



> eaN) sc C 5 e 5 



(9.36) 



by (6.2). The statement now follows from (9.29), ( |9.36 ) and the above-mentioned bound on 



Proof of Le mma 9.4 ■ By definition, the B"-BBM and C"-BBM coincide until Q n on the set G n - By 
Proposition 9.2 we have P(Gn) ^ 1 — ne 1+b for some numerical constant b > and by the coupling 



of (O n ) n ;>o with a Poisson process o f int ensi ty p Be A ir ~ e \ we have P(0 e _i_&/ 2 > e fe ' 4 ) — > 1, as 
yl and a go to infinity. Propositions 9.2 and 9.3 then readily transfer to the C"-BBM. □ 



10 Coupling the iV-BBM with other braching Brownian motions 

We will first establish a monotone coupling between the iV-BBM and a class of slightly more ge neral 
BBM with selection which includes the B -BBM from Section [sj and the C"-BBM from Lemma 19. 4 
In a second part, Theorem |1.1| is proven. 

10.1 A monotone coupling between iV-BBM and more general particle systems 

A selection mechanism for branching Brownian motion is by definition a stopping line 5£ , which 
has the interpretation that if (u, t) e „5f , we think of u being killed at the time t. The set of particles 
in the system at time t then consists of all the particles u £ ^{t) which do not have an ancestor 
which has been killed at a time s ^ t, i.e. all the particles u 6 ^(t) with Jzf ^ (u, t). 

Now suppose we have two systems of BBM with selection, the A^ + -BBM and the A^-BBM, 
whose selection mechanisms satisfy the following rules. 

1. Only left-most particles are killed, i.e. if (u, t) s Jzf, then X u (t) ^ X v {t) for all v s ^V{t) with 

2. Af + -BBM: Whenever a particle gets killed, there are at least A" particles to its right (but 
not necessarily all the particles which have A" particles to their right get killed). Af _ -BBM: 
Whenever at least N particles are to the right of a particle, it gets killed (but possibly more 
particles get killed). 

Let u^, v^ and v^ be the empirical measures of the particles at time t in A^ + -BBM, Af _ -BBM 
and A^-BBM, respectively. On the space of finite counting measures on R we denote by < the 
usual stochastic ordering: For two counting measures v\ and U2, we write v\ < v<i if and only if 
ui([x, oo)) < V2(\x, oo)) for every x e R. If x±, . . . , x n and yi, ■ ■ ■ ,y m denote the atoms of v\ and V2 
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respectively, then this is equivalent to the existence of an injective map : {1, . . . , n} — » {1, . . . , m] 
with Xj ^ y<t>{i) f° r all i e {1, . . . , n}. Furthermore, for two families of counting measures {y\(t))t^o 
and (i*2(i))t>o> we write (^i(t))t^o ^ (^2(i))t>o if ^i(i) ^ ^2 CO for every i > 0. If (i/i(t)) t>0 and 

St 

(^2(*))t;s0 are random, then we write (^i(i))t^o ^ (&*2(i))t^o if there exists a coupling between the 
two (i.e. a realisation of both on the same probability space), such that (^i(t))t^o ^ (^2(*))t>o- 

Lemma 10.1. Suppose that Vq < v$ < v$ . Then (f t ~jt>o ^ (^t )t>o ^ (^ Jt^o- 

Proof. We only prove the second inequality z//^ < f t + , the proof of the first one is similar. By a 
coupling argument and conditioning on &q it is enough to show it for deterministic Vq and v^. 

Let n+ = f(j"(R) and let n = ( n(1) > n(2) > • • • ,^ {n+) ) be a forest of independent BBM trees with 
the atoms of Vq as initial positions. We denote by ^V u (t) the set of individuals alive 7 ] at time t 
and by X^{t) the position of an individual u 6 ^V n {t). Denote by jV + {t) a ^Y u (t) the subset 
of individuals which form the iV + -BBM (i.e. those which have not been killed by the selection 
mechanism of the iV + -BBM). We set f 4 + = Yjue,yV+(t) ^x n (t)- 

From the forest n we will construct a family of forests I Sj = (Sj , . . . , E t ) I (not necessarily 

comprised of independent BBM trees), such that 

• if t\ ^ £2, then the forests 3^ and Et 2 agree on the time interval [0,£i], 

• the initial positions in the forest Ho are the atoms of Vq , 

• for every t ^ 0, the iV-BBM is embedded in St up to the time t, i.e. for < s < t, if 
^ a (s) denotes the set of individuals^! from Ej alive at time s and Xjf(s) the position of 
the individual u e jV~(s), then there is a subset ^V N {s) c jV~(s) such that (^)o^s^t = 
(2 u6i/ fiv( s ) <^x H (s))o<s<t * s e< l u al in law to the empirical measure of iV-BBM run until time t, 

• for every t ^ 0, there exists a (random) injective map 4>t : Jf^ (t) — * ^V + (t), such that 
**(*)< Xj (u) (t) for every ue^(t). 

We will say that the individuals u and (f>t{u) are connected. If at a time £ an individual v e r yV + (t) 
is not connected to another individual (i.e. v $ (f)t(^V N (t))), we say that v is /ree. 

The construction of the coupling goes as follows: Since Vq < z/q" , we can construct Ho, ^V (0) 
and 0o, such that for every u e ^^(0) we have X£(0) ^ Xp, ..(0) and the subtrees Eq and 
tjOpoWiO) are the sarne U p t translation. We now define a sequence of random times (t n ) n ^o 
recursively by to = and for each n, we define t n+ \ to be the first time after t n at which either 

1. a particle of the iV-BBM branches, or 

2. the left-most particle of the A r+ -BBM dies without a particle of the iV-BBM branching. 
Note that between these times, there is nothing to do, i.e. if Hj and <pt are defined for all t < t n , 
then we simply set Hj = Ej n and 0t = 0£ n for all t e [£ n ,£ n +i). This is clear if no branching 
occurs between the times t n and t n +i, but even if a particle of the iV + -BBM branches at a time 
t E (i n >£n+i); since no branching occurs in the iV-BBM between t n and t n+ \, this particle must be 
free at the time t and therefore no change in Ej and 0j are necessary. At the times t n however, we 
may need to "rewire" particles. We show that this is always possible in such a way that the above 
conditions on the family St are satisfied. Figure [6] is a graphical description of the key step. 

Let n e N and suppose that 

a) Si and cj>t have been defined for all t < t n+ \ and are equal to 5 <n and 4>t n , respectively, 

b) for each u £ jV (t n ), the subtrees S 4 ' and U t tn(u> ' n > are the same, up to translation. 
Note that this is the case for n = 0. We now distinguish between the two cases above, starting with 
the second: 



7 The term "alive" has the same meaning here as in Section 
8 Note that this does not depend on t. 



3.1 
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n-bbm qn o o o 

iV+-BBM fee) O QDO 

w' v' 

Figure 6: The connection between the particles w of the 7V-BBM and w' of the 7V+-BBM breaks. By 
definition of the 7V+-BBM, there exists a free particle v' to the right of w' and w is rewired to that particle. 



Case 2: The left-most particle w' of the iV + -BBM gets killed without a particle of the iV-BBM 
branching. If w 1 is free, nothing has to be done. Suppose therefore that w' is connected to a particle 
w of the iV-BBM. Then, since there are at most TV — 1 remaining particles in the iV-BBM and there 
are at least N particles to the right of w' in the iV + -BBM (otherwise it would not have been killed), 
at least one of those particles is free. Denote this particle by v' . We then "rewire" the particle w to 
v' by setting (j) tn+1 {w) := v' and define H 4n+1 by replacing the subtree H^ n in E tn by IL; ' ™ , 
properly translated. Note that we then have 

^4>t n+1 (w) (Wl) = ^»'(Wr) > X w ,(t n+ i-) Js X~(tn+i—), 

where the first inequality follows from the fact that w' is the left-most individual in iV + -BBM at 
time t n+ i and the second inequality holds by hypothesis. 

If more than one particle of the A r+ -BBM gets killed at the time t n +i, we repeat the above 
procedure for every particle, starting from the left-most. 

Case 1: A particle u of the iV-BBM branches at time t n +i- By the hypothesis b), the particle 
4>t n {u) then branches as well into the same number of children. We then define 4>t n (uk) = 4>t n (u)k 
for each fee {1, . . . ,C U } (recall that C u denotes the number of children of u), i.e. we connect each 
child of u to the corresponding child of 4>t n {u). Now first define 4 > t n+1 to be the restriction of 4>t n to 
the surviving particles, i.e. to JV (t n +i)- Then continue as in Case 2, i.e. for each particle w' of the 
iV + -BBM which gets killed and which is connected through <j)' t 1 to a particle w of the iV-BBM, 
rewire w to a free particle v'. In the end, we get 4>t n+1 - 

In both cases, we then set E t = S in+1 and <pt = <t>t n+ i for each t 6 [t n+ i,t n+ 2), as explained 
above. Note that each time we are rewiring a particle, we rewire it to a particle whose subtree is 
independent of the others by the strong branching property, whence the particles from JV (i) and 
^V + {t) still follow the law of iV-BBM and iV + -BBM, respectively. Furthermore, we have for every 
oj e Q,: vf (u) < f t + (w) for every t ^ 0. This finishes the proof. □ 

10.2 Proof of Theorem ITTT1 

Let (vt)feo be the empirical measure process of iV-BBM starting from N particles independently 
distributed according to the density proportional to sin(7rx/aAr)e _:E l( ;ce (o )ajv )), where a^ = (log iV + 
3 log log N). We wish to show that for every a e (0, 1), the finite-dimensional distributions of the 
process (M^(tlog N)) t ^>o (defined in the introduction) converge weakly as iV —*■ oo to those of the 
Levy process (Lj)tj>o stated in Theorem |1 . 1| starting from x a . We will do this by proving seperately 
a lower and an upper bound and show that in the limit these bounds coincide and equal the Levy 
process (L t )t^o- 

Lower bound. Fix a 6 (0, 1) and 5 > 0. We will let iV and in parallel A and a go to infinity (in 
the meaning of Section [71) in such a way that iV = 2-7re + a _3 e Ata and such that A goes to infinity 
sufficiently slowly such that the results from Section ^1 hold. We then have with c = log(27r), 

a = a N - (A + 5 + c + o(l)), and a = u N =(A + 6 + c + oil)), (10.1) 

a 6 
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where n is defined in ( |5,1[ ). Let {v\)t^o be the empirical measure process of B -BBM starting from 
the initial configuration (HBj_), i.e. Vq is obtained from [e ATJ particles distributed independently 
according to the density proportional to sm.{'Kx/a)e~i JjX \tQ a \{x). An easy calculation now shows 

St 



that z/q < Vq for large N, by (10.1). Now, if X\ denotes the barrier process of the B^-BBM, then 
(f 4 ~)t^o = { u t + fit + X^)t^>o is by definition an instance of the Af~-BBM, as defined earlier in this 
section. Here , for a measure v and a number x we denote by v + x the measure v translated by x. 

now gives for large N, (v~)t^o < (^t )t>0, which by definition implies 



Lemma 



10.1 



Given ^ t\ < . 
We then have 



St 



(qu>f ))^o > (qu^(^))oo- 



< t n , we now define t. 



N 



,-3 



U log 3 N, such that t 



N 



(10.2) 
tt for every i, as A" — > oo. 



rJV 



M^lUlog'N 



St 



i=l, 



> u ^*Wtv) 



/ijvti log 3 A^; 



i=i, 



by (10.2) 



M n a(»aH N ' 



X a 3 tf 



TT 2 (A + 5 + C + 0(l))tf) 



=l,...,n 



by (10.1) 



st 



L, 



0(5) 



'i=l, 



for large A/", 



where the last inequality follows from Propositions |8 . 2| and 8.3, with (Lfjt^o being the Levy process 
from the statement of Theorem 1.1 starting from 0. Letting first A^ — > oo, then S — > yields the 
proof of the lower bound. 



Upper bound. The proof is analogous to the previous case, relying on Propositions 9.2 and 

and 8.3. There are only two differences to notice: First, the B"- 



9.3 



instead of Propositions 



1.2 



BBM is not a realisation of the A^ + -BBM. Howe ver, the C"-BBM (defined in Lemma 9.4) is such a 



realisation and by that lemma, Propositions 9.2 and 



9.3 



hold for the C^-BBM as well. Second, if A 



st 



is distributed according to (HBjJ, we do not have z^q < v$. However, if v^ is obtained from v { 




by killing the particles in the interval [a^ — A 2 , a??], then by (10.1), a quick calculation shows that 

st 






y 



with high probability and z?q < vt as A" goes to infinity. This finishes the proof of the 



upper bound and of Theorem |1.1| 
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